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ABSTRACT 

To probe the distribution and physical characteristics of interstellar gas at temperatures T « 3 x 10^ K 
in the disk of the MiUcy Way, we have used the Far Ultraviolet Spectroscopic Explorer (FUSE) to observe 
absorption lines of O VI A1032 toward 148 early-type stars situated at distances > 1 kpc. After subtracting 
off a mild excess of O VI arising from the Local Bubble, combining our new results with earher surveys 
of O VI, and eliminating stars that show conspicuous locaUzed X-ray emission, we find an average O VI 
mid-plane density no = 1.3 x 10"^ cm"^. The density decreases away from the plane of the Galaxy in 
a way that is consistent with an exponential scale height of 3.2 kpc at negative latitudes or 4.6 kpc at 
positive latitudes. Average volume densities of O VI along different sight Unes exhibit a dispersion of 
about 0.26 dex, irrespective of the distances to the target stars. This indicates that O VI does not arise 
in randomly situated clouds of a fixed size and density, but instead is distributed in regions that have 
a very broad range of column densities, with the more strongly absorbing clouds having a lower space 
density. Line widths and centroid velocities are much larger than those expected from differential Galactic 
rotation, but they are nevertheless correlated with distance and A'CO VI), which reinforces our picture 
of a diverse population of hot plasma regions that are ubiquitous over the entire Galactic disk. The 
velocity extremes of the O VI profiles show a loose correlation with those of very strong lines of less 
ionized species, supporting a picture of a turbulent, multiphase medium churned by shock-heated gas 
from multiple supernova explosions. 

Subject headings: Galaxy:disk — ISM: clouds — ISM:kinematics and dynamics — ISM:structure — ultravioletrlSM 
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1. INTRODUCTION 

Very hot plasmas are recognized as important con- 
stituents in a broad range of astrophysical contexts, 
ranging from the solar and stellar transition layers and 
coronae up to the largest arrangements of matter in 
the universe that form web-like structures that bridge 
groups and clusters of galaxies. In all cases, the pro- 
cesses responsible for heating the gas give important 
insights on evolutionary and dynamical phenomena in 
the systems under study. On the largest scales, there 
is a recognition that hot gas in the intracluster me- 
dia and even the intergalactic media at low redshifts 
probably dominate the baryonic content of the universe 



(e.g.|Cen & Ostriker| 1999} [Tripp et al.|2000} |Danforth 



|& Shull'' 2005[|Bregman||2007| l. In this paper, we fo 
cus on the presence of hot gas within the interstellar 
medium (ISM) of the Milky Way, a system intermedi- 
ate in size between the extremes just mentioned. This 
is a topic that has been studied over the last 3 decades. 



and shortly (§ 1 .2 1 we will look back upon some of the 



key developments in this field to help define the moti- 
vations for our present study. 

1.1. Basic Principles of Detecting Hot Gas 

Atoms within plasmas at temperatures T > 10^ K 
have many of their outer electrons stripped by colli- 
sions. As a consequence, there is an upward shift in 
the energy of their resonance transitions, pushing im- 
portant spectroscopic information into the ultraviolet 
and X-ray bands. It is for this reason that progress 
in understanding the distribution and properties of hot 
gases in space have had to await the development of 
space-borne instruments that could observe radiation 
at energies above those in the visible band[^ Of all the 
transitions from highly ionized atoms at wavelengths 
longward of the Lyman limit, the AA1032, 1037 dou- 
blet of O VI is the most suitable for study. The ad- 
vantages of observing O VI are threefold; (1) the ion- 
ization potential of O"*"^ is high enough (114eV) to 
insure that there is no appreciable contamination by 
cooler material photoionized by starlight]^ as is the 
case with other partly stripped, Li-like atoms such as 



'Forbidden transitions in tlie visible from highly ionized atoms are 
generally too weak to observe ( Graney & Sarazin 19901, except as 
seen in emission for exceptionally bright sources such as supernova 
remnants jW oodgate et al. 1974 1975, Teske 1990 Sauvageot &| 
|Decourcherie^l995 ■ Sauvageot et al. J999 .^Szentgyorgyi et al. 2000^ ^ 
^Except for small contributions from very hot. He-poor white dwarfs 
jPupree & Raymond|19"83) . 



Si IV and C IV ( |Cowie et al.||1981b| l, (2) the transi- 
tions are strong [log/A = 2.138 for the 1032 A tran- 
sition ( |Morton|[2003| l], and (3) the cosmic abundance 
of oxygen is high. One mildly negative factor is that 
the peak abundance of the five-times ionized form of 
oxygen (22% at T = 3 x 10^ K under collisional ion- 
ization equilibrium conditions) is lower than represen- 
tative fractions that can be reached by other ionization 
stages just below or above that of O VI ([Shapiro & 



Moore 1976 a| [Shull & van Steenbergl|T982| [Suther- 
|land & Dopita 1993i Also, the temperature range over 



which O VI is the dominant species is smaller than that 
of many other ions, such as O VII A 2L6 seen in the 
X-ray band. 

The easiest way to observe the O VI transitions is 
to record them in absorption in the spectrum of a UV- 
bright source. In the disk of the Galaxy, the most suit- 
able stars are those hotter than spectral type B3 that 
have large projected rotation velocities vsin/. Diffi- 
culties with poorly behaved stellar continua sometimes 
arise to create ambiguities in the absorption strengths; 
this issue is one that we address later in §2.6| for our 
current study. For probing the Galactic halo, AGNs 
and quasars make excellent background sources. 

Energetic collisions in a hot gas not only ionize the 
atoms, but they also excite them to electronic states 
at energies well above the ground state. As a conse- 
quence, the atoms emit radiation in their strong reso- 
nance lines. O VI can be seen in emission (see, e.g., 
Table 1 of [Dixon et al.|2006tpte & Dixon|2006| l al- 
though the flux levels are very faint (~ 1 - 4 x 10~'^ 
ergs cm"^ s"' arcsec"^) and require long integration 
times to achieve a satisfactory detection. A principal 
difference between the absorption and emission mea- 
surements is that the former is proportional to the line 
integral of the electron density, while the latter senses 
the square of the electron density [but with some sen- 
sitivity to temperature beyond just the variation from 
the ion fraction cui-ve ( [Shull & Slavin[[T994j )]. This 
difference can be exploited to measure the filling fac- 
tor of the gas, (rtg)^/ (n^). One drawback in measuring 
the O VI emission is that self absorption and resonant 
scattering can complicate the interpretation if the op- 
tical depths are of order or greater than unity ( [Shelton[ 
[et al.|200T) . Also, the emission can be attenuated by 
the absorption caused by dust in the foreground. 

Line emission from atoms other than O VI also ra- 
diate from hot gases. As a rule, the characteristic en- 
ergy bands where the emissions are strongest and most 
plentiful increase toward higher energies as the tem- 
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perature rises. For gas at temperatures of 282,000 K 
where O VI has its maximum ion fraction in collisional 
equilibrium, the characteristic energy of the strongest 
emission-Hne features is below 300eV ( |Kato||T976l l. 
Unfortunately, there are a number of drawbacks that 
limit our ability to make good use of the X-ray back- 
ground measurements as we interpret the O VI ab- 
sorption data. First, at energies below 300 eV the ab- 
sorption by foreground neutral hydrogen and helium 
is strong (Mo rrison & McCammon|1983 i. Second, it 
is often difficult to know how much emission arises 
from beyond the target used for the O VI measure- 
ment. Third, over recent years there has been an in- 
creasing awareness that the soft X-ray background is 
contaminated by emission from charge exchange be- 
tween heavy solar wind ions and the local interstel- 
lar gas ( Cravens|2000[ [Cravens et al.|200l" i, a process 
identical to that which produces X-ray emission from 
comets (Kharchenko & Dalgamo 2000{ [Krasnopols^ 
[et al.|2002i|Kharchenko et al.,2003^ . This contamina- 
tion may be as much as half of the observed diffuse 
background level (Lallement 2004). Finally, with the 
exception of the lowest energy X-ray measurements, 
such as those that operate below the 111 eV absorption 
cutoff of a beryllium filter ( [Bloch et al.|1986p , the ob- 
servations have much stronger responses to gases that 
ai-e too hot (r > 2 x 10* K) to yield much O VI. All 
but the second reason listed above may help to explain 
why [Savage et al.[ ( [2003 ) found a poor correlation be- 
tween soft X-ray emission and the strength of O VI 
absorption toward extragalactic targets at high Galac- 
tic latitudes. 

1.2. Early Developments 

Prior to the advent of space astronomy, an ingenious 
insight on the probable existence of a large volume of 



space that contains hot gases was provided by Spitzer 
( [1956| l, who noted that cold clouds found in the Galac- 
tic halo had to have some means of confinement. He 
speculated that such clouds were stabilized by the ex- 
ternal thermal pressure of a hot Galactic corona. It 
was not until some 18 years later that firm evidence 
on the pervasiveness of this type of material within the 
Galactic disk was demonstrated by direct observations 
of absorption by the O VI doublet at 1032 and 1038 A 
in the spectra of 32 stars ( [Jenkins & Meloy | 1974| York 



19741, accompanied by a recognition by [Williamson 



et al 



( 1974| ) that the diffuse soft X-ray emission ob- 



served by instruments on sounding rockets (Bowyer 



let al.||1968[[Bunner et al.|197T]|Davidsen et al.|1972| 



Yentis et al.|1972[ l probably came from the same type 
of material. 



Following the initial discoveries, Jenkins (1978a I 
(hereafter J78) conducted a survey of O VI absorption 
which increased the accumulation of cases to 72 Milky 
Way stars. From this body of data, he offered some 
conclusions on the nature of the hot gas and how it is 
distributed in space ('J enkins[[197 8b'c'). In particular, 
the statistical fluctuations of absorption strengths and 
velocities from one sight line to the next indicated that 
approximately 6 hot gas regions, each with a column 
density A^(0 VI) « lO'-' cm"^, are distributed across 
a representative distance of 1 kpc. Later, inspired by 
evidence from the soft X-ray surveys that the Sun is 
surrounded by a volume of hot gas out to a distance 



of around lOOpc ( Sanders et al.|1977t McCammon & 
[Sanders|[T990l [Snowden et al.|[1990[ [T998 r 

a struc- 
ture now called the Local Bubble (hereafter, 'LB'), 
Shelton & Cox ( 1994[ l conceived of a slightly more 



elaborate model that accounted for a possible constant 
foreground contribution of O VI, with A^(0 VI) aver- 
aging about 1.6 X lO'^'cm"^, which reduced the ap- 
parent fluctuations of a much more patchy distribu- 
tion at greater distances. Their re-interpretation of the 
O VI data suggested that outside the Local Bubble the 
Galactic disk is populated by regions that have column 
densities in the range 2-7 x lO'-'cm"^ with separa- 
tions 450 - 1300pc. With the advent of the Far Ultra- 
violet Spectroscopic Explorer (FUSE), 'Oegerle et al.[ 



( 2005) looked at weak O VI lines towards nearby white 
dwarfs (WDs), and determined the average volume 
density to be 2.4 x 10"^ cm"^. They concluded that 
the local O VI contribution is much less than that esti- 
mated by [Shelton & Cox[ ( [1994| l unless it is confined to 
a thin shell at the edge of the LB, i.e., just beyond the 
white dwarfs. A much larger survey of nearby WDs 



was conducted by Savage & Lehner ( 2006 1, who mea 



sured the average O VI density to be 3.6 x 10 cm ^. 
A challenge that was presented to the early O VI 



investigators was a proposal by Castor et al. ( 1975 1 that 
all of the O VI arose from bubbles that were created by 



stellar mass-loss winds from the target stars. Jenkins 
( [1978c| l addressed this issue and concluded that some 
O VI could arise from the bubbles, but that most of the 
hot gas was truly interstellar in nature. We return to 



this topic once again in this paper (§3.4 1 



After the end of the Copernicus era, there was a 
long period where only the transitions of N V, Si IV 
and C IV could be seen in absorption using the In- 
ternational Ultraviolet Explorer (lUE) and the Hubble 
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Space Telescope (HST). Because these species can be 
created at lower energies (ionization potentials of the 
next lower stages are only 33, 48 and 77 eV, respec- 
tively), there was a persistent concern that photoion- 
ized gases could be responsible for these species, espe- 
cially in low density environments. Even so, important 
conclusions emerged about the distributions of highly 
ionized material in the Gal actic halo ([Savage & Massa] 



T987l|Sembach et al. |1991| [Sembach & Savage|1992|i, 



the Magellanic Clouds (de Boer & Savage 1980 Fitz 



patrick & Savage'l985), and within the Galactic disk 
ICowie et aL_198 lb ; Walbo m etal.|1984 HEdgar & Sav-| 
age 1992; |Spitzer& Fitzpatr ick[T992l [ Sembach 1993] 
"1994; Tripp et al.'1993; Sav age et al.[1994, ,Sembachl 
et al . 1994; Huang et al. 1995|, often using distant UV 



sources fainter than the sensitivity limit of Coperni- 
cus. Further studies of O VI, however, were confined 
to results that were more limited in scope (|Ferguson| 
eraL][T995l [Hurwitz et al. 1995; [Dixon F t al. 1996) 



Hurwitz & Bowyer|1996(,Widmann et al.|1998(|Sem"^ 



bach et al.||1999| l because they came from short-lived 
missions, such as the Hopkins Ultraviolet Telescope 
{HUT) (' Davidsen|1993| l and the Orbiting Retrievable 
Far and Extreme Ultraviolet Spectrograph (ORFEUS- 
SPAS) ( [Hurwitz et "aL| [TWH] ). With the availibility 



of FUSE, we can now return to large-scale surveys 
of O VI absorption. Already, interstellar O VI has 
been studied in the Galactic halo (Oegerle et al. 2000 
Richter et al.|200Tl [Howk et al.'2002a; Fox et al. 2004 



Savage et al.|2003| !Sembach et al. 2003; WakkeretaD 
2003| ), Large Magellanic Cloud (jHowk et alT||2002bt , 



Small Magellanic Cloud (Tloopes et al. 20 d2|i^ and the 



Local Bubble (Oegerle et al. 2005; Savage & Lehner 
[2006) ). Likewise, intergalactic O VI has been found 
for absorption systems at redshifts that are too low 
for coverage with HST ([Oegerle et al. 2000) [Sembach 



et al. 2001 ; Tripp et al. 2001 ; Savage et al. 2002; Jenk- 
ins et al.ii2003| frochaska et al. 2004, ,Richter et al. 



20041 [SembacheT al. 2004; Danf orth & S hull 2005|. 



With the current paper, we cover yet another important 
realm, the Galactic disk. 

Through the years, theoretical investigations kept 
pace with the observations, starting with initial assess- 
ments by [Shapiro & Field] ( |1976| on the nature of a 
possible Galactic fountain, and an early proposal by 



Cox & Smith ( 1974 1 on the origin of a pervasive net- 



work ("tunnels") of hot gas that is kept alive by the 
channeling of shock waves (from supernovae) within 
it. This picture was investigated in more detail by 
[Smith[ ( ]W77] l. [McKee & Ostriker[ ( ]T977] l presented a 



seminal treatise on the different phases of the medium, 
accounting for mass and energy balance between a 
hot interstellar medium created by supernovae and the 
more conventional phases at lower temperatures. This 
theory was a very important advance over an earlier 
model for ISM phases developed by [Field et al.[( [1969] ) 
that overlooked the hot phase and its secondary ef- 
fects. Some noteworthy consequences that might arise 
from the hot gas production include mechanical energy 



transport to the embedded cool clouds (Cox 1979 1 and 
a feedback mechanism that may regulate star forma- 
tion in disk galaxies ( [Cox|198l) . Enhanced accumula- 
tions of O Vl-bearing gas will undoubtedly arise from 
hot gas at the edges of, and within, the interiors of su- 
pernova remnants, which is discussed in the models 
presented by [ Chevalier j([19741l, jCowie et al] ([1981all, 
ISlavin & Cox (1992 ), Edgar & Cox[ ( [1993| l, |smith~& 
,Cox ( 2001 ), and Shelton (2006). 

More recently, exciting realizations of how the ISM 
may be configured from supernova explosions have 
been constucte d by . de Avillez, (2000), jde Avillez~& 
[Breitschwerdt[ ( l2004| and de Avillez & Breitschwerdt 
( 2005a| l using hydrodynamic and magnetohydrody- 
namic simulations. Their models predict global equi- 
hbrium in the disk-halo circulation when the explo- 
sions in the ISM are traced over several hundred Myr, 
and demonstrate how effective turbulent diffusion can 
be in mixing hot and cold gas. [de Avillez & Breitschw-[ 
[erdt[ ( [2005 b') predict specifically the expected O VI col- 
umn densities along random lines of sight over dis- 
tances of up to 1 kpc, and find good agreement with the 
FUSE WD and Copernicus data. Their models predict 
a clumpy turbulent ISM in which O VI vortices and 
filaments are interecepted every ^100 pc, each with 
A^(0 VI)~ 1-2 X lO'^cm-^. 

1.3. Current Topics 

A persistently debated unknown quantity is the vol- 
ume filling factor of the hot gas in the Galactic disk 
( [Ferriere[|l995 ). Equally important is the topology of 
the gas: is it confined to isolated, compact bubbles or 
does it comprise a matrix of multiply connected re- 
gions within which the rest of the ISM is confined? For 
observers, it is easy to measure the average amount of 
O VI in the Galactic disk, and with reasonable assump- 
tions we may relate this quantity to the total amount of 
gas in the temperature regime 10^ - 10^ K. However, a 
much more difficult task is to understand the structure 
of the gas along the viewing direction, i.e., informa- 
tion that is lost by the projection of everything along 
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the sight Hne, which is compounded by the restriction 
that not all directions in the sky can be sensed. |Jenk-| 
|uis| ( |1978b ) placed an approximate upper limit of 20% 
for the filling factor from the apparent lack of an an- 
ticorrelation between the amount of O VI and denser 
(mostly neutral) hydrogen, as revealed by the color ex- 
cesses of the target stars (however, this argument can 
equally well be used to say the filling factor is greater 
than 80%). Models that propose that the hot gas has a 



very large filling factor [e.g. that of McKee & Ostriker 
1977 1] must be restrained by observations by'Heiles 



1980| l that indicate that the neutral, warm medium 
(gas that emits 2 1 -cm radiation but has a high spin tem- 
perature and thus negligible absorption) has its own 
significant filling factor One feature of the McKee & 
Ostriker model is that it assumes that supemovae occur 
at random locations and time. Actually, supernovae 
ai-e highly coiTelated ( [McCray & Snow|[T979[ |Heiles 
|1990[ ), and this should influence theories on the per- 
vasiveness of O Vl-bearing material, leading to fewer, 
but larger volumes whose interiors hold the hot gas. 
The existence of such large, coherent structures seems 
to be validated by observations of enormous shells of 
gas seen in 21 -cm line emission (Heiles 19791 11984{ 
IMcClure-Griffiths et al.|2002| l. 

The kinematics of O VI and its relationship with 
other forms of gas, those that are highly ionized and 
those that are not, can provide important clues on the 
nature of the hot gas on a more microscopic level. 
One of the first conclusions to emerge from the ve- 
locity profiles of the O VI features is that they were 
broader than those produced by less ionized species, 
such as Si II and Si III ( |York||1974"l ), indicating that 
the O VI was not simply a more ionized counterpart of 
the lower excitation gas. The absorptions exhibited by 
York were of very good quality, which enabled him to 
place upper limits to the Doppler motions that corre- 
sponded to a range T = 4x 10^ to 2 x 10^ for differ- 
ent cases. The most striking conclusion from the sur- 
vey by I Jenkins] ( | 1 97 8b| l was that the scatter of velocity 
centroids, amounting to only 26 km s"', is far below 
the shock speeds needed to create the post-shock tem- 
peratures that can strip 5 electrons from the oxygen 



atoms, i.e., at least 130- 160 km s ('Raymond] 1979 



|Shull & McKee|1979| l A similarly small disp er sion of 
21 km s~ was reported recently byjSavageetaL ( 2003 1 
for sight lines through the Galactic halo. The assertion 
that O VI does not have extraordinarily large veloci- 
ties must be viewed with caution however. We can not 
rule out the presence of some additional O VI that is 



moving very rapidly, for if the absorption profiles of 
such high velocity material were very broad and shal- 
low, they would be very difficult to detect against the 
undulating stellar continua. 

Simply put, O VI is not created in zones that im- 
mediately follow the shocks that heated the gas to high 
temperatures. Instead, the low velocities indicate that 
what we see arises either from some interaction with 
the normal, low-ionization material or, alternatively, 
perhaps ions whose velocities are regulated by mag- 
netic fields in the disk (Roger & Landecker'1988) that 
might have been compressed and amplified by the mo- 
tions of expanding supernova remnants ( |Spitzer|1990| l. 
For the former possibility, we may consider two alter- 
natives. One is that the O VI originates in transition 
layers where heat flows by conduction across inter- 
faces established between regions that hold v ery hot 



] (r > 10''K) gas and denser, cooler mater ial (jCowie 



& McKee 1977| |McKee & Cowie 


1977 




Begelman 


& McKee 1990'; 'Bertoldi & McKee| 


1990 




McKee & 



Be gelman 1990 ). Material is transported between the 
regions: at early times cool gas evaporates into the 
hot region, whereas later, after radiative cooling be- 
comes important, hot gas condenses onto the cool gas 
( |Borkowski et al.|1990| l. If there is a substantial differ- 
ence in velocity between the cool and hot material, an 
additional means of transporting heat and matter op- 
erates: instabilities in the transition zone create tur- 
bulence which acts as a source of mechanical mixing 
of the cool and hot phases (Begelman & Fabian, 199(){ 
ISlavin et al.|T993] l. 

Observational support for links between O VI- 
bearing gas and the cooler material was provided by 
|Cowie et al.| ( |1979] ), who noted correlations in velocity 
between O VI and the ions Si III and N II. Addi- 
tional support was provided by the good velocity cor- 
respondence between O VI and C II absorption along 
short path lengths through the local ISM ( [Savage &| 
|Lehner|[2006) . One way to gain insights on the rel- 
ative importance of the different processes discussed 
above, as well as to see if the gas is simply radiatively 
cooling ( Edgar & Chevalier]! 198 6) instead of being 
cooled by contact with low temperature material, is to 
test whether or not various other highly ionized atoms 
(N V, Si IV and C IV) agree with theoretical predic- 
tions for the different regimes ( Spitzer|1996 Fox et al. 
[20031 ). 

In addition to the possibilities listed above, we must 
be watchful for another physical condition, that, under 
certain circumstances, may apply to the regions that 
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are responsible for O VI absorption. In a discussion 
of hot gas within the Local Bubble, 'Breit schwerdt &| 
|SchmutzIer| ( |1994] l considered the possibility that an 
explosion in a dense region followed by expansion in 
a lower density medium would lead to adiabatic cool- 
ing of the hot gas. The lag in recombination as the gas 
cools would be even more severe than for a radiatively 
cooling gas. They envision the possibility that we are 
surrounded by gas at T ^ 50, 000 K that still has an 
appreciable concentration of ions that normally appear 
in equilibrium at much higher temperatures. It is im- 
portant to note that we cannot depend on identifying 
this gas on the basis that the lines have a velocity dis- 
persion characteristic of the Doppler broadening at the 
lower temperatures, since the expansion process itself 
creates a broadening that replaces the thermal broad- 
ening. 

In the discussion presented up to now, we have 
considered a number of contemporary problems in 
the study of gases responsible for O VI absorption in 
the disk of the Galaxy. More comprehensive discus- 
sions are provided by review papers writte n by|Cox & 



Reynolds' (T987ll, |Savage| ( fTM?) !, [Jenkins] ( |1987| l, and 
'Spitzer(1990j . 



1.4. Outline of This Paper — A Summary 

The capability of FUSE to record spectra of stars 
much fainter than the magnitude limit of Copernicus 
(V < 7.5), allows us to probe regions well beyond the 
original O VI survey of J78. In this paper, we present a 
new FUSE sample of disk stars which has about twice 
as many stars as the older survey. These two consid- 
erations allow us to overcome some shortcomings of 
the earlier work and address such issues as the possi- 
ble presence of very large scale structures containing 
either more or less than the usual amount of O VI, and 
whether or not the apparent decrease of O VI at mod- 
erate distances from the Galactic plane agrees with an 
exponential scale height of ^ 2.3kpc as determined 
from sources well outside the Galaxy ( [Savage et al.| 
2003 1. Also, with the larger column densities that arise 



from the greater lengths of the sight lines, we should be 
less sensitive to distortions of the results arising from 

i) the presence of bubbles around the target stars and 

ii) an over-representation of the region surrounding the 
Sun that arises from the constraint that all sight lines 
must emanate from a single point. 

This paper focuses on three principal topics. First, 
we describe in detail the data collection and reduction 



of the FUSE spectra, and the analysis of the O VI 
absorption line profiles. Second, we present a re- 
evaluation of the distances to the targeted stars, which 
is an important parameter for understanding how O VI 
is distributed in the Galactic disk. Third, we discuss 
the global properties of O VI absorption in the disk, 
which, of course, is the main science goal of this sur- 
vey. 

As far as the first of these topics is concerned, the 
data reduction of FUSE spectra is not unlike that of 
any other data taken with modern photon-counting de- 
vices. We were concerned, however, that we should 
be able to precisely calibrate the wavelength scale of 
the spectra, since FUSE does not take calibration spec- 
tra contemporaneously with an object. This led us to 
compare the data from many of our sight lines with 
high quality, high resolution HST data, from which we 
found that the standard FUSE wavelength scale for the 
pipeline calibration (CalFUSE 2.0.5) required a cor- 
rection of 10 km s"'. This analysis is introduced in 
§2.4 and described in detail in Appendix [A] 

In analyzing the O VI lines, we faced two main ob- 
stacles. First, J78 realized that any O VI A1032 line is 
contaminated to some degree by the HD 6-0 R(0) line 
at 1031.9 A. Before measuring the physical parame- 
ters of an O VI line, the contribution from this HD 
line must be removed. Fortunately, in the FUSE data, 
there exist other HD lines that are free of the dense H2 
'forest' and which can be used to model the contam- 
inating HD line. We found that removing the model 
using Voigt profile line-fitting techniques was highly 
successful and actually presented few problems for the 
subsequent analysis of the O VI A1032 line. Details 
are given in ^2.5| and [2.7[ ^3. 3 [ also discusses a consis- 
tency check on the efficacy of the HD removal using 
multi-epoch data. 

A more problematic obstacle in analyzing the O VI 
lines arose in determining the shape of a star's contin- 
uum. With the data reduced, we quickly discovered 
that the gradient of the continuum was probably vary- 
ing on scales less than the width of an O VI line. The 
most extreme variations were prevelant among the 09- 
Bl type stars. It became clear that we needed to track 
three possible continuum fits: the adopted best fit, and 
two possible variations above and below this fit — 'up- 
per' and 'lower' continuum fits, as we will refer to 
them throughout this paper Although measuring con- 
tinuum errors in absorption-line studies is far from un- 
usual, in many cases, the interstellar lines under study 
are quite narrow, and continuum errors are largely ir- 
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relvant. For our FUSE data, however, we found that 
continuum errors actually dominated the O VI column 
density errors, and not the usual Poisson counting sta- 
tistical errors. We discuss how we chose our continua 



in 5 2.6 With the spectra normalized, we were able to 



derive column densities, Doppler parameters and line 
velocities using both Voigt profile fitting ( ^2.8. l| l and 
the Apparent Optical Depth (AOD) method ( ^2.8.2| i, 
both of which gave very similar results, as expected. 
For the subsequent analysis, we assumed that the ab- 
sorbing gas had no small-scale velocity structures hid- 
den within the broad profile that might invalidate ei- 
ther analysis method ( ^2.9[ ). We favored the results 
from the Voigt profile fitting procedure, but used AOD 
measurements when no line was clearly present and an 
upper limit to the O VI column density was required 
( ^2.8.2 and Appendix [D|. The spectra themselves are 
presented in ^ 

The second principal goal of this paper is to re- 
evaluate the distances to the stars used in our survey. 
There are several databases now available (the Hip- 
parcos, Tycho, and 2MASS catalogs) which provide 
a more reliable estimate of a distance than was possi- 
ble a decade ago. We introduce the problem first in 
|3.2| where we show the distribution of stars in the 
Galactic plane. We expound a full discussion in Ap- 
pendix |B] since the methods used to derive distances 
are technical and outside the main scope of our discus- 
sion on O VI in the Galactic disk. A reader disinter- 
ested in the mechanisms of deriving stellar distances 
can certainly skip this Appendix, and instead simply 
note that the most egregious error in determing a dis- 
tance arises from the difficulty in determing the abso- 
lute magnitude of a star. My- This is not because of 
a problem in knowing My for the particular spectral 
type or luminosity class of a star, but because of the 
problems inherent in defining the type and class cor- 
rectly in the first place. The error in My is not only 
a function of the type and class, but depends on the 
quality of the stellar spectra used to define these pa- 
rameters, and on the skill and experience of the ob- 
server making the measurements. Despite the signifi- 
cant improvements in measuring various stellar param- 
eters with modem techniques, we believe that errors in 
stellar distances range from 10-30%, and that even 
larger deviations are possible. A full explanantion is 



given in Appendix B.3 



introduced in the previous sections, but we note the 
following. With the huge amounts of data available 
from satellites that cover many different regions of the 
electromagnetic spectrum, it is now possible to com- 
pare the results from absorption-line surveys with re- 
sults from other types of studies. In this paper, we con- 
centrate on X-ray emission data taken with the Ront- 
gen Satellite (ROSAT). Although X-ray emitting gas is 
nominally at a higher temperature than we would ex- 
pect for O VI absorbing gas, one can imagine a natural 
confluence of hot gas in structures that might be de- 
tected in both O VI absorption and X-ray emission. In 
P.4| and [3.5| we describe the types of X-ray emission 
in ROSAT maps seen towards our FUSE stars. In many 
cases, these maps indicate regions of hot gas along the 
line of sight from bubbles of gas around stellar asso- 
ciations or supernova remnants. These observations 
lead us to classify the morphology of X-ray emitting 
patches (or absence thereof) around or near our target 
stars. 



Using these classifications, we show in 5 3.7.3 that 



the average line of sight density A^(0 VI)/c/ is some- 
what higher for sight lines which pass through X-ray 
enhanced regions, leading us to conclude that a small 
fraction of A^(0 VI) towards some stars does indeed 
arise in circumstellar material or wind-blown bubbles. 
This information is important, because by selecting 
sight Unes with no clear X-ray emission, and by cor- 
recting O VI column densities for absorption by the 
Local Bubble using recently published data, we can 
derive a less biased value of the mid-plane density of 
O VI, no. We discuss this value at some length in 
P.7| concluding with our best determination of «o and 
its likely errors in ^3.8.3| In calculating these val- 
ues, we quantify the significance in the difference in 
A^(0 VI) towards stars in the north and south Galactic 



hemispheres originally found by Savage et al. ( 2003 1 
( |3.8.1| l. The difference has no effect on the value of 
«() we derive, but does produce two different values of 
the scale height of the O VI gas layer in the plane of 
the Milky Way. 



The final topic of this paper is, of course, the global 
properties of O VI absorption in the disk of the Milky 
Way. Our motivation for the study has akeady been 



In 5 3.9 we show that A^(0 VI) is directly correlated 
with the distance to a star in the plane of the disk. 
This correlation can most easily be understood as the 
effect of longer sight lines passing through more ab- 
sorbing clouds. This shows that O VI absorption can- 
not arise primarily from the circumstellar environment 
of the observed stars, but must be a truly interstellar 
phenomenon. The data also reveal a discontinuity in 
A^(0 VI) with distance beyond a few hundred pc, in 
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that column densities appear to be too high for a given 
distance. As we show, however, this arises from excess 
contributions to the measured A^(0 VI) from the Local 
Bubble. 

However, the data also suggest that the hot O VI- 
absorbing ISM is neither a smoothly distributed homo- 
geneous gas layer, nor a population of randomly dis- 
tributed clouds with similar properties. First, there is 
no acceptable fit to the decrease in O VI density with 
height above the plane of the Galaxy unless an addi- 
tional 'error' is introduced, a term which represents 
an intrinsic 'dumpiness' in the interstellar medium 
( ^3.8.2| i. Second, this dumpiness does not decrease 
with distance, as would be expected if a sight line sim- 
ply passed through more a greater number of generic 
clouds (5 3.9 1. To explain this, we propose a picture 



of the O VI absorbing ISM that contains a mixture 
of cloud sizes, with small, randomly distributed small 
clouds and sparsely distributed large clouds. We also 
note an alternative explanation, that we may be seeing 
the clustering of many small clouds. 

There are other lines of evidence that suggest that 
the O VI absorbing 'clouds' are actually one compo- 
nent in a complex, evolving, turbulent ISM. The ve- 
locities of the bulk of the O VI absorption do not fol- 
low the general pattern of differential co-rotation in the 
disk of the Galaxy ( ^3.11 1. Moreover, if we compare 



the extremes of the O VI absorption line profiles with 
the extremes of other lower ionization species (O I, 
C II, C III and Si III) we find that they match each 
other quite closely (5 3.12i. This apparent 'coupling' 
of at least some of the components in the absorption 
line profiles (at least those with the highest peculiar 
velocities along a sight line) suggests that the struc- 
tures giving rise to absorption from hot gas are multi- 
ionized. 

A natural explanation for all these phenomena has 
been discussed above. The numerical simulations of 
|de Avillez & Breitschwerdt] ( |2005b| and references 
therein) vividly realize how the ISM may look when 
driven by supemovae explosions. Their two dimen- 
sional density maps (for example) in the plane of the 
Galaxy show how hot gas arises in bubbles around su- 
pernovae, which is then sheared through turbulent dif- 
fusion, destroying the bubble and stretching the hot 
absorbing gas into filaments and vortices that dissi- 



pate with time. In 5 3.14 we confirm the relationship 
between A^(0 VI) and Doppler parameter found by 
other authors, and determine that the O VI lines are too 
wide to arise from differential Galactic rotation alone. 



Clearly, the velocity fields that govern the O VI com- 
ponents along a sight line are much larger than can 
be reproduced by the simple motions of the Galactic 
disk, and of course, expanding supernova shells natu- 
rally explain these large velocities. When considering 
an ensemble of such regions along long sight lines, it 
seems plausible that the correlation of b with A^(0 VI) 
could arise from the existence of exceptionally large 
internal motions inside large clouds, or that smaller 
clouds within clusters are moving more rapidly than 
the gases within individual clouds. 

|de Avillez & Breitschwerdt| ( |2005b| l found good ap- 
parent agreement between their models and the FUSE 
WD data and the results from Copernicus over a path 
length of ^ 1 kpc. Our results are also likely to be 
consistent with these models, although a detailed com- 
parison is beyond the scope of this paper. The data 
presented herein should provide sufficient information 
to test these, and any other detailed simulations, over 
larger scales. 

2. FUSE OBSERVATIONS, DATA REDUCTION, 
AND ANALYSIS 

2.1. Sample Selection 

The core of our sample of sight lines comprises 
stars observed under the auspices of the FUSE PI 
Team-time programs PI 02 and PI 22. Stars for these 
programs were originally selected based on a variety 
of criteria: they were more than 1 kpc distant, in or- 
der to probe regions of the Galaxy beyond those stud- 
ied by Copernicus; they had low values of redden- 
ing [E{B-V) < 0.3] to minimize contamination by 
H2 lines; they lay within ±10 degrees of the Galactic 
plane to confine the study to the Milky Way disk; they 
had spectral types of 09 or earlier, or 09 to B2 if their 
projected rotational velocities (vsin /) were known to 
be > 100 km s"'; they had luminosity classes of III 
to I for 09 spectral types (or earlier), or class V if 
vsin / was high; their spectra were not known to show 
emission line characteristics; and there existed ancil- 
lary data of the star from previous UV satellites such 
as lUE or ORFEUS-SPAS, or from ground-based opti- 
cal observations. 

In addition to the stars selected for the PI Team- 
time programs, we added data from the FUSE Multi- 
mission Archive at the Space Telescope Science Insti- 
tute (MAST). To investigate a range of ISM environ- 
ments, we took data from PI Team programs P116 and 
P216 ("Molecular Hydrogen in Translucent Clouds"), 
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PI 17 ("Hot Stars"), and PlOl ("The Properties of Hot 
Gas in the Milky Way and Magellanic Clouds"). In 
all cases, we selected only sight lines within \b\ < 10° 
of the Galactic plane, again, in order to focus primar- 
ily on absorbing gas in the disk of the Milky Way. We 
only included data that had signal-to-noise (S/N) ratios 
at least as high as those obtained for stars in the core 
program and we rejected stars that had highly irregular 
continua (see §2.6[ ). 

We also searched the FUSE MAST archive to find 
stars from Guest Investigator (Gl) programs that might 
also provide suitable targets. We adopted criteria very 
similar to those used for selecting targets for the PI 02 
and PI 22 programs; we did, however, exclude regions 
inside the Vela supernova remnant region, at Galactic 
longitudes 260° < I < 270°, and the Carina Nebula, 
280° < / < 300°. We also avoided any FUSE pro- 
grams whose science goals were obviously oriented 
towards studying O VI absorption along especially se- 
lected sight line^ We eventually collected all suitable 
data available prior to April 2003. Our sample com- 
prises 148 stars in total, whose parameters are Usted in 
Table[T] To introduce this table, we note the following: 
column 1 is a unique identifier ("unique ID") assigned 
solely to help cross-reference stars in different figures 
and tables in this paper. Column 2 gives the name of 
the star, while column 3, labelled 'Prog ID' is the com- 
bination of FUSE program ID (the first three numerals) 
and a number assigned within that program (the last 
two numerals) for a particular star. Columns 4 and 5 
list Galactic coordinates / and b; for many of the tables 
in this paper, the sight lines are sorted by Galactic lon- 
gitude, in order to map any variations in the physical 
characteristics of the O VI absorption along the plane 
of the Milky Way. Table |2] lists all the stars sorted by 
HD number instead, and includes the unique ID num- 
ber (as well as I and b) which can be cross-referenced 
with entries in Table [T] and any other tables or figures. 
The remaining entries for Table[T]are discussed in later 
sections below. 

2.2. FUSE Observations 

Of the 148 stars in our sample. 111 were observed 
as part of Pi-Team programs P102 andP122. Ajournal 



of observations is given in Table [2j which includes the 
observation date and exposure time. For some stars, 
more than one day passed between return visits to a 
star, in which case we list the second (or third) date 
of the exposure in column 6. There are also exam- 
ples where a PI Team program star was reobserved by 
a Guest Observer ('GO'; the first letter of these pro- 
grams run from 'A' to 'D' for Cycles 1-4); these are 
also listed after the Pi-Team program entry. In a few 
cases, some stars were observed with very large time 
intervals between satellite visits. These data are use- 
ful for determining how strongly our results are in- 
fluenced by possible time variations in the underlying 
stellar spectra. This in turn gives us an indication of 
how accurately we are able to normalize a star's con- 
tinuum around an O VI absorption line. We discuss 
this topic more fully in ^3.3| 

For each star, the data obtained usually consisted of 
several individual spectra (hereafter "sub-exposures") 
which needed to be co-added. The exposure time Tgxp 
in column 7 represents the total time actually spent on- 
target after the co-addition of useable sub-exposures; 
so, if individual sub-exposures were discarded (due, 
e.g. to loss of signal in a channel) the sub-exposure ex- 
posure time was not included in column 7. In total, just 
over 754 ksec was expended on programs P102 and 
P122. All stars were observed through the 30"x30" 
LWRS aperture, except HD 175754, HD 041161 and 
HD 069106, which were observed with the 4"x20" 
MDRS aperture. 

Details of the FUSE satellite and instrumentation 



^As we shall see, our Team-time programs P102 and P122 did include 
stars toward both the Vela and Carina regions, and those data are 
presented in this paper; however, by the time we searched the FUSE 
MAST archive it had become clear that the absorption towards these 
regions is far more complex than that seen in the rest of the Galactic 
disk (see ^3.7.2| and requires a separate analysis. 



are given by Moos et al. (2000i and Sahnow et al. 
( |2000) . The wavelength regions of interest for detect- 
ing O VI are covered by four separate mirror/grating 
and detector combinations, namely the LiFlA/2B and 
the SiClA/2B channels. However, we made measure- 
ments mainly using only LiFlA data, because that 
channel provides the highest effective area, the astig- 
matism correction point was designed to be at 1030 A 
in the LiF channels, and the LiFlA data were ini- 
tially the best calibrated. There are two exceptions: 
the HD 041 161 and HD 060196 spectra came from the 
SiClA and LiF2B channels, respectively, as no LiflA 
data were available. 

2.3. Sub-Exposure Corrections and Co-Addition 

The raw data for all the stars in our survey were 
reduced using version 2.0.5 of CalFUSE, which was 
the most advanced version available during the data 
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reduction. Since the stars observed in this study are 
bright, the S/N ratios of individual sub-exposures 
were usually high enough for shifts between sub- 
exposures (caused by thermal motions in the spec- 
trograph) to be seen. To determine possible shifts, 
we used the sub-exposure with the highest S/N ratio 
near the O VI A 1032 line as the template spectrum 
against which other sub-exposures were compared. 
An initial shift of a sub-exposure was determined by 
cross-correlating a 10 A region of the data (centered 
at 1032 A) with the template. A more refined shift 
was then made visually, in order to determine a shift 
to within ~ ±0.5 pixels (~ 1 km s~'). Both flux and 
error an^ays were shifted using sine interpolation if a 
shift of more than ± 0.5 pixels was measured. Ideally, 
we would like to have moved spectra only by integer 
numbers of pixels to properly preserve the noise char- 
acteristics of a spectrum. However, shifts of less than 
a pixel were clearly evident in data with good S/N, so 
fractional shifts were required. (As we will discuss 
later, continuum uncertainties dominate the errors in 
our data analysis, so noise errors are only of secondary 
interest.) 

With individual sub-exposures mapped to the wave- 
length scale of the sub-exposure with the highest S/N, 
the wavelength and flux arrays of all the data were re- 
binned by a factor of 3, in order to increase the signal- 
to-noise of the spectra but still properly sample the 
data. The error arrays were rebinned by the same 
amount and reduced by a factor of \/3. The initial dis- 
persion of the data (from the CalFUSE calibration) was 
~ 0.0067 A pix"', 0but the resolution of the LiFlA 
channel is between 15 and 20 km s"' (0.052-0.069 A 
at 1030 A) FWHM. After rebinning, the final disper- 
sion of the data was 0.020 A pix"', close to that re- 
quired for optimally sampling the data. Individual sub- 
exposures were then co-added, with the flux weighted 
by the inverse of the variance formed from error arrays 
which were usually smoothed by 7 pixels. 

2.4. Wavelength-scale Accuracy and Zero-Point 
Correction 

An important quantity to measure in the FUSE 
spectra is the velocity of the detected O VI A 1032 
lines, the precision of which depends in large part on 



Note that this dispersion was less than the dispersion of 
0.013 A pix"' used by later versions (> 3.0) of CalFUSE. This new 
number samples approximately one-quater of a FUSE spectral reso- 
lution element (Dixon et al. 12007). 



the accuracy of the wavelength calibration. Unfortu- 
nately, FUSE has no on-board facilities to provide such 
calibration in situ. Instead, CalFUSE assigns a wave- 
length scale to the data, accounting for distortions in 
the dispersion and spatial directions, and correcting for 
instrumental distortions. 

Even with these corrections, errors arise in the zero- 
point of the wavelength calibrations, because, e.g., the 
source is not centered in the aperture, or, for channels 
other than LiFl, from misalignment of the primary 
mirrors. The LiFl channel has the most stable wave- 
length calibration, since for the observations in this pa- 
per, the Fine Error Sensor used light from this channel 
to acquire an object. However, since we found small 
offsets between sub-exposures taken during the obser- 
vations of a single star, we expected to find small errors 
in the absolute wavelength scale. (Indeed, although we 
used the sub-exposure with the highest S/N as a tem- 
plate against which to shift all the other sub-exposures, 
there was no reason to believe that the zero-point for 
this spectrum would be any more reliable than that of 
any other exposure.) 

Fortunately, CalFUSE adequately corrects for ex- 
treme non-linearities when mapping pixels to wave- 
length, particularly over the region of interest around 
the O VI A 1032 line. The measured velocities of the 
two molecular hydrogen (H2) lines which flank the 
O VI A 1032 line always agreed with each other (to 
within a Icr dispersion of ~ 2-3 km s"'), which val- 
idates this assumption. We therefore needed only to 
find and apply a zero-point offset to the spectrum to 
correct the wavelength scale. 

To derive a zero-point for a subset of the FUSE 
stars, we compared the velocities of H2 lines near the 
O VI A 1032 fine, with that of a CI I A 1347 line in 
Archival HST data of the same stars, where the wave- 
length calibration is more precise than that for FUSE. 
In all, we were able to correct 55 spectra out of the 
total 148 spectra using this method. A comparison 
between the CalFUSE wavelength scale and the cor- 
rected value showed that for the remaining spectra, 
a shift of 10 km s"' would correct the zero-point to 
within ± 5 km s"' of the true value. Full details of this 
analysis are given in Appendix [A| 

2.5. The O VI Wavelength Region and Contami- 
nating Lines 

An important step in the processing of the FUSE 
data was to normalize the spectrum of each star. Be- 



10 




1029 



1030 



1031 1032 
Wavelength (A) 



1033 



1034 



Fig. 1. — : The spectrum of HD 101190, show- 
ing wavelengths of possible contaminating absorption 
Unes around and within the O VI A 1032 line. 



fore describing this process, however, it is helpful to 
first discuss the identification of contaminating lines 
in the wavelength region around the O VI A 1032 ab- 



sorption line itself. As we show in 5 2.6 understanding 
the expected strengths of these lines was important for 
correctly determining the true shape of the continuum. 

The positions of possible contaminating lines 
are shown in Figure [Tj where we plot the the un- 
normalized spectrum of HD 101190 over the wave- 
length range 1029-1034 A. The wavelengths of the 
fines shown in the figure are given in Table |3] This 
table also lists other lines that were used in our study 
(often at wavelengths far from the O VI wavelength 
region discussed in this section) and which will be 
introduced in later sections of this paper. 

The most serious contaminant is the Lyman HD 
6-0 R(0) line at 1031.909 A, which is coincident in 
wavelength with the O VI A 103 1.927 fine if the ve- 
locity differences are only moderate. For strong O VI 
lines, the contamination is largely negligible, but for 
weaker O VI lines, the contribution of the HD line can 



be significant. We discuss in 5 2.7 how we removed the 
contribution from this HD line. 

CI I A 1031 is expected to arise at 1031.507 A. This 
line was often seen in the FUSE spectra, always nar- 
row and weak. We attempted to use the CI I A 1004.678 
fine to estimate the CI I column density and thereby 



model the line at 1031 A, since the 1004 A line is 
relatively free of other contaminating lines. Our at- 
tempts were unsuccessful, for two principal reasons. 
First, defining the rapidly changing continuum around 
both the 1031 and 1004 A lines proved quite difficult, 
a problem that made the depths of these intrinsically 
weak lines highly uncertain. Second, even in the few 
(rare) cases where the continua were well behaved, we 
could not fit both lines consistently with a single cloud 
model of column density and Doppler parameter, lead- 
ing us to suspect that problems may exist with pub- 
lished oscillator strengths. Recent work by |Sonnen-| 
[trucker et al.|p006| l confirms these suspicions, at least 
for the 1004 Aline. 

Hence, in our analysis of the O VI A 1032 line, we 
note where CI I A 1031 is expected, but no attempt has 
been made to model it. Fortunately, the position of the 
CI I line is usually sufficiently blueward of the O VI 
feature that its contamination is unimportant. There 
are sight lines where the O VI feature can become 
blended with the CI I line; however, as we discuss 



more fully in §2.8.1 the fitting of Voigt profiles to the 
O VI line is unaffected by the presence of CI I, and the 
contribution to the computed O VI column density is 
negligible. 

Close to the CI I line is the H2 8-0 R(8) line at 
1031.557 A. To predict if we should see this line, we 
looked at the 7-0, 6-0, and 4-0 R(8) fines at 1042.745, 
1054.520 and 1079.932 A, respectively. These lines 
are relatively free of contamination by other interstel- 
lar species, and have similar oscillator strengths (see 
Table [3]) compared to the 8-0 line, thereby providing 
suitable comparisons. We found absorption from these 
three R(8) transitions only along the line of sight to 
HDE 303308, which fies in the Carina Nebula. For 
those data we modelled the H2 column density for the 
8-0 line from the other three transitions and removed 
the profile from the spectrum. Weak R(8) lines may 
also be present in the spectrum of HD 187282, but 
there is significant blending with other H2 features. 
Fortunately, the velocity of the 8-0 line is well enough 
removed from the O VI that it poses no threat of con- 
tamination to the O VI absorption. 

As shown in Figure[T] O VI A 1032 is flanked by the 
two strong H2 fines at 1031.192 and 1032.354 A, the 
6—0 P(3) and 6-0 R(4) Lyman lines. These often blend 
with the O VI feature, but can be modeled and removed 
quite accurately, as we discuss in the following section. 
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2.6. Continuum Fitting 

Even with a complete inventory of the lines ex- 
pected near the O VI A 1032 line, normalization of 
the stellar spectra proved challenging. We adopted 
the method of fitting continua as described by |Sem-| 
|bach & Savage| ( [T992| l, whereby continuum levels are 
defined by least-square fits of Legendre polynomials 
to intensities either side of a line. Unfortunately, for 
many stars, we found that the stellar continuum varied 
in shape over widths comparable to that of the O VI 
features. 

No interstellar absorption lines are expected be- 
tween the H2 7-0 R(6) Hne at 1030.08 A and the first 
of the H2 lines flanking the O VI feature, the 6-0 P(3) 
hne at 1031.192 A. This gap of 1.1 A should provide 
~ 50 (rebinned) pixels from which the blue side of 
the continuum can be well defined. Except for O VI 
with high negative velocities, there was indeed a small 
amount of continuum available between the P(3) fea- 
ture and the O VI line in our spectra, again helping 
to tie down the continuum fit on the blue side of the 
O VI feature. There was, however, often contamina- 
tion from the CI I A 1031 line, potentially leading us 
to underestimate where the flux should lie between the 
P(3) line and the O VI. Fortunately, the CI I line was 
always weak and narrow, and its effect on causing con- 
tinuum misplacement was small. 

As we shall see, the O VI A 1032 lines themselves 
are wide, and there was often good reason to believe 
that the continuum was changing shape at 1032 A in 
many of our spectra. Further, at the red edge of the 
O VI line, absorption from the second flanking H2 fea- 
ture — the 6-0 R(4) line — is expected to set in. In 
many cases, we found that the O VI and the R(4) line 
were clearly beginning to blend together. By the time 
the continuum had recovered redward of the H2 ab- 
sorption, it was often a very different flux than that 
predicted simply by extrapolating from the (well fit- 
ted) continuum blueward of the P(3) line. 

We attempted to find some consistency checks for 
the continua that we adopted. As a first step, we com- 
pared the P(3) and R(4) lines with their profiles from 
modeling other H2 lines in the FUSE spectrum. For 
the R(4) line, we used the 5-0 and 4-0 R(4) lines 
at 1044.542 and 1057.376 A, respectively. These 
fines are relatively free of contaminating features, 
and have similar oscillator strengths to the 6-0 lines. 
For the P(3) line, we used only the 4-0 P(3) line at 
1056.472 A. In both cases, we fitted Voigt profiles 



(convolved with the instrumental line spread function 
— see |2.8.1| l to the data to derive column densities 
and Doppler parameters for the R(4) and P(3) transi- 
tions, calculated the expected profiles for the 6-0 P(3) 
and R(4) lines, and compared those with the observed 
lines. 13 

In Appendix|E|we show all the spectra in our survey 
(see Fig.[24l a full description of this Figure is given in 
AppendixjS. The top panel (the first of three) for each 
star shows the adopted continuum, along with the ± Ict 
error 'envelopes' which express the positive and nega- 
tive deviations possible in the polynomial coefficients 
used to define the continuum. (We will refer to these 
as the 'upper' and 'lower' continuum fits throughout 
this paper.) The middle panel shows the predicted pro- 
file fits to the H2 lines, blended with O VI hne fits 
( ^2.8. l| l, for the normalized spectra. Given the inher- 
ent uncertainties in producing the H2 models from the 
other H2 lines in the FUSE spectra, these fits appear 
to predict the observed features quite well. We notice 
a tendency for the theoretical profiles of the flanking 
R(4) and P(3) lines to be slightly stronger than the data 
in many cases. This probably arises from an incom- 
plete knowledge of the FUSE Line Spread Function 
(LSF) at the wavelengths used for the H2 models as 
well as the lines near O VI, since the LSF is largely 
uncharacterized for the spectrograph and is known to 
change with wavelength. We should say, however, that 
these models did not provide particularly rigorous con- 
straints on the continuum fits. Since the H2 lines were 
always quite strong, even fairly large continuum devia- 
tions around the lines produced very small differences 
in the depths of the line profiles. 

The most difficult continua to fit occurred when the 
O VI absorption arose in a region where the stellar flux 
is reaching some local maximum. For example, the top 
panel of Figure [2] shows the flux of HD 178487, which 
has one of the most difficult continua to characterize of 
all the stars in our sample. The solid line in the figure 
shows a plausible fit, but it is clear that the dotted lines 
could also represent the true underlying continuum. As 
we describe below, we can calculate column densities 
from the data normalized by our adopted fit (second 
panel of Figure |2]i but we can also re-calculate the val- 

'in fact, the range of oscillator strengths for all the R(4) lines avail- 
able in the FUSE channels are too small to allow a unique determi- 
nation of the R(4) H2 column densities, given the FUSE resolution. 
We base our use of the 6—0 line on the fact that we can correctly 
recreate the equivalent width of the line given the equivalent width 
of the other R(4) lines. 
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ues assuming that the upper or lower fits are correct 
(third and fourth panels). For HD 178487, the O VI 
lines are strong, and even a continuum as uncertain as 
the one indicated in the figure introduces errors of only 
~ ±0.1 dex in the total O VI column density. For- 
tunately, for all the stars with these type of continua, 
the O VI A 1032 line was indeed strong, and we were 
able to characterize the likely range of column densi- 
ties through the errors adopted for the continuum fits. 

A far more common problem in fitting continua can 
be seen for one particular star, HD 094493, in Figure[3] 
There are two possible continuum fits: one which con- 
tinues upward from an extrapolation of the fit blueward 
of the H2 P(3) line, before having to turnover beyond 
the R(4) (labelled 'Fit-2' in the Figure); and one which 
quickly drops downwards to meet the small amount of 
continuum flux left between the O VI line and the H2 
feature ('Fit-l')- 

The difference depends on whether Fit- 1 is the cor- 
rect continuum, despite the unsettling sudden devia- 
tion around the O VI A 1032 line, or if the stellar flux 
is instead depressed by a high velocity extension to the 
main O VI absorption. The second panel of Figure [3] 
shows the data normalized by Fit-l, and the resulting 
O VI absorption and H2 R(4) are easily reproduced 
with a single component fit. If Fit-2 is adopted instead, 
then a single component fit is no longer satisfactory. 
The third panel of Figure [3] shows the data normalized 
by Fit-2, with a single component fit to the O VI line 
(along with the predicted line profile for the H2 R(4) 
Une). There is clearly excess absorption between the 
O VI and H2 line. The blue side of the H2 line is also 
poorly matched. To account for this absorption, a sec- 
ond component is required at higher velocity. Such a 
fit is shown in the bottom panel of Figure |3] 

This type of problem occurs for 26 stars in our sam- 
ple. Perhaps significantly, the difficulties are confined 
to stars of a particular stellar type. Of the 26 problem 
stars, 21 (81 %) are of type 09 to Bl (inclusive). (In 
contrast, for the remaining stars, only 52% are of type 
09 to Bl.) In many cases, it seems that O VI absorp- 
tion arises in the blue wing of a narrow O VI emission 
line feature which predominates in these types of stars. 
There is no evidence of excess high (positive) veloc- 
ity O VI absorption in the majority of earlier spectral 
types. Nor are these problem stars confined to a partic- 
ular region of the sky, where, e.g. Galactic rotation (or 
some other coherent kinematical process) might be ex- 
pected to produce additional O VI absorption. We also 
note that when we compare the expected profile of the 
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Fig. 2. — : An example of the worst case for fitting 
a continuum to a star when the O VI and H2 6-0 
R(4) lines are strong and placed in a wavelength re- 
gion where the continuum reaches a local maximum. 
The solid curve in the top panel shows the adopted fit 
to the observed flux, while the dotted curves mark two 
extremes arising from the ilcr errors on the Legendre 
polynomial coefficients. All three are plausible fits to 
the data. The second panel shows the profile fits for 
the O VI and the molecular lines, given the adopted fit. 
Two components are needed to fit the O VI profile in 
this example, and the centroid of each is shown with 
a gray triangle. The lower panels show the fits for the 
upper and lower continua. However, because the ab- 
sorption is strong, the errors in the total O VI column 
density are only 0.1 dex. (See Appendix |E] for full de- 
tails of how this plot is marked.) 
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Fig. 3. — : A small subset of stars in our survey have 
continua whose shape can be fit with one of two possi- 
ble choices. Here we show one such star, HD 094493. 
In the top panel, Fit-2 is based on a simple extrapola- 
tion of the continuum from the blue side of the spec- 
trum, while Fit-1 assumes that the flux between the 

VI and R(4) line is real continuum. Adopting Fit- 

1 means that the O VI line is easily fit with a single 
component (panel 2; the centroid of the component is 
indicated by a gray triangle). Using Fit-2 results in 
additional O VI which cannot be fit with a single com- 
ponent (panel 3) but requires a second, broad, higher- 
velocity component (panel 4; the centroids of both 
components are indicated with gray triangles). Many 
of the stars in our sample have this ambiguous contin- 



uum shape, but for reasons discussed in 5 2.6 we adopt 
continua with shapes similar to those of 'Fit-l' above. 
(See Appendix |E] for full details of how this plot is 
marked.) 



R(4) line with the data normalized by Fit- 1 type con- 
tinua (marked as line '5' in the second panel of Fig- 
ure [3]l, we find that the predicted profile agrees well 
with the data every time. 

For these reasons, we adopt continuum levels of the 
type shown by Fit-1 for the 26 problem stars. Can we 
rule out the possibility that we are excluding real addi- 
tional, high velocity O VI components? If such com- 
ponents exist, then the absorbing clouds must have a 
rather unique set of properties: they must always be of 
lower column density than the lower velocity compo- 
nent; they must always be about twice as broad; and 
they must always be ~ 50 km s"^ redward of the lower 
velocity component. Such a fixed set of characteristics 
from sight lines scattered across the plane of the Milky 
Way seems contrived, and we believe that continua 
fits similar to Fit-2 in Figure |3] are simply incorrect. 
In many cases, the difficult continua are likely to be 
the result of irregularly shaped P-Cygni profiles which 
arise from the winds around the stars themselves. 

Table |4] details the results of our measurements of 
O VI lines towards all the stars in our survey except 
those for which we could not detect O VI absorption 
(see |2.8.2| for a discussion of these stars, which are 
Usted in Table |5]l. Of particular relevance to the prob- 
lem of continuum fitting is the information given in 
column 15 (labelled 'Cont'), a flag that simply warns 
of possible problems in the continuum fits. For the 
stars discussed above, where we could have higher 
continua and a second, broad, positive velocity O VI 
component, we assign a value of '2' for the contin- 
uum flag. There are also a few examples where the 
continuum could actually be drawn in such a way that 
additional O VI components might exist both redward 
and blueward of the O VI line. These are assigned a 
value of '1'. We consider it unlikely that these stars 
really have continua so very different from those we 
chose for them, but given the complex nature of stellar 
continua, we cannot rule out the possibility. In both 
cases, it is important to note that values for the contin- 
uum flag greater than zero are only assigned if a differ- 
ent observer could have adopted a continuum different 
from the one we adopted, by an amount which would 
not be covered by the continuum error envelopes intro- 
duced above. 

2.7. Removal of the HD 6-0 R(0) Line 

As indicated in ^23] Figure [T] and Table [3] O VI 
absorption occurs at nearly the same wavelength as the 
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HD 6-0 R(0) line at 1031.909 A. The line is a serious 
contaminant for weaker O VI lines, and must be re- 
moved before the properties of an O VI line can be 
measured. 

To remove the HD line, we first constructed a pro- 
file of the 6-0 absorption line based on a fit (see 
pXT) to another HD Une, the 7-0 R(0) feature at 
1021.453 A. This line has a similar oscillator strength 
to the 6-0 line (see Table |3]l, and is usually well iso- 
lated from other interstellar features. The resulting 
profile created for the 6-0 line should therefore be, at 
the very least, an accurate reflection of the amount of 
equivalent width needed to be removed from the O VI 
profile. 

We also checked to see if our fits to the 7-0 line 
were consistent with the 4-0 HD line at 1054.286 A, 
which also has an /-value close to that of the 6-0 
fine and is relatively free of contaminating interstel- 
lar lines. We often found that although the equivalent 
width of the 4-0 line (as predicted from the 7-0 line) 
had roughly the correct equivalent width, the profile 
of the 4-0 feature sometimes departed from a single 
component Voigt profile. We did not therefore use the 
4-0 line to constrain the HD model. 

Having derived a plausible model for the HD ab- 
sorption, we created a synthetic HD 6-0 line profile 
and removed it from the O VI absorption. The results 
of this removal can be seen in the top panels of Fig- 
ure |24] where the difference in the flux before and 
after the subtraction is shown as a gray shaded region. 
Towards several stai's (e.g. HD 210839, HD 005005 A, 
HD 012323, etc.) there is no apparent O VI absorption 
after the subtraction of the HD line. Significantly, for 
these cases, our HD subtraction never leads to the pro- 
duction of extra flux above the predicted continuum 
(i.e., an 'emission' feature), at least beyond that ex- 
pected from Poisson fluctuations. It seems likely there- 
fore that our HD models are sufficiently accurate to 
properly remove the HD line. 

2.8. Measurement of O VI Absorption Lines 

To measure the column density of the O VI absorb- 
ing clouds, we would, ideally, like to use both mem- 
bers of the O VI doublet. Both lines were recorded for 
most of the stars observed with Copernicus, since the 
H2 absorption along the sight lines was usually weak. 
For the FUSE stars, however, H2 absorption was al- 
ways too strong for us to use the 1037.6 A line. An 
example is shown in Figure |4j which shows the region 
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Fig. 4. — : Example of a stellar spectrum around the 
expected position of the O VI A 1037 line. The star is 
HD 67659, and shows the typical strengths of the H2 
features which make it impossible to measure the 1037 
line. 

around 1037 A for the spectrum of HD 167659. These 
data are quite typical for the stars in our survey . 

The O VI A 1037 could sometimes be seen in the 
wing of the 5-0 R(l) line when the H2 column den- 
sity was low. In principle, it is possible to fit 7 = 1 
fines in other parts of a spectrum, then calculate and 
remove the H2 line profiles around 1037 A and thereby 
recover the O VI A 1037 line. We performed such an 
exercise for several sight lines where the H2 column 
density was low, and the H2 lines weak. The results 
were disappointing. Although reasonable fits could be 
made for the set of H2 lines at, e.g, 1078 A and 1063 A, 
it was difficult to reproduce the H2 complex at 1037 A 
with sufficient accuracy (particularly without knowing 
the continuum a priori) to recover the O VI A 1037 line 
with anywhere near the precision needed to make the 
resulting absorption line reliable enough to validate 
our measurements of the 1032 A line. Indeed, lines 
of sight with a low H2 column density along the sight 
line tended to be at relatively short distances; the O VI 
was itself relatively weak in these cases, making an 
accurate determination of the 1037 A line even more 
difficult. 

Consequently, all the analysis of O VI in this pa- 
per uses the O VI A1032 line alone. In the following 
sections we summarize the two methods used to an- 
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alyze the O VI A 1032 lines toward our sample stars: 
the fitting of theoretical (Voigt) line profiles, and the 
Apparent Optical Depth (AOD) method. 

2.8.1. Column Densities: Profile Fitting 

As described above, a significant impediment to 
measuring O VI column densities towards some stars 
arose when the line was blended with contaminating 
CI 1 A1031 and H2 lines. The most straightforward 
way to circumvent this problem was to fit theoretical 
Voigt profiles to the regions of an O VI line which were 
clearly not blended with a contaminant. 

This method has been used by researchers for many 
years. Details of our specific line-fitting procedures 



are given in Bowen et al. ( 1995 1. For the FUSE spec- 
tra, we generated (one or more) theoretical Voigt line 
profiles from an initial guess of the Doppler parameter 
b, line velocity v, and column density A^. These were 
convolved with the instrumental LSF. A definitive LSF 
is not available for the FUSE channels, so we adopted 
a Gaussian of width 20 km s"' FWHM. Fortunately, 
the widths of most O VI lines are wider than this, so as 
long as the shape and width of the true LSF is close to 
the adopted value, the errors introduced are negligible. 

For example, using an LSF with a width of 
20 km S-' for the data of HD 134411 (#131), we find 
that the line has a Doppler parameter of 20.0 km s"' 
and column density 3.63 x lO'^* cm"^. This line is rep- 
resentative of lines with the smallest Doppler param- 
eters in our survey (see ^3.13 1 and would be most af- 
fected by changes in the LSF. In fact, if we had instead 
used an LSF of width 15 km s"', we would have mea- 
sured a Doppler parameter of 2 1 .5 km s"' — an increase 
of 7% — and a column density of 3.59 x lO'-' cm"^ — 
a decrease of 1%. For O VI lines which are wider than 
that seen towards HD 134411, the relative changes 
in column density and Doppler paramater are even 
smaller. 

Our fitting procedure allowed the OVIN, b and v to 
vary until a minimum was reached in the fit between 
profile and data using the POWELL minimization rou- 
tine ( Press et aLl|1992 1. Data points lying within the 
wavelength region predicted for the CI I line were ex- 
cluded from the fit. The O VI profile was often blended 
with the predicted profiles of the two flanking H2 lines 
— predicted in the sense that and b of the H2 lines 
were fixed at the values derived from fitting the rel- 



evant lines in other parts of the spectrum (see 52.6i, 
although V was allowed to vary for the H2 lines in or- 



der to take into account any small non-linearities in the 
wavelength scale of the data. 

For each star, we first tried to fit only a single O VI 
component to the absorption. Additional components 
were added only if single components were clearly in- 
adequate. It is important to note that this does not pre- 
suppose that absorption arises from a single cloud. As- 
suming such a model is likely to be erroneous, given 
the possible explanations for the origin of O VI ab- 
sorption discussed in ^and the results discussed later 
in this paper. However, in most cases the O VI line 
profiles appeared quite simple and shaped like a single 
component. 

Line fits were made three times for each star, and 
the results are given in Table |4] once for the data nor- 
malized by the best continuum fit (giving column den- 
sities and Doppler parameters b — columns 10 and 
5, respectively), and twice more for the upper (A^„, /?„ 

— columns 1 1 and 6, respectively) and lower (Ni, bi 

— column 9 and 4, respectively) error envelopes dis- 
cussed in ^2.6| For all sight lines, changes in v for each 
of the three continua fits are negligible (< 0.1 km s"' 
in most cases). This procedure quantifies how the fit 
to the O VI changes from errors in the continuum, but 
does not address errors arising from Poisson statistics. 
To estimate these values, we performed a Monte-Carlo 
simulation ( Bowen et al.|1995j ); for each sight line, we 
used the best fit values A^, b and v to generate a the- 
oretical line profile and added the amount of Poisson 
noise given by the FUSE error arrays (wherein each 
pixel is assigned a deviation cr,). This new 'synthetic' 
spectrum was then refitted, and new values of A^, b and 
V produced. This procedure was repeated 300 times, 
resulting in 300 different values of A^, b and v. In most 
cases, the distributions of these values were Gaussian 
in shape and the values of a for each were taken to rep- 
resent the errors from noise, <t(A^), (T(b) and cr(v). The 
values for the first of these two quantities are given in 
Table |4] (columns 12 and 7, respectively). We do not 
list values of a{v) since, along nearly all sight lines, 
they are only ^ 1 km s"', less than the error in the 
wavelength scale. 

In principle, these Monte-Carlo runs should have 
been repeated twice more for each star using the best- 
fit profiles derived from fits to the data normalized by 
the upper and lower continuum errors (i.e. [A^;, bi] and 
[A^„, buY). In practice, however, since the values of cr, 
at each pixel were the same, and the three profiles for 
each continuum fit were quite similar for most stars, 
the differences in (j(N), a{b) and a{v) were negligible. 
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For the O VI lines, which are much wider than 
the LSF, (j{b) and (j{N) are uncorrelated, and can be 
treated independently. The total errors in the Doppler 
parameters, (T{b)T are given by two values [bu-b,b- 
bi] added in quadrature to (T{b). These values are given 
in column 8 of Table [4] In a few cases, the posi- 
tive and negative values of (j(b)T in column 8 appear 
quite asymmetric. This usually occurred if the upper 
and lower continuum fits for a star were very differ- 
ent in shape from each other. The final errors in A^, 
a(N)T, are formed in the same way, since the two 
sources of error in are independent of each other, 
with [Nu-N,N-Ni] added in quadrature to cr(N). Col- 
umn 13 of Table |4] shows these errors in the form 
N[+a(N)T-a(N)T]. 

Towards some stars, more than a single component 
was needed to fit the O VI line. For convenience, col- 
umn 13 of Table |4] lists the sum of the column den- 
sities, AfroTAL^ along with their errors. In the case of 
a multicomponent fit, the final errors in Njojal are 
formed in a similar way to those of single component 
fits. Since the continuum errors affect all components 
in the same way, the error in AfpoTAL from the contin- 
uum fit is [{EiNi„-j:,Nd, (EiNi-EiNu)l where 
i is the component number. For the noise error, the 
column densities for all the components are added for 
each Monte-Carlo run, and a is calculated from the 
distribution of 300 values of X^/W)- '^gain, the con- 
tinuum error and noise error are added in quadrature to 
give the errors listed in column 13 of Table |4] 

Some stars from our survey are not listed in Table|4] 
They have no appreciable O VI absorption that can be 
fitted with a theoretical line profile. The limits to their 
column densities are tabulated in Table|5]and discussed 
below. 

2.8.2. Column Densities: Apparent Optical Depth 
(AOD) and Upper Limits 

Clearly, the Voigt profile fitting method described 
above is only applicable when an absorption line is 
clearly detected. In some of our spectra, there was 
good evidence that an O VI line was not detected. 
In order to rigorously define either the presence of 
very weak O VI hnes, or meaningful upper limits, 
we used the AOD method of analysis, which yields 
a straightforward formal measurement of A^(0 VI) and 
its associated error. A desire to accurately characterize 
marginal detections, or the lack of any detections, with 
a well-posed quantitative calculation, was the motiva- 



tion for carrying out AOD measurements. 

The AOD method is described in detail by e.g. |Sav-| 
age & Sembach (^1991^ and ,Jenkins, ( ,1996j . To sum- 
marize: the optical depth, t,, at a given pixel / (usu- 
ally in velocity space, where the width of a pixel is 
5v) is given by r, = -In where Fj is the normalized 
flux. If cr, is the error in the normalized data at the /'th 
pixel, then the corresponding error in t is cTt- = cr,//^ 
provided that cr, ^ f). The AOD column density as- 
sociated with this single pixel, A^,, in units of cm"^, is 
then given by 
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■5v= 3.768 X 10' 
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where / is the oscillator strength of the absorption line, 
Ao the rest wavelength of the line in A, 5v is the pixel 
velocity width in km s"' , and the other symbols have 
their usual meaning. Similarly, in the absence of con- 
tinuum errors, the error at the i'th pixel, cr(Ni), is sim- 
ply 



t7(A?,) = 3.768 X lO''* ^(5v 



(2) 



The final column density and its variance, integrated 
over the whole line (n pixels) are then 



A^: 



(3) 



For the FUSE data, cr, was taken from the error 
arrays supplied by CalFUSE, which measures fluctu- 
ations from Poisson counting statistics, and includes 
a contribution from the (modeled) background sub- 
traction. The error arrays do not take into account 
fixed-pattern noise (FPN) since no stable flat-field is 
available. Fortunately, much of the FPN in the one 
dimensional spectra is averaged out when individual 
sub-exposures are co-added ( ^2.3| l, provided that some 
drift was found in the spectra along the dispersion di- 
rection between sub-exposures. Our re-sampling of 
the data by a factor of three also ensured that the 
FPN was smeared out while preserving the shape and 
strength of an absorption line, at least when the FPN 
features were narrow. 

As we discussed, the O VI profiles towards our 
sample stars were often contaminated by CI I or H2 
lines. Our approach to reconstructing the A^(v) profiles 



17 



without these interfering Hnes was to fit theoretical 
Voigt profiles to the contaminants and remove them. 
In most cases, there was sufficient information (i.e., a 
large enough number of pixels) to define the shape of a 
line quite precisely. For this procedure, we do not need 
to know the true column density or width of a line to 
obtain a good fit; we need only assume that the LSF is 
a Gaussian. If this assumption is correct, we need only 
fit the profile and remove it to recover the true O VI 
N(v) profile. 

With the contaminating lines removed, we mea- 
sured the AOD column density as follows. For each 
sight line, we measured using the number of pix- 
els that corresponds to a velocity interval of -120 < 
Av < +120 km s"' . A measurement was also made of 
a column density limit over the same velocity range, 
defined as 
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If the measured value of was greater than a(N)T, 
we considered the line to be detected (see below). Al- 
though Av was the same for all stars, in some cases the 
A^(v) profile clearly extended beyond this range, so Av 
was increased accordingly to produce the correct value 
of A^. We again measured three times, once for the 
data normalized by the best continuum fit, and twice 
more for the upper and lower error envelopes to give 
A^,A^„ andA^,. 

In principle, we could list these AOD column den- 
sities when O VI was detected in Table |4] However, 
these quantities turned out to have very similar values 
to those measured using the line profile fitting proce- 
dure described in ^2.8. 1| We noted, however, a slight 
systematic difference in A^(0 VI), in that column den- 
sities measured with the AOD method were slightly 
over-estimated at the lowest column densities. This ef- 
fect, caused by the difference in the behavior of r for 
positive and negative fluxes, is well understood, and 
discussed, e.g., by |Fox et al.|p005] l. For this reason, 
we do not list the AOD values in Table g] The AOD 
values are given, however, in the bottom panel of the 
figures which show the FUSE spectra, as described in 
Appendix |E] 

In some cases, the values of A^(0 VI) measured 
from the AOD method can be far less than the er- 
rors, and sometimes even negative values can arise. 
For these, one has some freedom in how the numbers 
should be treated, depending on the problem at hand. 



For example, in a study that combines results from 
many sight lines, it is perhaps best to combine the for- 
mal numbers and their errors in some optimum way 
to arrive at some improved global value. In this case, 
noise uncertainties in the individual cases are allowed 
to cancel one another, which allows one to achieve an 
answer that is more meaningful than a large collection 
of upper limits. However, it may be more useful in- 
stead to give a more precise statistical interpretation of 
the outcome for an individual sight line, such as stat- 
ing to a given confidence level an upper limit for A^(0 
VI). For the latter, we propose that one could make 
use of a limit calculation for nonsignificant measures 



in the presence of Gaussian noise defined by Marshall 
1992)1. Full details of this procedure are given in Ap 



pendix|D| 

As mentioned above, in deciding whether only an 
upper limit exists for an O VI line, we adopted the fol- 
lowing procedure. If the measured value of A^ was less 
than (t{N)t, we considered O VI undetected, and in- 
stead listed the star with the measurements of A^(0 VI) 
in Table [5] Obviously, adopting a Icr cutoff between 
detection and non-detection may seem somewhat ar- 
bitrary. Another way we could have decided whether 
data yielded only an upper limit might have been to 
argue that a limit should be adopted when no Voigt 
profile could be fit to the data. In fact, these two pro- 
cedures are equivalent. There exists only one sight line 
where we could have tried to fit a profile, but where the 
measured value of A^ was less than (t{N)t- [That star 
was HD 052463 (unique ID 70), and the depression at 
the expected position of O VI has a reality that indeed 
appears ambiguous.] 

Hence Table[5]contains all the pertinent information 
available for the O VI column densities towards the 
stars (except, perhaps, forHD 052463). Columns 3-5 
give the measured values of A^(0 VI) and the values 
measured using the upper and lower continuum fits, Nu 
and Ni. Column 6 lists the error from the noise, a(N), 
while (t{N)t is the quadratic sum of all the errors, as 
defined in equation]?] Column 8 gives the 2<j limit as 
defined b y jMarshall, (1992) and discussed in detail in 
Appendix |D] Ultimately, whether one uses the AOD 
measurement of A^(0 VI) or A^umit depends on the situ- 
ation, and we discuss these issues later in the paper. 
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2.9. Misleading Effects from Unrecognized Satu- 
ration? 

A final consideration centers on the question of 
whether or not the column densities we measure for 
the FUSE stars might be underestimated due to satura- 
tion effects that are not evident in the apparent optical 
depths. Since the O VI A 1037 line is unobservable 
( ^2.8| l, it cannot be used to check for line saturation. 
If the lines are resolved, then the Voigt profile fitting 
(when a suitable LSF is convolved with the theoretical 
profile) accounts for any possible saturation. The more 
problematic situation occurs when the line is actually 
composed of many components, some of which are un- 
resolved but blended (and perhaps hidden by broader, 
resolved components). Profiles that have a mixture of 
clouds with high and low optical depths instead ap- 
pear as if they are composed of only a few clouds with 
moderate optical depths after being smeared by the in- 
strumental LSF. [For more details, see ,Jenkins| ( |1986| l.] 

There are several observations that suggest that nar- 
row unresolved O VI lines do not exist. O VI lines 
can arise from clouds that are either photoionized or 
coUisionally ionized. If the absorbing gas is in col- 
Uonal ionization equilibrium, the thermal width of the 
fine is given by = 2kT / mo, (where mo is the atomic 
mass of oxygen) which for oxygen gives T = 969b^ or, 
equivalently, b = 0.032ir'''^. The ionization fraction 
of O VI peaks at logl^iax = 5.45, which would give 
b = 16.6 km s"', and decreases sharply F] as the tem- 
perature falls, declining to a fraction 10 of the total 
oxygen at log T = 5.2 (e.g., [Shapiro & Moore|1976b[ 



[Sutherland & Dopita|[T993| [Nahar & Pradhan|[2003] l. 
This lower value equates to a Doppler parameter of 
13 km s"'. Clearly, the Doppler parameters corre- 
sponding to either the peak temperature or this lower 
value are of the same magnitude as the FUSE LSF at 
1030 A. 

It remains possible that O VI arises in collision- 
ally ionized gas that is not in equilibrium (due to rapid 
cooling, for example), or in photoionized gas. In these 
instances, the line may be narrower. We consider 
the existence of such narrow O VI components un- 
likely, for the following reasons: i) in high resolution 
(~ 3 km s"' FWHM) Interstellar Medium Absorption 
Profile Spectrograph (IMAPS) data, there is no evi- 
dence for narrow multiple sub-components in O VI 



profiles (Jenkins et al. 1998 1; ii) along path lengths 
through the entire Galactic halo, where A^(0 VI) ranges 
from ~ 1 - 4 X 10'"^ cm"^, comparisons of the AOD 
profiles for both the O VI A1032 and the A1037 lines 
show that unresolved saturation effects are relatively 
small ([W akker et al.|2003| [Savage et al.|2003"] l; Hi) no 
narrow components have been found in High Veloc- 
ity Clouds, where the high velocity of the absorbing 
gas might make such features distinct from (broader) 
absorption by the Galactic halo ( Sembach et al.|2003 



Collins et al. 2005 i; iv) similarly, in cases where an 
O VI component has a high enough velocity to be sep- 
arated from the bulk of the absorption, such as Ca- 
rina and Vela, no very narrow components are seen. 
Although the environment in which these clouds ex- 
ist (circumstellar or SNR outflows) are likely differ- 
ent from the rest of the ISM, none have widths b < 
20 km s"' ; vj the widths of the O VI lines seen towards 
stars within only 100 pc of the sun — although hard to 
measure due t o the weakness of the lines — are rarely 
< 20 km s"' (Oegerle et al.|[2005 Savage & Lehner 



2006 1. We also detect few lines with b <20 km s' 



in the FUSE sample (see ^3.13[ ). So, even when the 
number of absorbing clouds along a sight line is a 
minimum, no particularly narrow lines are seen. For 
these reasons, we conclude that although the number 
of clouds intercepted in the disk may be large, no sig- 
nificant column density errors arise from narrow (per- 
haps saturated) O VI components. 

It is interesting to note that in a recent study of ex- 



*For temperatures T which are different from Tmax, h 
16.6(r/r,^ax)i kms-'. 



tragalactic O VI absorbers, [Tripp et aL] ( [2007[ ) found 
evidence that approximately one-third of the extra- 
galactic O VI lines in low-redshift QSO spectra arise 
in relatively cool gas with logT < 5. However, those 
O VI lines are frequently located in truly intergalac- 
tic regions that have much lower densities and larger 
sizes than the gas clouds in the Milky Way disk. The 
cool, photoionized clouds found in the IGM cannot ex- 
ist in the higher density, shorter-path regions that we 
are probing in the Galactic plane. 

2.10. The Final FUSE Sample 

All the spectra belonging to the FUSE survey de- 
scribed in this paper are presented in Appendix [E] and 
Figure |24] A detailed explanation of the material 
shown for each sight line is also given. 

Although observed as part of programs PI 02 and 
P122, several stars were excluded from the final 
dataset. For completeness, these are listed in Table |6] 
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and plotted in Figure [5] Most of these sight Hnes were 
excluded because we could not choose reliable con- 
tinua. HD093204 (P10235) was rejected because no 
reliable fit to the HD line could be obtained. 

3. INTERPRETATION 

3.1. Inclusion of Other Data Sets 

The sight lines in our FUSE survey were selected to 
cover distances beyond ^ 1 kpc in order to study the 
distribution of O VI absorbing gas over a significant 
fraction of the Galactic disk. However, there are sev- 
eral other sets of O VI absorption line data that probe 
gas over different scales, from a few tens of parsecs 
in the local bubble ('LB'), to halo gas in the outer re- 
gions of the Milky Way. As we will show, we need to 
include these in our analysis to fully understand the na- 
ture of the O VI absorption. In this section, we briefly 
discuss each of the datasets used and note any modifi- 
cations we made to the published data in order to avoid 
systematic differences between quantities used for the 
FUSE sample and those used by other authors. 

3.1.1. White Dwarf Sight Lines 

Few O- and B-type stars can be found within a few 
tens of parsecs of the Sun. The nearest OB associations 
he in the Sco-Cen association, and are 118-145 pc 
away ( |de Zeeuw et al.||1999[ ). Investigations with 
FUSE therefore used nearby white dwarf stars (WDs) 
to study O VI absorption over very short path lengths. 
The first FUSE survey to study LB O VI absorption 
was completed by Oegerle et al. ( |2005|l. This was 
superseded by the larger survey of 'Savag e & Lehrier] 
( [2006, hereafter SL06), whose data we use in this pa- 
per. We adopt most of the values given in Tables 1 
and 3 of SL06, who claim that distances to the WDs 
are known to ±20-30%. In our analysis we therefore 
use a distance range defined by ± 25% of the distance. 
O VI column densities and Doppler parameters are 
taken from measurements made using line profile fit- 
ting procedures similar to our own. SL06 define O VI 
to be absent when the equivalent width Wx < 2a(Wx), 
as measured over an interval w 100 km s~' (see col- 
umn 11 of their Table 1). This is somewhat different 
than the method we adopted in ^2.8.2| but since SL06 
do not measure A^(v) in these cases we simply use col- 
umn density upper limits derived from the 2(t{W\) lim- 
its. We do not include a WD sight line if SL06 suggest 
that detected O VI may be photospheric in origin. 



3.1.2. Copemicm Sight Lines 

The Copernicus satellite was able to obtain O VI 
measurements towards some of the nearest (bright- 
est) stars lying a few hundred parsecs away, most of 
which are too bright to be re-observed with FUSE, 
unfortunately. The data were originally presented in 
[Jenkins & Meloy| ( |1974| l, ( |1974[ ), and [Jenkins 

( |1978a| l (J78)~lnthis paper, we have re-derived many 
of the stellar parameters for the Copernicus stars us- 
ing the same procedures used for the FUSE sample 
stars, so that the two samples can be properly com- 
bined. The changes made are discussed more fully 
in Appendix C.l Two stars observed by Coperni- 



cus (HDs 041161 and 186994) were re-observed with 
FUSE; these latter data supersede the Copernicus re- 
sults and the information for the stars can be found in 
Table [1] (unique IDs 52 and 27, respectively). 

O VI column densities for the Copernicus sight 
lines are taken directly from Table 1 of J78. That table 
also Hsts a value = ((y- {v) f) for the O VI A1032 
line, which is the second moment of the line profile 
after the second moment of the Copernicus LSF has 
been removed. We convert this to a Doppler parame- 
ter in the usual way, b=V2A. No errors were given 
for the derived column densities due to the difficul- 
ties in fitting the stellar continua over the short wave- 
length regions scanned by Copernicus. We have, how- 
ever, made some simple estimates of the likely errors 



in A^(0 VI); these can also be found in Appendix C. 1 



3.1.3. Sight Lines to the Vela SNR 

Our sample contains three stars which lie behind the 
Vela supernova remnant (SNR), HD 074920 (unique 
ID #77), HD 074711 (#78) and HD 075309 (#79). 
These are discussed below in |3.5[ and |3.7.2[ We 
added to these the stars discussed by Jenkins et al. 
( |1976| l and [Slavin et al.[ ( [2004 ). Distances to these stars 
were recalculated according to the prescription laid out 
in Appendix [B] 

3.1.4. Distant Halo Stars and Extragalactic Sight 
Lines 

For some of the analysis presented later, it is useful 
to compare the O VI absorption seen in the Galac- 
tic disk with that observed by FUSE over very long 
path lengths towards distant halo stars and extra- 
galactic sight lines. Zsargo et al. (2003 hereafter 



Z03) observed a small sample of stars in the Galac- 
tic halo. Three of these were retrieved by us from 
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Fig. 5. — : Spectra from programs P102 and P122 for which a reliable normalized O VI A1032 profile could not be 
analyzed. 



the FUSE MAST archive (HD 088115, HD 148422, 
and HDE 225757) before Z03's results were available, 
because they lie at latitudes \b\ < 10°. The O VI to- 
wards these sight lines is strong, and we would expect 
errors in HD subtraction and continuum placement 
to be largely insignificant (particularly since Z03 re- 
jected stars with difficult continua). In fact, we find 
that logA'^(0 VI) from our measurements agree well 



with ZOB's analysis. The difference in logA^(0 VI) is 
< 0.08 dex for all three stars, and identical within the 
la errors. 

For the rest of the stars in Z03, we have used the 
O VI column densities, Doppler parameters, and ab- 
sorption line velocities given in their Tables 4 & 5. We 
have, however, re-derived some of the distances to the 
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stars using the same procedures used for all the other 
stars described in this paper. Details are given in Ap- 
pendix [C2] 

Finally, we have included the extragalactic sight 



lines observed with FUSE by Savage et al. (2003 ) and 
[Wakker et al. (2003; ). We only use data that have S/N 
ratios > 9 per 20 km s"^ bin [i.e., those with codes 
'Q=3' or '4' in Table 2A of Savage et aT] ( [2003] )]. The 
distances to these objects can be regarded as infinite 
for the purpose of studying Galactic O VI, and no new 
distances need to be computed except for one halo star: 
PG 0832+675, for which we use c/ = 8.1 kpc. 

3.2. Distribution of Stars in the Galactic Plane 

An important parameter in the study of the O VI 
distribution in the Galaxy is the distance to any par- 
ticular background star in our sample. Unfortunately, 
for stars beyond a few hundred parsecs from the Sun, 
the only way to measure stellar distances is through 
the use of "spectroscopic parallaxes", which convert 
the observed magnitude of a star to its distance as- 
suming that the spectral type and luminosity class of 
a star are known, along with a calibration of the true 
absolute magnitude of a star of that particular type and 
class. Knowledge of the interstellar extinction along a 
given sightline, E{B-V), is also required. Although 
distances for many of our stars have been calculated 
before with this procedure, in this paper we have re- 
derived stellar distances using more modern databases. 
Since the details are complicated and not directly rele- 
vant to the analysis of interstellar O VI absorption, we 
document our procedures at the end of this paper, in 
Appendix [B] 

The distribution of our FUSE survey sight lines in 
the plane of the Galactic disk is shown in Figure [6] 
To illustrate the regions of the Galaxy in which the 
stars lie, we have plotted in gray the distribution of 
the electron density in the Milky Way, as modeled 
by Cordes & Lazio| ( [2002] l, using their NE2001 pro- 
granj^ Overplotted are the spiral arm models from the 
meta-analysis of Vallee ( 2002 1. Vallee et al. used a 
value of 7.2 kpc for the distance of the Sun to the cen- 
ter of the Galaxy; for consistency with the rest of the 
analysis in this paper, we have changed this distance to 
8.5 kpc, and altered the pitch angle of the arms to be 
11.3°. The structures shown are only meant to be illus- 
trative of where in the Galactic plane the stars lie. For 



'http://astrosuii2. astro. Cornell. edurcordes/NE2001. See also 
|2504) for a summary of NE2001. 



Cordes 



more precise details about the likely structure of gas in 
the Milky Way, see |Levine et al.| ( [2006| l and references 
therein. 

In Figure |7] we show the distribution of stars in the 
plane of the Galaxy superimposed on a map of the 
diffuse X-ray emission constructed from the ROSAT 
All Sky Survey (RASS). This map includes all the 
FUSE stars shown in Figure [6j as well as stars from 
the other datasets discussed above. The map is a 
composite of data from three ROSAT bands, the Rl 
(0.11-0.28 keV) and R2 (0.14-0.28 keV) bands, and 
the sum of the R4-R6 (0.44- 1.56 keV) bands. In 
the figure, the Rl, R2 and R4-R6 bands are assigned 
the colors red, green and blue, respectively. The pre- 
cise interpretation of the emissivity at different ener- 
gies is somewhat complicated, and is discussed in de- 
tail by'S nowden et al.] ( |1997 V However, in ^3.4| and 
3 3.5 below, we use the X-ray data to examine in more 
detail the circumstellar environments around the stars 
observed for this paper. 



3.3. Variability of O VI Absorption Lines Deter- 
mined from Multi-epoch Data 

O VI absorption from the Galactic ISM is super- 
imposed on stellar continua which may show broad 
O VI P-Cygni profiles from stellar winds. Although 
the widths of these profiles can span several thousand 
km s"' , their shapes are often highly irregular — a fact 
which likely leads to the problems we encountered in 
^2.6| in fitting many difficult continua. Moreover, the 
shapes of these wind profiles are known to vary on 
many different timescales, and can show 'discrete ab- 
sorption components' (DACs) which vary in velocity 
as well as width and strength ( York|1977t |Snow et aT] 
[T980HMassa et al.|2000||Lehner et al.|2003l ). In princi- 
ple, O VI A 1037 DACs from outflowing gas at veloc- 
ities ~ -1650 km s"' could become blended with the 
interstellar O VI A 1032 absorption studied in this pa- 
per A DAC which contaminated the interstellar O VI 
line at one epoch might only be shown to exist by tem- 
poral variations in its properties, as recorded in multi- 
epoch data. 

To understand whether these complicated continua 
variations might adversely affect accurate measure- 
ments of interstellar O VI absorption, |Lehner et al.| 
(2001) studied a small sample of stars which had been 
observed repeatedly over time intervals of either a few 
days or several months. Lehner et al. concluded that, 
at least on these timescales, stellar wind variabily had 
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Fig. 6. — : Distribution of sight lines in our FUSE survey, looking down into the Galactic plane. For illustrative 
purposes, the grayscale regions reproduce the electron density distribution modeled by'Cordes'('2004'), while the dashed 
lines show the positions of the spiral arms from models by ,Vallee ( 2002 ). [The electron densities shown are the sum 
of components measured for the inner thin disk, the s piral arms, and the local ISM; the grayscale for the spiral arms 
represent electron densities of « 0.01-0.1 cm"^ — see Cordes (20041 for more details.] There is insufficient room to 



label all the stars observed towards Carina (/ = 287.7, d = 2.5 kpc), so we indicate the region in the Sagittarius arm 
where the majority of the observed stars lie with a hght-gray circle. The gray box close to the Sun represents the Gum 
nebula; the small dark-gray circle next to it shows the position of the Vela supernova remnant (I = 263.9, d = 0.25 kpc). 
The number for each star refers to the unique ID number given in Table [2] Sight lines for which only upper limits to 
the O VI column density are found are shown as straight dashed lines. 



little influence on the measurement of O VI column 
densities. 

We can perform a similar experiment for our FUSE 
sample of stars, over slightly larger timescales. Given 



the assumptions we made in 5 2.6 for fitting the con 



tinua of stars, a comparison of O VI column densities 
measured along the same sight line at different epochs 
— particularly when the continuum is a different shape 
at those epochs — is a useful indicator of the reliability 
of our data analysis. We can also search for contam- 
ination from DACs, if only along a few specific sight 
lines. 



Six stars in our FUSE sample were observed at 
times more than a few months apart and for which 
spectra with good S/N ratios were obtained at all 
epochs; these are listed in Table |7] The four spectra 
with the highest S/N are shown in FigurejS] As the left- 
hand panels of Figure[8]show, the continua changed to 
some degree for most of the stars. To search for vari- 
ability in the O VI line, we compare the normalized 
data at each epoch in the right-hand panels of Fig- 
ure |8] By accident, several of the stars in this small 
sample had low O VI column densities compared to 
the sample as a whole. To avoid detecting differences 
in absorption line profiles arising from difficulties in 
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Fig. 7. — : O VI results superposed on a depiction 
of the diffuse X-ray emission recorded by the ROSAT 
All Sky Survey (RASS) over the longitude range a) 
£ = 310° to 50°, b) £ = 40° to 140°, c) £= 130° to 
230°, and d) £ = 220° to 320°. Ellipses are centered 
on the positions of target probes; their vertical axes 
are proportional to A^(0 VI) (totaled over all veloc- 
ity components), and their horizontal axes signify the 
distance to the targets, both of which can be gauged 
by the scales in the lower left corner of the figure. A 
circle corresponds to an average sight-line density of 
1.62 X 10"^ cm-3 [i.e., A^(0 VI) = 5 x 10^^ cm'^ over 
a distance of Ikpc], which is a representative aver- 
age for the whole survey. It follows that tall, skinny 
ellipses indicate more O VI than usual for a particu- 
lar distance, while the converse is true for short, fat 
ellipses. O VI absorption towards stars at distances 
< 100 pc are not plotted since the size of their ellipses 
in both the horizontal and vertical direction would be 
too small to be seen using these scaling rules. The 
direction of a small arrow on the top of each ellipse 
signifies the average velocity of all of O VI absorp- 
tion, with an angular deflection to the left (negative 
velocities) or right (positive velocities) whose mag- 
nitude in degrees corresponds to the radial LSR ve- 
locity in kms~'. Cases with only upper limits are 
shown as dotted ellipses, where the A'(0 VI) upper 
limit defines the height of the ellipse in the same man- 
ner as described above. White ellipses mark stars in 
the FUSE survey reported in this paper, red ellipses 



denote Copernicus measurements Jenkins & Meloy| 
[19741 |^rfe|[T974l |Jenkins|[T978a| ) and light blue el- 
lipses identify targets in a FUSE survey of stars in 
the Galactic halo (Zsargo et al.'2003). Yellow ellipses 
show the data from .Jenkins et al., (1976) and |Slavin| 



et al. ( 2004 1 for the Vela supernova remnant. Green 



arrows represent O VI absorption towards nearby 
WDs from the data of |Savage & Lehner| ( [2006| l. Ver- 
tical dark blue bars at high Galactic latitude show 
the results from the FUSE survey of extragalactic 



objects by Savage et al. (2003 1 and Wakker et al. 
( 2003| l, and have no width because they he outside 
the Galaxy. Colors in the X-ray sky are keyed to the 
ROSAT Rl (0.11-0.28 keV) band (red), R2 (0.14- 
0.28 keV) band (green), and the sum of the R4-R6 
(0.44- 1 .56 keV) bands (blue), with scale factors that 
yielded a good differentiation of colors across differ- 
ent regions. The RASS results were obtained from an 
ftp download site http://www.xray.mpe.mpg.de/rosat/ 
survey/sxrb/12/fits.html maintained by the Max- 
Planck-Institut fiir Extraterrestrische Physik. 
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Fig. 7. — : continued 
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subtracting the HD R(0) line, we show only the O VI 
profile before the removal of the HD line. (To see the 
strength of O VI after the HD removal, see Fig. [24]) 

Figure |8] shows that there is little evidence for any 
variation in the O VI(+HD) absorption. Any differ- 
ences seen can be attributed to problems in fitting the 
continuum. The absorption towards HD 103779 is a 
particularly interesting example: the continuum is very 
hard to fit at both observed epochs, yet the absorption 
at 1032 A is remarkably similar in observations made 
four years apart. Note that the line towards HD 000108 
is entirely HD absorption (again, see Fig. [24] ID #43). 

Although we find no evidence for variations in ab- 
sorption over differences of 3-4 years, the spectra 
do allow us to test how robust our measurements of 
A^(0 VI) really are. Although the O VI hne in Fig- 
ure [8] is blended with HD absorption, we can proceed 
with the analysis described in the previous sections 
by removing the HD line and measuring the resulting 
A^(0 VI). The results are given in the last column of 
Table |7] Despite having to contend with difficult con- 
tinua (e.g. HD 103779), and removing HD lines that 
are a significant fraction of the O VI column density 
[because A^(0 VI) is relatively weak], the column den- 
sities measured for each epoch agree to within their 
errors. (The upper limits for the spectra of HD 000108 
observed a year apart are obviously more a reflection 
of the S/N of the data, which were similar due to the 
similar exposure times. For both spectra, however, the 
HD could be modeled and removed to show that no 
O VI was obviously present.) This implies that our 
measurements of A^(0 VI) are quite insensitive to tem- 
poral variations, continuum changes, and HD removal. 

3.4. Accounting for Circumstellar Bubbles 

The target stars and their association members 



produce stron g stellar winds (|Garmany et al. 1981 
Lamers 1981 Garmany & Conti 1984| l which inject 



substantial amounts of energy into their surrounding 
gaseous media ( |Abbott|198"2{ |McKee|1986| . From an 
initial theoretical investigation of the probable con- 



sequences of such wind-driven shells. Castor et al. 



( 19751 proposed that nearly all of the low-velocity 
O VI features seen in the spectra of early-type stars 
recorded by Copernicus and reported by |Jenkins &| 
Meloy ( 1974 1 probably arose from a thin conduc- 



terstellar medium. In a more refined investigation of 
these circumstellar structures, [Weaver et al.| ( |1977[ ) 
predicted that a typical bubble interface could produce 
A^(0 VI) « 2 X 10" cm-2. While this amount of O VI 
is small compared to the typical column densities reg- 
istered in our FUSE survey, we must recognize that our 
inventory of truly interstellar O VI could be elevated 
by the contributions from circumstellar bubbles. 

Many factors can influence the structure of a bub- 
ble, such as the star's age, its wind energy, its motion 
through the medium, and the density of the surround- 
ing gas. These elements, together with the added ef- 
fects from neighboring stars, make it difficult to make 
trustworthy predictions about the O VI contributions 
for specific cases. Nevertheless, we know that the in- 
ternal hot gas within a bubble should emit soft X-rays. 
Thus, stars with well established bubbles should be de- 
tectable in the RASS, although target regions for which 
there is strong foreground absorption by neutral hydro- 
gen may show up in only the highest energy bands. 
While we do not expect a simple one-to-one corre- 
spondence between the X-ray emitting properties of a 
bubble and its ability to contribute O VI absorption, we 
can still try to establish an empirical relation between 
the two and thus gain some insight on the probable ef- 
fects from bubbles in general. 

Figure |9] shows selected RASS image field^ that 
include groups of stars in our survey. These exam- 
ples illustrate the fact that some stars are surrounded 
by clearly defined X-ray emitting zones, while others 
are not. Sometimes there are several stars that reside 
within a common envelope that is bright in X-rays. 

Ultimately, we need to determine the approximate 
magnitude of the bubble contributions to the O VI col- 
umn densities that we measured, so that we can ar- 
rive at a clearer understanding of how much O VI 
is distributed in the general regions of the Galactic 
disk. However the distinction between "bubble" and 
"non-bubble" O Vl-bearing gas is somewhat arbitrary, 
since there is probably a continuum of bubble prop- 
erties that extend from well formed, identifiable bub- 
bles with distinct boundaries to highly disrupted struc- 
tures that begin to blend into the general medium. De- 
spite this limitation, we felt it was important to make a 



tion zone situated between an interior region of hot, 
shocked wind material at T ^ 10'' K and the sur- 
rounding medium, rather than from the general in- 



The ROSAT images were obtained from the ROSAT Data Browser 
at the web site http://www.xray.mpe.mpg.de/cgi-bin/rosat/dataJ 
browser using the default settings except for the full-scale limits 
being set to Band 1 = 5, Band 2 = 6, Band 3 = 4 counts deg"^ s"' 
for Figs.[9^-d and f, while Band 1 = 5, Band 2 = 12, Band 3 = 8 was 
implemented for Fig.|9fe. 
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Fig. 8. — : Four stars (of six) for which multi-epoch data are available. Black/blue spectra represent the earlier/later 
epoch data. The left hand panel shows the variation in absolute flux at two different epochs, along with the continua 
used to normalize the data. Fluxes are in the usual units of ergs cm"^ s"' A"'. The right-hand panel shows the 
normalized flux before subtraction of the HD 6-0 R(0) line, while the dates of the two epochs used are given in the 
top left of each panel. Small shifts (< 0.06 A) in the wavelength scale have been added where necessary to align the 
features. 



distinction between stars that had visible bubbles, ac- 
cording to some specific criterion, and those that did 
not. To accomplish this, we assigned each target to 
one of three categories according to appearances in the 
RASS fields, based on the concepts of having no ev- 
idence of surrounding emission, having an intermedi- 
ate or doubtful level of emission, or having a clear, 
surrounding bubble. For convenience, we assign the 
symbol R to represent a category. We then have: 



the sky, well free of any emission; 

• R=\: the target is in or near a general diffuse 
enhancement of X-ray emission or is on the edge 
of a bright spot; 

• R = 2: the target is clearly situated in projection 
in a bright X-ray spot that appears to have been 
created by the star and/or its association neigh- 
bors 



R = Q: the target is in an X-ray-dark portion of 



(An additional category will be assigned in § 3.5 ) 
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Fig. 9. — : X-ray emission observed in the ROSAT All Sky Survey (RASS) within fields centered on close groups 
of stars selected as targets for the O VI survey. Intensities for the low, middle and high energy bands are depicted 
with varying brightness levels of red, green and blue, respectively. The positions of O VI survey targets are identified; 
following each star's name are parentheses that contain our designation of ROSAT category followed by the value for 
A^(0 VlVlOi"* cm-2. 
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The immediate surroundings of virtually all early- 
type stars exhibit some X-ray emission arising from 
shocked stellar wind s (Hillier et al.||1993t |Berghoe 
fer & Schmi"ttl|T994| [Cassinelli et al. 1994; Berghoe^ 
fer et al.|1996 1 1997V Absorption features from O VI 
should arise from such material, but the large wind 
speeds allow these contributions to be separated from 
the interstellar features ^Lehner et al.|2001|[2003l l. For 
this reason, we refrained from classifying stars asR = 2 
if the emission was not clearly extended in the sky be- 
yond the smearing by ROSAT point-spread function. 
As an aid to determining whether or not the emis- 
sion was extended, we examined entries in the ROSAT 
bright and faint source catalogs (Voges et al. 1999| l. 



Sources coincident with the star positions that were 
not clearly extended (extent likelihood < 1) were con- 
sidered to arise from fast winds, which we considered 
grounds for rejecting the assignment of R = 2. 

3.5. Accounting for Supernova Remnants 

In addition to bubbles around the target stars, ex- 
traordinary enhancements of O VI can in principle 
arise if any of the lines of sight penetrate a supernova 
remnant. This is clearly the case for 3 of our stars, 
HDs 074920, 0747 11 and 075309, which lie in the di- 
rection of the conspicuous Vela SNR [see Fig. [Tj^j]. 
The near edge of this remnant is at a distance of only 
about 250 pc ^ Cha et al.|1999[ ), a value that is consis- 
tent with the slightly greater distance 280 pc) to the 
nebula's center as determined from the distance to the 



Vela pulsar ( Ogelman et al. 1989 Caraveo et al. 2001 
podson et al. 20031 and X-ray observations of the neb 



ula (Bocchino et al .|1999 1. Thus, there is no question 
that the remnant is in front of these stars with distances 
ranging from 1.1 to 2.7 kpc. 

It is important to identify less obvious cases where 
sight lines may penetrate SNRs. A comparison of our 
target locations with those in a compilation of identi- 
fied remnants ( |Green|2001| l reveals a number of coin- 
cidences, which are hsted in Table |8] Once again, in 
order to exclude special cases from our general assay 
of O VI, we have added one more category to the ones 
defined in 33.41 above: 



• R = 3: stars are situated inside the boundary of 
a supernova remnant and are not clearly in the 
foreground. 

To be conservative, when the relative distances are 
in doubt, we assign a Category 3. Unlike the situa- 



tion for the Vela SNR, the distances to the remnants 
presented in Table [8] may have large errors. For this 
reason, we sometimes experienced difficulties in de- 
termining with much confidence whether or not the 
targets are behind the SNRs. Two stars that fall into 
this category ai-e HD 168080 and HD 047417. For 
HD 124314, we could find no distance estimate for 
the remnant in the literature. For other stars, such as 
HD 185418 and HD 199579, the differences in dis- 
tance between the stars and the remnants are so large 
that we can be more certain that they are indeed in the 
foreground. Without the knowledge that HD 185418 
was by chance situated in front of G53.6 -2.2 (with 
angular dimensions of only 33 x 28'), we would have 
assigned a rating of R = 2 for this star (i.e., by mis- 
taking the X-ray emission as coming from a bubble 
surrounding the star) instead of our adopted R = 0. Al- 
though one might have said the same for HD 199579, 
we have retained the R = 2 assignment because there is 
a bright X-ray spot centered at the position the star that 
appears to be on top of the more extended emission 
from the SNR. While it may seem that HD 13441 1 Hes 
far behind its remnant, we must recognize that the dis- 
tance estimate for this SNR is probably much less ac- 
curate than the others, since it is based on the empirical 
(z-adjusted) surface-brightness — distance (E-D) rela- 
tion for estimating remnant distances (Mi lne|1979| l. 

A very young remnant G266.2 -1.2 has a diame- 
ter of 2° and is located in projection inside the bound- 
ary of the much larger Vela SNR ( |Aschenbach„ 1 998| 
lyudin et al.|1998[ ). [Redman et aL] ( 2002^ suggest that 



a nearby nebula RCW 37 (the Pencil Nebula) is cre- 
ated by a collision between a collimated flow of gas 
ejected from G266.2 -1.2 and a boundary of the Vela 
SNR. If this interpretation is correct, the distance of 
G266.2 -1.2 should be nearly the same as that of the 
Vela SNR. If not, we can still rely on an upper limit 
of 0.5 kpc based on the present-day strength of ga mma 



ray emission from "''*Ti ( Aschenbach et al.|l999 ). The 



two stars listed in Table |8] have total column densities 
(in units of lO'"* cm'^) of 1.96 (HD 074920) and 3.37 



(HD 075309), along with 3.0 for HD 075821 (Jenkins 



et al. 1976| l, which is also inside G266.2 -1.2. By com- 



parison, stars that are still inside the main Vela remnant 
but outside G266.2 -1.2 include HD 074711 in our 
survey along with HDs 074455, 072108 and 074753 
in the study by Jenkins et al. ( 1976| l. These stars have 
A^(OVI)/10i4 cm"^ = 0.7, 4.5, 1.3, and 0.8, respec- 
tively. It is thus not clear that the smaller remnant con- 
tributes much extra O VI. Perhaps most of the gas in- 
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side this young remnant is too hot to show O VI; its 
X-ray spectrum indicates t hat gases exist at a temper a- 
ture in excess of 3 x 10^ K ( Aschenbach et al. 1999 k 



There are a number of targets in the vicinity of the 
Cygnus Superbubble, a horseshoe-shaped X-ray emit- 



ting region with a diameter of about 13° ( jCash et al. 
|1980[ ) (see Fig. [tJjJ). Recent investigations indicate 
that this may not be a single structure, but rather a 
superposition of unrelated shell and bubble-like ob- 
jects at different distances ( Bochkarev & Sitnik|1985 
Uyamker et aL][200r 



We therefore refrain from as- 
signing /? = 3 to stars that appear near the edge of this 
emission. 

The stars HDs 045314, 047417, 047088, and 



047360 overlap with the Monogem Ring ( ,Plucinsky 



etal.|1996| ), but Figure |9}:j shows that HD 046202 and 



HD 046150 are just outside the edge of this remnant 
(visible in the top left corner). The yellow color of 
the Monogem Ring depicted in Figure [7|:j indicates 
the presence of a strong low-energy X-ray flux that 
has not been appreciably absorbed by foreground neu- 
tral hydrogen [i.e., compare with the blue color of the 
Cygnus Superbu bble in Fig.|7|7j], and indeed [Plucin-| 



sky et al. ( 1996 1 estimate that only 5 x 10 



cm ^ of 



neutral hydrogen is in front of the nebula. Our stars 
in this vicinity all have A^(H I) > 10^' cm"^ (Diplas & 



Savage|1994) , so they are clearly at greater distances. 
As with the stars behind the Vela SNR, these stars 
show more O VI than expected for their distances. 

Stars in the general vicinity of £ = 240° are be- 
hind a lai-ge superbubble GSH 238h-00h-09 identified 
by |Heiles ( 1998) l, who estimated 0.8kpc for a distance 
to this structure. Figure |7]:/j shows no appreciable de- 
partures from the general average A^(0 VI) per unit dis- 
tance for the stars behind this superbubble, although all 
of the O VI velocities are consistently positive. Per- 
haps this velocity pattern arises from O VI in the back 
edge of an expanding bubble. 

3.6. ROSAT Categories For Other Datasets 

For completeness, we also assigned ROSAT cat- 
egories to the stars observed with Copernicus, and 



those observed by |Zsarg6 et al. (2003 1, using the same 
system we used for the FUSE stars described above. 
These values are given in Tables [12] and [13] in Ap- 



pendix C.l and C.2 



ROSAT categories were also derived for the WDs 
of SL06. The majority of these stars were assigned to 
R = 1 or R = 2 categories — there were very few WDs 



with no detectable (R = 0) X-ray emission, largely due 
to their small distances. The WDs give us an interest- 
ing opportunity to look for a correlation between X-ray 
flux and A^(0 VI). Since WDs are very close by, inter- 
stellar O VI absorption is likely to be much less impor- 
tant than it is over longer path lengths, and circumstel- 
lar O VI absorption might in principle be the dominant 
source of O VI absorption. Using WDs also provides 
the advantage of using sources with uniform proper- 
ties. Although not shown here, we searched for a cor- 
relation between X-ray flux and A^(0 VI), but found 
none. The circumstellar environment of a WD may, 
of course, be different from that of a stellar associa- 
tion, but for the WDs, the detected O VI absorption is 
unrelated to the WD X-ray flux. 

3.7. Average O VI Line of Siglit Volume Density 

5. 7. 7. Sample Definitions and Basic Statistics 



Savage et al. ( 2003 1 demonstrated that O VI absorb- 



ing gas can be traced several kpc above the plane of 
the Milky Way. In the Galactic halo, the volume den- 
sity falls with height, but the derivation of the mid- 
plane density no requires fitting an assumed model to 
the observed data (see next section, §3.8| l. The stars 
in our sample (excluding the halo stars of Z03) have 
a distance above the Galactic plane z = dsinb with a 
mearj^ distributed around z = and with 98% of the 
stars below |z| = 1 kpc. We would therefore expect a 
simple average of the densities n = N{0 Nl)ld to be 
close to the mid-plane density since the scale height is 



~ 2.3 kpc according to Savage et al. (2003 1. (Again, 
we discuss this assumption in ^3.8[ ) With the sight 
lines coded by the presence or absence of X-ray emis- 
sion, it should also be possible to see if n is differ- 
ent for lines of sight with and without enhanced X-ray 
emission. 

The first problem in determining n involves decid- 
ing how to treat sight lines when only upper limits are 
available for A^(0 VI). We initially employed various 
survival statistics, including the well known Kaplan- 
Meier (K-M) estimator, which calculates the true 
distribution function of randomly censored data and 
yields the distribution's mean and median ( |Feigelson] 
|& Nelson|[T985) |Isobe et al.|[T986l ). Given measured 
values of, and upper limits to, A^(0 VI)/t/, we should 
in principal have been able to recover (n) . However, 



as Isobe ( 1992 1 makes clear, the K-M technique only 



'The precise mean is (z) = -44 pc, witli a dispersion cr = 98 pc. 
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works when the upper Umits are distributed randomly 
with respect to the overall distribution ("Type 2 cen- 
soring")- Although the method worked well for SL06, 
who had upper limits with values comparable to the 
distribution of detected values of n towards WDs, all 
of our upper limits for n were clustered below the ap- 
parent mean of n. 

Instead, we simply defined a sample for which the 
column densities were set at the upper limits given 
in Table |5] or as stated in the other surveys when no 
Une was present. We discuss the possible errors intro- 
duced by defining the sample with these upper limits 
in the next section. Table |9] summarizes the statisti- 
cal tests employed, and the results obtained. We cal- 
culated two statistics: (1) the average value of n de- 
rived from the ensemble of log(n) value^ ; and (2) 
the distance weighted average ^A^(0 Vl)/^(i. The 
data were further subdivided into two subsets: stars at 
distances > 2 kpc and those with distances between 0.2 
and 2 kpc. These distance categories partition the data 
into two halves with approximately equal numbers (the 
median distance to the stars is 2040 pc). Of course, the 
two subsets contain data from different surveys in dis- 
proportionate numbers: the d = 0.2-2 kpc subset con- 
tains a mixture of Copernicus and FUSE data (45 and 
43 sight lines, respectively) while the d > 2 kpc con- 
tains almost exclusively FUSE data (97 of 102 sight 
lines). Although the Copernicus sight lines are essen- 
tial for mapping O VI absorption at distances less than 

are 



1 kpc, the errors in A^(0 VI) (Appendix C.l 



very much larger than the FUSE N{0 VI) errors. The 
distance-selected subsets have the advantage, however, 
of subdividing the data according to whether or not the 
sight lines are strongly influenced by O VI from the 
Local Bubble. Finally, within each of these categories, 
the samples were subdivided further into subsets ac- 
cording to their ROSAT classifications. 

3. 7.2. The O VI Volume Density n: ( i) Differences for 
Near and Distant stars. 

Examples of the distributions of n values are shown 
in Figure 10 In each case, the distribution of n for 
sight lines regardless of their ROSAT classification 
("All-T?") is given in gray, while R = sight lines are 
shown as white histograms. All the histograms are 



normalized by the total number of stars used to create 
the distribution. The top panels include stars at "far" 
distances d > 2 kpc; the bottom panels show stars at 
"near" distances of 0.2 < d < 2 kpc. The mean of 
a distributiorj^ is shown by an inverted triangle, and 
a Gaussian profile for the All-/? distribution using this 
mean and the standard deviation given in Table |9] is 
shown with a dotted line. A cross marks the distance- 
weighted average log[^A^(0 VI)/^ii]. As Figure [To| 
shows, the histograms are reasonably well approxi- 
mated by log-normal distributions, enabling us to de- 
fine means and standard deviations consistent with a 
theoretical normal distribution. 

In parallel with these measurements, we con- 
structed a synthetic dataset of absorption line profiles 
in order to better understand the distributions shown 
in Figure [TO] and the values tabulated in Table |9] We 
were particularly interested in how well we might re- 
produce the observed distributions if we based a sim- 
ulation on a single value of n, and hence whether the 
dispersion in the data might really be some artifact 
of the quality of our data. We therefore took a set of 
random distances, used a fixed value of n to derive 
A^(0 VI), created line profiles (adopting a fixed value 
of ^7 = 40 km s"'), added a fixed amount of noise 
(similar to that obtained for the majority of the FUSE 
sight lines), and re-measured A^(0 VI) and the appro- 
priate errors. We also adopted the same procedures 
for defining an upper limit to the column density as 
described in ^2.8.2 when lines were barely detected, 
and included the upper limits to n determined from 
the upper limits to A^(0 VI). Finally, we calculated the 
'measured' values of n and determined how close they 
were to the original value. 

The simulation provided several interesting results. 
How well we could reproduce (n) (using the same 
methods described in the previous section and shown 
in Table |9]l depended primarily on the S/N we adopted 
when adding noise to the theoretical profiles. With a 
high S/N, all lines were recovered and n was indis- 
tinguishable from the original value. With very low 
S/N data, many lines had measured column densities 
which were either negative, or were less than the upper 
limit determined from the noise. As a consequence, 
the mean value of n became larger than the original 
value (since an upper limit is larger than the true orig- 



'"This is a simple average, without weiglits tliat account for tlie dif- 
ferent magnitudes of the en'ors. The analysis presented in j3.8.2| 
evaluates minima in which automatically recognizes differences 
in the size of the errors. 



"jc = S(logn,)/A', where N is the total number of i points used. 
The standard deviation entries in Table [9] are then the usual = 
Ti(\ogni—xf- /N; they are not the uncertainties of the listed aver- 
ages. 
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Fig. 10. — : Left: Distribution of O VI line of sight volume densities, log n = log[A^(0 YI)/d], divided into two sets, 
one for stars beyond 2 kpc (upper panel), the other for stars in the interval 0.2 < d <2 kpc (lower panel). The gray 
histogram shows the distribution of stars regardless of ROSAT class, while the white histogram shows only R = sight 
lines. The dotted line shows a Gaussian with a mean and standard deviation for the All-/? samples given in Table [9] 
while the inverted triangles show the means of each distribution. A cross marks the distance weighted average of 
A^(0 VI). (More details are given in Table[9]) Note: distant high-latitude halo stars, extragalactic sight lines, and sight 
lines towards the Vela SNR and Carina Nebula are not included. Each histogram is normalized by the total number 
of objects making up the histogram. Right: The same distributions, only with the Local Bubble (LB) column density 
[A^lb(0 VI)= 1.11 X 10'^ cm"^] removed for each sight line, and distances reduced by the radius of the LB, which we 
take to be 100 pc. 



inal value). 

Translating these findings to our FUSE and Coper- 
nicus data sets is not so straightforward, since the real 
data have differing S/Ns (more precisely, the spectra 
all have different A^(0 VI) upper limits). However, our 
simulations show that while the error in the true mean 
of n depends on the S/N of the data, it is equivalently 
dependent on the relative number of lines which are 
lost (lower S/N means more n upper limits). Hence it is 
possible to simply count the number of non-detections 
(upper limits) we find in our data and use the results of 
the simulations as a guide to the likely over-estimation 



of («>. 

In the near sample (0.2 < d < 2.0 kpc) 14% of our 
sight lines show no O VI absorption lines and have up- 
per limits for A^(0 VI). For the far sample (d > 2 kpc) 
only 7% have upper limits. In our simulation, we find 
that these percentages of non-detections cause an over- 
estimate of ^ 12% and ~ 6% for the near and far sam- 
ples, respectively. Hence, i/the underlying O VI col- 
umn density can be represented by a constant value, 
then the numbers shown in Table [9] are likely to be 
over-estimated by only a small amount after including 
column density upper limits. 
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Another result from the simulation was that the 
largest values of n (in cm"^) recovered arose from the 
weakest A^(0 VI) lines, and hence towards the closest 
stars. This was because weaker lines were more easily 
lost in the noise, and were more likely to give erro- 
neously large values of n, just because of noise fluctu- 
ation^Given this possibility, we considered whether 
large values of n measured in our data were actually 
real, and not just a consequence of selecting lower 
quality data. In fact, values of log(n) > -7.5 are found 
at all distances and column densities and so are real 
over-densities in the global distribution of O VI ab- 
sorption. Moreover, seven sight lines have log(«) > 
-7.0, all of which arise towards the Vela SNR. We also 
noticed that most of the sight lines towards the Carina 
Nebula also had large values of n. 

Sight lines through the Carina Nebula (NGC 3372, 
[Ub] = [287.7,-0.8]), of course, are not typical probes 
of the Galactic ISM. The nebula is a giant H II region 
(1-2 degrees in radius in optical emission) powered 
by more than a dozen star clusters, and although most 
of the gas has been blown away by the clusters, low- 
and intermediate-mass star formation remains at the 
peripheries (Smith et al.|2003) . The distance to Carina 
is somewhat ambiguous, but appears to be in the range 
2.5-3 kpc ( |Feinstein|[T995| |Walborn|[T995] ). Most of 
the stellar distances we derive to stars towards Carina 
(Appendix [B] and Table [T} are consistent with this dis- 
tance, as expected. The extent of the nebular in X-ray 
emission can be seen in Figure |9^, which also shows 
the positions of the stars observed by FUSE. 

It has been known for a long time that the spec- 
tra of Carina stars show complex, multicomponent 
absorption line profiles, with individual components 
distributed over several hundred km s"' at both posi- 
tive and negative velocities ( jWalborn & Hesser|1975 



19821. Absorption from ions in both high and low 



ionization states is seen ( |Walborn et al.|[T984l |2002a| 
[2007, and references therein). Most of the stars in our 
FUSE sample in this region are of very early spectral 
type and have reasonably well determined continua, 
which makes the detection of interstellar O VI straight- 
forward. They generally show clear evidence for ab- 
sorption at large negative velocities away from the bulk 
of the absorption at v ~ 0, high velocity features which 



'^These inherently weaker lines were also more likely to give erro- 
neously small or negative values of n as well; however, since weak 
lines deemed below the noise were set to upper limits, the distribu- 
tion of logn values discussed above does not include these small- 
log(;i) outliers. 



are not seen along other sight lines through the ISM 
(except towards the Vela SNR — see below). In some 
cases, a component can be clearly resolved from the 
main complex; in other cases, the additional negative 
velocity gas produces a discernable asymmetry to the 
V ~ complex. It seems reasonable to link this ad- 
ditional absorption to gas directly associated with the 
nebula. Detecting high positive velocity gas in absorp- 
tion is somewhat harder, since the H2 6-0 R(4) line 
hes at a wavelength that corresponds to any O VI ab- 
sorption offset from rest by -1-121 km s"^ potentially 
masking any high positive velocity components. Nev- 
ertheless, it is clear that there is no preponderance of 
very strong O VI absorption lines at « 100 km s"', 
since the profiles of the 6-0 R(4) lines are always close 
to the shape predicted by fitting other H2 lines in dif- 
ferent regions of the spectra (see §2.6i. Only towards 



HD 093146 (#92) is there good evidence for excess 
O VI absorption at positive velocities. 

As noted in ^ the Vela SNR ([l,b]= [263.9,-3.3]) 
is much closer than Carina, at a distance of only 
250 ±30 pc (|Cha et al.|[T999l). Sight lines which 



pass through Vela again show complicated multi- 
component absorption line complexes (see, e.g. |Slavin| 
et al.|200"4 and references, therein). An initial study 



of O VI absorption was made using four stars situated 



behind Vela (Slavin et al. 2004); our FUSE sample 
includes an additional three stars which lie behind 
the SNR: HD 074920, HD 075309, and HD 074711 
(#77-79). The first two of these show evidence for 
multicomponent absorption. HD 075309 in particu- 
lar appears unusual, with very weak O VI absorption 
at V0 = 33 km s"', and much stronger absorption at 
-74 km s~' . Again, the kinematic structure of the O VI 
absorption seen towards HD 074920 and HD 075309 
seems quite atypical of that seen in the rest of the 
Galactic ISM. 

A full analysis of the conditions within the Carina 
Nebula and Vela SNR infeiTed from absorption line 
data is beyond the scope of this paper Nevertheless, 
the sight lines towards these regions present us with 
a dilemma. Although they present a particularly in- 
teresting opportunity to study phenomena local to ac- 
tive star formation and destruction, they also present 
particularly extreme examples of O VI absorption and 
probably do not represent the more general distribu- 
tion of O VI absorbing gas in the ISM. For this reason, 
we excluded Carina and Vela sight lines from all the 
measurements presented in Table |9] 

We make one final note concerning the possible 
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effect that the non-detections of O VI might have on 
the values of (n) given in Table |9] For the distance- 
weighted statistic Y.N(0 Vl)/Y.d (column 4 of Ta- 
ble [9]l, it is also possible to set all upper limits to zero 
and see whether (n) is very different than that calcu- 
lated using the upper limit value^ When calculated 
this way, we find that (n) is only 7% and 14% smaller 
for the near and far ("All-7?" — see below) sample, 
respectively. 

IShelton & Cox| ( |1994| ) first showed that the LB has 
a higher average density than the rest of the Galactic 
disk. (We see this effect in Fig. [15] in p?T4| ) To 
better reflect the value of («) in the ISM beyond the 
LB, we recalculated all the statistics after subtracting 
a column density representative of that from the LB, 
for which we used a recent value of the LB volume 
density of O VI found by SL06, «lb = 3.6 x lO"** cm'^ 
For a bubble with a radius of 100 pc, this corresponds 
to a column density A^lb(0 VI) = 1.11 x 10'^ cm"^. 
We also subtracted 100 pc from all stellar distances. 
Values of (n) corrected in this way are given in Table|9j 
these are labelled as "-LB" measurements, compared 
to the values for which no LB contribution is removed, 
"H-LB". 

The corrected distributions of log(«) are shown in 



the right-hand stack-plot of Figure 10 and tabulated in 
Table |9] For the distant d > 2 kpc stars, it is unsur- 
prising that there is little difference (0.05 dex) between 
the H-LB and -LB means. For the nearer (0.2 < d < 
2.0 kpc) stars, the difference is far more pronounced 
(0.13-0.18 dex), as would be expected. Moreover, 
the near and far H-LB samples (for the same ROSAT 
classes) also show large differences (0.09-0.14 dex); 
yet with the contribution from the LB removed, the 
means between the near and far samples are nearly 
identical. We conclude that after correction for the Lo- 
cal Bubble, the average line of sight density of O VI is 
the same for all distances beyond 200 pc. Note that 
in the bottom-right panel, the width of the near, -LB 
sample is much wider than in the other three panels. 
This is because of the increase in the relative errors in 
A^(0 VI) — and hence log(«) — that arise after sub- 
tracting Ni^BiO VI). This can also be seen in the stan- 
dard deviations listed in column 3 of Table |9] which 
are much larger for the near, -LB sample. 



3.7.3. The O VI Volume Density n: (ii) Differences 
for Different ROSAT Classes. 

We repeated the statistical tests described above, for 
stars separated by ROSAT class. The results can be 
seen in Table |9] where average densities are broken 
into categories defined by samples of sight lines which 
include "All-/?", or which have R = OotR>0. Table|9] 
shows that within each of the four near/far, +LB/-LB 
classes, the differences in the means between the R = 
and the R > sub-samples are between 0.05 and 0. 1 
dex. Are such differences significant? 

To test if these differences are real, we used two 
well established tests on the +LB sample of stars. 
We first considered the results from a Student's f-test, 
which tests the null hypothesis that the means of two 
populations are equal. For the d > 2 kpc stars, we 
found a f -statistic off = -1 .6, and a probability that we 
might incorrectly reject the hypothesis that the means 
are equal of only P = 0.1 1. For the 0.2 < d < 2 kpc 
stars, we found t = -0.8, and P = 0.44. We would con- 
clude from these numbers that there is a small differ- 
ence in n between R = Q and R > Q sight lines beyond 
2 kpc, but that no significant difference can be found 
nearer than 2 kpc. 

Student's f-test is a parametric test whose accu- 
racy relies on the fact that the two populations being 
tested are normally distributed. To avoid problems in 
comparing two distributions when the underlying dis- 
tributions are not known, we used the standard non- 
parametric Kolmogorov-Smirnov (K-S) test. For the 
distant stars, the K-S D parameter was measured to be 
D = 0.3, with P = 0.014, while for the nearer stars, 
we found D = 0.13 and P = 0.84. These numbers 
confirmed the conclusion that the R = and R > 
have different distributions of log(n) for distant stars, 
but very similar means for the nearby stars. The num- 
bers from the K-S tests and from the Student's f-test 
changed little if we used the -LB instead of the +LB 
sample. 



'^Obviously, log[iV(0 VI)]/rf cannot be calculated if ^(O VI) is set to 
zero, which is why we could not calculate (n) from a distribution 
using such a sample of column densities. 



Figure 1 1 shows the distribution of log(n) values 
for the two +LB ROSAT classes with the two distance 
subsets combined, i.e., for all stars beyond 200 pc. 
For the R = and the R > samples, the means 
and standard deviations are (log(n)) = -7.85 ±0.32 
and -7.77 ±0.31, respectively. A f-test shows that the 
likelihood of a difference between the two sub-classes 
has increased by combining the near and far samples: 
t = -1.9, f = 0.07. A K-S test gives D = 0.20, 
P = 0.03. As we mentioned above, these tests were 



34 



distances > 200 pc 
0.15 - □ -R = 

: n R>o 

Gnnssinn profiles: 
0.10 - <log n>=-7.S5±0.32 
<log n>=-7.77±0.31 

0.05 - r- 



0.00 [ 



-9.0 



-8.5 -8.0 
log [n (cm"')] 



-7.5 



-7.0 



Fig. 11. — : Distribution of O VI line of sight volume 
densities for all stars in the Galactic disk (thereby ex- 
cluding halo stars) beyond 200 pc, divided into R = Q 
(white histogram) and R > Q (gray histogram) subsets 
for the +LB samples. Each histogram is normalized by 
the total number of objects making up the histogram. 
Gaussian profiles for the R = (dotted line) and R > 
(dashed line) samples are drawn with means and stan- 



dard deviations given in 5 3.7.3 



performed with Vela and Carina sight lines excluded. 
If we include those sight lines, most of which have 
^ > 0, the two distributions become very different 
(t = -3.6, P = 3x 10-4; D = 0.26, P = 1 x 10"^) be- 
cause the lines of sight clearly intercept regions with 
higher densities. Again, these values change little if 
we use the -LB instead of the +LB sample. 

We discussed the likely origins for higher (on av- 



erage) values of n towards R > sight lines in 5 3.4 



and ^3.5| Bubbles visible in X-rays and the volumes 
within SNRs probably contain some gas in the tem- 
perature range that favors the production of O VI, and 
although the effect is small, sight lines towards these 
structures show elevated O VI column densities. 

We note a few interesting exceptions, however 
HDs 199579 and 110432 have R = 2, and yet they 
both show virtually no O VI absorption at our level 
of sensitivity. A more spectacular example of this ef- 
fect is 7 Cas observed by J78. The upper limit of 
A^(0 VI) < 10'^ cm-2 is considerably below the FUSE 
detection limit, and yet this star is centered on an ex- 
tended source that is 19" in diameter, according to 
the ROSAT Bright Source Catalog ( jVoges et al.fT999 1. 



However 7 Cas is only at a distance of 0.2 kpc, which 
means that the detected X-rays are emitted within 
0.02 pc of the star. This dimension is far smaller than 
a classical, well developed bubble. Thus, the emis- 
sion could arise from the more extended parts of the 
hypersonic wind from the star, rather than an estab- 
lished interface region. It is also likely that if this star 
was ten times further away, at the typical FUSE sam- 
ple distances, the X-ray emission would probably be 
undetected. 

The reason for being able to see the differences be- 
tween the R = and the R > sight lines to distant 
stars and not nearby ones is likely due to selection ef- 
fects. X-ray fluxes fall with distance as l/d^, and neu- 
tral gas absorbs X-rays, so for the more distant stars 
we probably only see the very largest, brightest inter- 
stellar bubbles which have the largest enhancements 
of A^(0 VI). Nearby stars are less likely to be members 
of dense star-forming regions, and their X-ray fluxes 
are probably easier to detect; so, for these, we assign 
aR > category even if the circumstellar bubbles are 
not particularly bright and there is little enhancement 
in N{0 VI) from the bubble. 

Finally, we note one other experiment we per- 
formed with the data. We attempted to subtract a 
column density that might be representative of cir- 
cumstellar O VI absorption, A^ciic(0 VI). Our aim was 
to find a value of A^circ(0 VI) which could be sub- 
tracted from the R = 2 sample, leaving a sample with 
a mean which was the same as the R = Q sample. Un- 
fortunately, we were unable to find a plausible value 
of A^circ(0 VI); the distribution of n for the d > 2 kpc 
sample was largely unaffected by subtracting a small 
A'circ(0 VI) from each sight line, but it greatly changed 
the 0.2 < d < 2 kpc sample, sending the majority 
of the n values to values less than zero. This proba- 
bly reflects the selection effect mentioned above: the 
value of A^circ(0 VI) towards the more distant stars 
is probably higher than A^circ(0 VI) towards nearby 
stars, because our R> Q categories include very large, 
very X-ray-luminous bubbles at large distances, and 
smaller, less luminous (but stiU detectable by ROSAT) 
bubbles at smaller distances. 

3. 7.4. Bias Toward Stars with Small Reddening 

In large part, our selection of targets is governed 
by the limitation that the stars must not be so faint 
at 1032 A that a prohibitively large amount of observ- 
ing time is needed to obtain spectra with respectable 
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values for the S/N. One consequence of this is that 
more distant stars tend to have much lower than nor- 
mal amounts of reddening per unit distance. For closer 
stars, this selection is not as strong. It is fair to 
ask whether or not this bias in the sampling of the 
Galactic disk volume distorts our conclusions on the 
general distribution and properties of the regions that 
hold O VI, and in particular, on the volume densities 
we have derived above. This could be an important 
consideration if O Vl-bearing gas is strongly corre- 
lated or anticorrelated with the cold gas that causes 
reddening. For instance, surveys of the soft X-ray 
background radiation exhibit an anticorrelation, which 
probably arises from the effect that the X-ray emit- 
ting gas displaces the cold material ( |McCammon & 
|Sanders|1990l l. 

To examine the issue of whether or not we must 
acknowledge and attempt to correct for a bias in our 
survey, we must determine the likelihood that the char- 
acter of the O VI, either its average volume density or 
velocity dispersion, is affected somehow by the red- 



dening per unit distance. Figure 12 shows the distri- 
bution of all stars in our FUSE survey on a diagram of 
E{B-V) vs. distance, with the points having a color 
coding that represents either the O VI volume density 
n (top panel) or the second moment of the absorption 
profile (v^) (bottom panel). The expected selection ef- 
fect in our survey is clearly evident from the fact that 
there are no points in the upper right-hand portion of 
the diagrams. Nevertheless, from the distribution of 
colors shown in the plots, we gain the impression that 
internally within our collection of targets, there is no 
tendency for the results to be driven much by the red- 
dening per unit distance. We can easily imagine that if 
the diagram were populated more evenly, i.e., without 
the exclusion of cases where color excesses increase 
roughly in proportion to distance, the outcome would 
not be appreciably different from what we have already 
found. 

Looking at the problem more quantitatively, we 
find a median reddening per unit distance of 0.10 
mag kpc~', represented by the dotted line in each panel 
of the figura^ For the half of our measurements 
that are above this value, the median value of n is 
1.64 X 10"** cm""*, while for the half of the sample 
below the median E(B - V)/d, the median value for 



'''Our median reddening per unit distance is considerably lower than 
the general average of 0.61 mag kpc"' iSpitzer 1978, Fitzgerald| 
|1968^ . Whether or not this is important depends on the volume oc- 
cupation fractions of gases at different densities. 
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Fig. 12. — : The distribution of targets with regard 
to their color excesses E(B-V) and distances d, ex- 
cluding targets towards the Carina Nebula or the Vela 
supernova remnant (which are atypical). In each panel, 
the dotted line indicates the median value of E(B- 
V)/d = 0.10 mag kpc"' in our survey, and the dashed 
line represents the general average of 0.61 mag kpc~' 
in our part of the Galaxy. Points with circles around 
them are targets that have ROSAT classifications of 
R = 0, and the crosses represent cases where only up- 
per limits for A^(0 VI) could be obtained. The colors 
in the top panel indicate n(0 VI), while those in the 
bottom panel indicate the second moment of the ab- 
sorption profile (v^). 



n is 1.30 X 10"^ cm""*. This outcome is opposite to 
our expectation that regions with less reddening than 
usual might be more likely to hold more than the usual 
amount of O VI. However, an apparent correspondence 
between ^(B-V)/!^ andA^(0 VI)/c/ could arise purely 
from the fact that errors in a common denominator, i.e. 
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the distances, might make the quantities appear to be 
mildly correlated when, in fact, they are not. A K-S 
test on whether the two separate samples of n could 
arise from the same population indicate that there is 
an 8% probability that differences this large or larger 
could arise simply from chance, given that the two 
samples are drawn from the same population. (Such a 
test cannot be run for the velocity dispersions because 
they are indeterminate when O VI is not detected.) 

Finally, we note that even if the selection of sight 
lines is biased against highly reddened sight lines, the 
clouds which cause the reddening are likely to occupy 
only a small volume compared to the volume that the 
O Vl-absorbing gas occupies. If so, then the value of 
n we calculate along a line of sight will be unaffected 
by whatever the extinction happens to be. 

3.8. O VI Mid-Plane Volume Density 

The correlation between A^(0 VI)sin|Z7| and the 
height z above the plane of the Galaxy has been de- 
termined before ( |Jenkins||I978b[ [Savage et alT][200? 



However, the increased size of our sample of disk 
stars, and a well determined column density for the 
LB, makes it worthwhile revisiting the problem. 

The O VI volume density as a function of z is usu- 
ally defined such that n = nof{\z\). The column density 
for any individual sight line at a height |z| = t/sin|/7| 
above or below the Galactic plane is then 

N(OWl) = [ n(r)dr= /' n(|z|)csc(|fe|)t/(|z|) 
Jo Jo 

r\zsu,\ 

= noCsci\b\) f(\z\)d(\z\), (5) 
Jq 

where r is the distance along a sight line. Convention- 
ally, it has been assumed that the O VI volume density, 
n, falls off exponentially with height above the disk, so 
that /(|z|) = e"l^l/* which means that 



n = n^e 



-\z\lh 



(6) 



where h is the scale height of the O VI absorbing gas. 
Equations [5]and[6]then lead us to an expression for the 
column density projected on a line perpendicular to the 
Galactic plane. 



A?(0 VI)sin|/7| =«(,/!(! -e"!'!/'') 



(7) 



include are those from our own FUSE survey, data 
from Copernicus, the WD sample of SL06, the halo 
stars observed by Z03, and the extragalactic sight lines 
measured by Savage et al. ( 2003| i and Wakker et al. 
( |2003) ). However, including the latter dataset forces 
us to consider an important problem uncovered by 
these authors, which we must discuss before we at- 
tempt to derive «o and li, namely the difference be- 
tween sight lines in the northern and southern Galac- 
tic hemispheres. Further, we must also introduce the 
problem that, although the density of O VI does de- 
crease in the way suggested by equation]?] the layer of 
O VI absorbing gas in the disk is not smoothly strat- 
ified, but instead is rather clumpy. We discuss these 
two problems below. 

3.8.1. The North-South Divide 



Obviously, to measure «o and h, we can fit the func- 
tion given in equation |7] to the data. The data we 



In their survey of extragalactic sight lines, |Savage| 
etal](j2003) reported an excess of 0.20— 0.30 dex in the 
average value of log[A^(0 WYjsinb] for the north Galac- 
tic Polar region towards sight lines at b > 45° com- 
pared to lower latitude directions in the north and to the 
entire southern Galactic sky. To investigate the asym- 
metry in more detail, we examined whether the extra- 
galactic O VI column densities measured at latitudes 
b > — the 'N' or 'North' sample — were statisti- 
cally different from those measured atb <Q — the 'S' 
or 'South' sample. In all cases, column densities were 
reduced by a factor sin \b\. We first considered the re- 
sults from a Student's f-test; for the N and S samples of 
log[A^(0 WYjsinb], we found a f-statistic of f = 4.1, and 
a probability that we might incorrectly reject the hy- 
pothesis that the means are equal of only P = 1 x 10"^. 
A K-S test yields D = 0.42, P = 4x lO"! Such re- 
sults clearly bear out the assertion of Savage et al. that 
a significant difference between A^(0 VI) exists in the 
northern and southern Galactic hemispheres towards 
extragalactic sight lines. 

As described in ^3.1.4[ we decided to use only 
the higher-quality data from [Savage et al.| ( |2003| l and 
Wakker et al. ( 2003[ l (Q > 2) in our analysis. If we re- 
peat these two tests for this sub-sample alone, the dif- 
ference between the N & S samples is smaller, but still 
significant. Student's f-test gives t = 2.4 and P = 0.02, 
while the K-S test gives D = 0.39 and P = 0.05. 

Despite these variations, which occur in the Milky 
Way halo, we see no difference in our disk stars for 
similar N and S samples. For stars beyond 200 pc, 
there is no clear difference between the two samples 
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when measured by a f-test (f = 0.5, P = 0.62) or a 
K-S test (D = 0.11, P = 0.62). This conclusion holds 
whether or not a correction is made to the observed 
column densities for the Local Bubble. 

For this reason, in fitting the theoretical relationship 
of A^(0 VI) sin I I with |z| given in equation |7] to the 
data, we use all our sight lines regardless of whether 
they lie in the northern or southern Galactic plane, ex- 
cept for the extragalactic sight lines: we derive no and 
h twice, once each using the northern and southern ex- 
tragalactic sight lines. 

3.8.2. dumpiness in the O VI Absorbing Gas and 
Errors Used for N( O VI) sin \b\ 

To measure no and h, we fitted the function given 
in equation |7] to the data by varying no and h and min- 
imizing the statistic between the theoretical func- 
tion and the data. The first obvious result we found 
was that although the data were clearly correlated in 
the way described by equation [vj the reduced was 
much greater than the expected value of unit^p^ 

At face value, we might say that large values of 
xt arise either because the assumed errors in the data 
are inadequate or that the underlying model must be 
wrong. We discuss the second of these conclusions 
below. As far as the first conclusion is concerned, we 
have taken special care to derive realistic errors for the 
quantities measured. A similar problem was found by 
[Savage et aL] ( |1990| l for the distribution of Al 111 in 
the Galactic plane. Their solution was to introduce an 
additional "error", which takes into account the patchi- 
ness or dumpiness of the absorbing gas that is ignored 
by using equation]?] Although the density of O VI fol- 
lows, primarily, an exponential decrease with height 
above the plane, the distribution is not smooth and pla- 
nar, but is, at some level, clumpy. Adding an additional 
"error", ccl, allows for deviations caused by this addi- 
tional dumpiness when seeking an acceptable value of 
X^. The tJcL term is applied assuming that the quan- 
tity is independent of distance, an assumption which 
we demonstrate to be true in ^3.9| 

Given the introduction of acL, we can combine the 
natural deviations attributed to this term with the ex- 
perimental errors in determining A^(0 Vl)sin|/7|. Or- 
dinarily, this operation would be problematic because 
variations of the parameter of a fit would respond dif- 



ferently to the two different kinds of errors. In our 
case, we are fortunate to have most of our measure- 
ments made at heights |z| < h. When this is true, the 
slope of the dependent variable x in the fitting function 
in relation to the data points y is fixed to unity. With 
this constraint, errors in x are equivalent to errors in y. 
For each data point, therefore, we describe a total error 
exclusively in terms of an error in just the y-direction 
equal to a value 



.-r^ _ I I 



(8) 



These individual errors are all relative errors: ap- 
plies to the measurement error in the column densitjPj 
cr[A^(0 VI)]/A^(0 VI); Ud is the relative error in the dis- 
tance, aid) I d; and acL is the additional fractional er- 
ror we derive by requiring xt to be unity when min- 
imizing the fit of A^(0 VI)sin|/7| with |z| to the data. 
The absolute error in each point is then simply 

cr [A^(0 VI) sin |fe|] = N(OVl)sm\b\xaT. (9) 

It should be understood that ctcl is introduced not 
simply to justify a model whose functional represen- 
tation is wrong. In principle, we could propose some 
other trend of decreasing A^(0 VI) with distance from 
the Galactic plane, but we would stiU require the use of 
(TCL to account for an intrinsic scatter in the observed 
O VI column densities that exceeds the measurement 



3.8.3. The Mid-Plane Density and Scale Height of 
O VI 

A plot of log[A^(0 VI) sin \bW with log(|z|) is shown 
in the top panel of Figure [13] where we distinguish 
between R = Q (black points) and R > Q stars (gray 
points). In ^3.7.2| we found evidence that the density 
of the ISM beyond a few hundred pc was better repre- 
sented with the removal of a contribution from the LB. 
Hence, to derive no and h, we used the R = Q sample 
of stars, and subtracted a column density A^lb(0 VI) = 
1.1 1 X lO'^ cm~^ for the LB along each sight line. Ac- 
cordingly, we also subtracted the LB radius — 100 pc 
— from all distances. Sight lines with only upper lim- 
its to A^(0 VI) have no at^ errors, yet values are needed 



'^Here, we use the conventional definition of = wliere u 

is the number of data points used in the fit minus the number of 
variables used, in this case 2. 



'*As discussed in 5 2.8.1 the column density errors are asymmetric 
due to the different continuum errors. Here, for on, we simply take 
the average of the two column density errors. As we will show, 
however, (jn "C cFd, so the exact errors in the column densities have 
little effect on the derivation of and h. 
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for calculating y^. For these points, we used the me- 
dian (Tai (0.08 dex) for the sample with detected O VI. 

The fit of A^(0 VI) sin \h\ with |z| must be made in 
log-log-space, for the following reasons. In a ^ fit, 
the errors in each point determine the weight that that 
point gives to the value of x^. Unfortunately, the ab- 
solute error in the column density is proportion al to 
A^(0 VlQ and because (as we will show in § |3.9| ) 
A^(0 VI) is proportional to the distance to a star, the ab- 
solute error in the column density is also proportional 
to distance. Similarly, aid) is more or less a fixed con- 
stant of « 20-30 % of (see Appendix [B]i, so the ab- 
solute error in d also increases with d. The result is that 
the final absolute error is proportional to distance, so 
that more distant stars would have had lower weights 
in the fit had we chosen the linear representation. If, 
instead, we fit log[A^(0 VI)sin|/7|] against log|z| (i.e., 
a fit in log-log-space), then the error is the same rela- 
tive error as that given in equation [8] and is no longer 
proportional to the distance. 

With the data divided into two samples containing 
separate northern and southern extragalactic sight lines 
( ^3.8. l| l we found the following: 



N : no 

h 

S : no 
h 



1.33 X 10"** cm"^ 
4.6kpc, (TcL = 0.25 dex 

1.34 X 10"** cm-^ 
3.2kpc, (TcL = 0.28 dex 



Unsurprisingly, the value of no does not depend on 
whether the N or S samples are used, since no depends 
primarily on the data below |z| ~ 1 kpc. Conversely, 
the values of h do depend on which extragalactic sight 
Unes are used. As described in ^3.8.2| the values of 
<JCL are set to give x^ = 1 ■ We note that these values 
of log no = -7.89 are almost identical to the centers 
of the Gaussian fits made to the distribution of n val- 
ues discussed in § |3. 7. 2| (column 5 of Table [9] for the 
= 0, -LB sample). 

The results from this fit are shown in the bottom 
panel of Figure [13] (where we only plot -LB, = 
data points). Note that we only show the original col- 
umn density errors, and not the values of uj used for 



'^This is because strong lines are wider at the continuum, and have 
higher continuum errors. Since a(N) in the FUSE data is domi- 
nated by continuum errors, cr{N) and N are roughly con'elated, even 
though there is no correlation between distance and the quality of 
the data. 



minimizing x^- In fact, the derived value of ctcl dom- 
inates the errors, with (Jn << (Jd < <Jcl- For example, 
increasing the errors in the distances to the stars by as 
much as 1.5 times their assumed value has no effect on 
the values of no and h; it is the inherent dumpiness of 
the interstellar medium that dominates the dispersion 



in A^(0 VI) sin 1^7 1 seen in Figure 13 



The new values of h given above differ from the 
single scale height calculated by Savage et al. ( 2003| l, 
h = (log[A?(0 VI)sin|fe|])/no = 2.3 kpc. Most of 
the difference in h for the southern hemisphere comes 
from the reduction in no, for which Savage et al. used 
a preliminary value of 1.7 x 10 



(e.g. 



Jenkins 



2002), a value which included all sight lines through 



the Galactic disk (including those toward Vela and Ca- 
rina as well as the R > Q sight lines) and for which no 
correction was made for the LB. 

Savage et al. noted that A^(0 VI)sin \b\ towards ex- 
tragalactic objects was enhanced by ^ 0.25 dex in the 
high latitude directions 45° < b < 90° compared to 
lower latitudes in the north and compared to all south- 
ern hemisphere directions. Rather than reporting two 
scale heights as we have, they instead described the 
distribution of O VI in the halo with a single scale 
height combined with a high latitude northern hemi- 
sphere 0.25 dex enhancement. (Savage et al. consid- 
ered a number of possible explanations for the origin 
of the enhancement.) With the smaller mid-plane den- 
sity found in this paper, the single scale height used 
in Savage et al.'s model would increase from 2.3 to 
2.9 kpc. 

What are the likely errors in these values of no and 
hi It is straightforward to calculate x^ for values of 
no and h away from those that gave the minimum x^- 
Figure 14 shows contours of x^ for different values of 



no and h, with each contour level drawn at a particular 
limiting value of x^ = xl - The value of x? is calcu- 
lated as Xl = Xmm~^Xp(^) (Lampton et al. 1976 1. Here, 
Xp(Q;) is the value of x^ which would be exceeded a 
fraction a of the time in a set of random trials, where 
the degrees of freedom is set to the number of free pa- 
rameters p that we allow to vary (2 in our case). Given 
some contour xl set by adopting a significance a, we 
would expect the true value of no and h to lie within the 
Xl contour in a 1 - a fraction of similar O VI surveys. 

we adopt values of xi for confidence lev- 
0.68, 0.95 and 0.997, which were 



14 



In Figure 
els I- a of 0.38 



chosen because they resemble the conventional inter- 
pretation of the 0.5, 1, 2 and 3a confidence levels that 
apply to a normal distribution. It should be understood. 



39 



Fig. 13. — : Top: Plot of O VI column density re- 
duced by the sine of the Galactic latitude b, against 
the height of a star above (or below) the plane of the 
Milky Way, \z\ = dsin\b\. Point shapes, explained 
by the key in the top left corner, refer to datasets ob- 
tained by various authors, and include the measure- 
ments made in this paper ("FUSE"), those obtained 
by Copernicus (Jenkins'lQTSa ), local WDs discussed 
by Savage & Lehner (2006), sight lines towards the 



Vela SNR (Slavin et al. 



20041 |Jenkins"eraL][T976ll, 



those towards halo stars gs argo et al.|2003jl, and the 
extragalactic sight lines observed by [Wakker et al.| 
( |2003] l and |Savage e"taL] ( [2003] l. Points colored black 
refer to stars which have ROSAT class R = 0, while 
gray points have R > Q (see 5 3.4 and 5 3.5i. The 



extragalactic sight lines are shown as crosses at an 
arbitrarily large distance; they have have no ROSAT 
classification, and are separated into sight lines above 
(Northern, or 'N' sight lines, left group of crosses) 
and below (Southern, or 'S' sight lines, right group 
of crosses) the Galactic plane. Their distances are 
slightly offset from each other so that the errors in the 
column densities can be seen. Two-sigma upper lim- 
its are shown by downward-pointing arrows attached 
to open symbols. (The type of symbol again indi- 
cates which dataset the limits are taken from.) The 
size of a point indicates the relative uncertainty in 
the distance to the star: a key at the lower right indi- 
cates how the size is related to the uncertainty (large 
points represent smaller uncertainties so that they are 
more influential in guiding the eye). An absolute er- 
ror of ±30 % in any distance is shown below this 
key. Bottom: The same plot, only with a contribution 
from the Local Bubble subtracted from every sight 
line [Ni^b(0 VI)= 1.11 x 10'^ cm'^], and a distance 
of 100 pc removed from every sight line. Only R = 
stars are plotted. The solid line shows the best fit 
to the data assuming that the O VI volume density n 
varies as n = no^^"'''/''' where no is the mid-plane den- 
sity and h is the scale height. Although the value of 
riQ remains unchanged (because all stellar sight lines 
were used in fitting n to the data), the value of h is 
different depending on whether N or S extragalactic 
sight lines are used. The two curves show the fit for 
h = 4.6 kpc (N) and h = 3.2 kpc (S). 
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Fig. 14. — : Contour plots of about the minimum 
found for fitting log[A^(0 VI)sin|Z7|] against log|z|. 
Contour levels are drawn at levels of which cor- 
respond to a significance of a = 0.62, 0.32, 0.046 and 
0.0027, or confidence levels l-a = 0.38, 0.68, 0.95 
and 0.997. Although the distribution function of is 
not described by a normal distribution, we use the anal- 
ogy to describe the contours as 0.5, 1, 2 and 3 sigma 
confidence levels. 



of course, that the distribution function for with a 
small number of free parameters is not normally dis- 
tributed. 

The values of ctcl given above seem to suggest 
that the dumpiness of the O VI absorbing gas may 
be slightly different between the northern and southern 
hemispheres. Are these differences real? The problem 
for estimating the errors in ctcl are somewhat differ- 
ent than deriving the confidence levels for no and h 
discussed above. (Tcl is not a variable in the function 
used to measure x^, so does not enter into the mini- 
mization of x^- Changing ctcl changes the numerical 
value of x^ when the minimum is found, but the val- 
ues of no and h do not chang^^ However, we can ask 
how much we can vary ctcl before we get 'unreason- 
able' values of Xmin- Given the number of degrees of 
freedom (number of data points - 2) for our data, we 
know how Xmin -should be distributed. We can there- 
fore select the two values of within which the y^. 
outcomes are expected by chance 68% of the time (i.e.. 



a ±1(7 probability interval), and then consider the lim- 
iting values of ctcl that make our evaluations of Xmm 
conform to these two x^ limits. 

For the values of no and h given above, we find that 
for values of x^ within ±1(T of the expected value, 
we obtain ctcl = 0.25(+0.02,-0.03) dex for the north- 
ern sample, and ctcl = 0.28(+0.02,-0.03) dex for the 
southern sample. Since the two values of <tcl are 
within ± Icr of each other for the N and S sight lines, 
it seems that there is no significant difference in ctcl 
in the northern or southern hemispheres. If we sim- 
ply average these two values then, we have acL — 
0.26 ±0.02 dex. 

We conclude this section by noting that in the sim- 
ulations of |de Avillez & Breitschwerdt] (|2005bj) this 
mid-plane density has a somewhat different meaning 
than the single value which equation|6]implies exists at 
|z| = 0. In their models, hot gas is generated by SN ex- 
plosions throughout the plane of the Galaxy, and O VI 
arises in clumps, vortices and filaments which arise 
from the mixing of hot (interior) and cold (exterior) 
gas associated with the SNRs. Inside these clumps, the 
O VI volume density can vary by as much as 6 dex, al- 
though typical clouds have densities of log n ^ -9 to 
-6. Although there is, therefore, no 'single value' of 
n(0 VI) in the disk of the Milky Way in these mod- 
els, the regions are sufficiently small (~ 100-200 pc) 
that over long path lengths, there is, obviously, an av- 
erage value of n. |de Avillez & Breitschwerdt] ( |2005bp 
find that the average is constant over time, and has a 
value of 1.8 X 10"^ cm"^ (during a time period of 100 
to 400 Myr and assuming solar metallicities). This 
number is surprisingly close to the values we list in 
Table |9] for c/ < 2 kpc stars when we make no dis- 
tinction in ROSAT class and subtract no contribution 
from the LB, as is appropriate for comparison with the 
simulations. 

3.9. Distribution of O VI Column Density with Ef- 
fective Distance 



Following Jenkins ( 1978b I we calculate a reduced 
distance to a star to take into account the decline in n 
with |z| discussed above: 



J, = /!(l-e-|'l/'') csc\b\ 



(10) 



'^This is not entirely true of course. With, for example, ctcl set to 
zero, the values of — and hence the weights for the fit — are 
different for each point, which leads to small differences in the min- 
imum value of riQ and h. These differences are very small, however. 



In fact, these distances are nearly identical to the regu- 
lar distances d for the stars studied herein — except for 
halo stars, whose long path lengths and high distances 
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above and below the plane of the Galaxy require cor- 
recting for the change in n with \z\. 



The top panel of Figure 15 shows the correlation 
of A'^CO VI) and dg, along with the expected value 
A^(0 VI) = node using the value of «o found above in 



5 3.8.3 We would expect this line to pass through the 
majority of the R = points, and this is the case for 
distances beyond a few hundred pc. Below these dis- 
tances, however, there is a clear departure from the lin- 
ear relationship. This apparent discrepancy can be re- 
moved by plotting the column density for each point 
reduced by the LB contribution. This is shown in the 
bottom panel of Figure [15] The excess of points above 
the node line below a few hundred pc now disappear 
below the y-axis limit shown in the figure. 

The fact that A^(0 VI) and d^ are correlated has 
some important consequences. First, it means that the 
O VI absorbing gas is interstellar, and cannot be pri- 
marily circumstellar in origin. Second, if the O VI 
lines arise in discrete structures, and not in a uniform 
smooth plane of gas (see below), then the lines must 
be composed of many individual components which 
are so close in velocity to each other that they blend 
together to form a 'single' line with no apparent ve- 
locity structure. Third, because the correlation arises 
using sight lines in all four quadrants of the Milky Way 
(Fig.[6]l, then the physical structures of the ISM which 
give rise to O VI absorption must be ubiquitous over 
all the Galaxy. 



Figure 16 shows how n(0 VI) (with a contribution 



from the LB removed) for stars at distances d > 1 kpc 
varies with Galactic longitude. Ignoring sighhtlines to- 
wards the Vela SNR and the Carina Nebular, we find 
little difference in n in whichever direction we look, 
with one exception: of the 23 FUSE sight lines in the 
region / = 80° - 140° and with d > I kpc, ten have only 
upper limits to A^(0 VI). This direction corresponds 
to the outer Perseus Arm (see Fig. |6]l. One could ar- 
gue that there exists a genuine deficiency of O VI in 
the anti-center direction. However, more than half the 
sight lines have «(0 VI) commensurate with the rest 
of the Galactic disk. We discuss below whether the 
scatter in n with Galactic longitude is significant. 

Figure[T5]shows that the dispersion of column den- 
sities around the best-fit line at a given distance is 
large. This in itself suggests that the O VI absorb- 
ing ISM is far from being a smoothly distributed in- 
tercloud medium — otherwise the points would more 
closely follow a A^(0 VI) = no d^ relationship. The fig- 
ure also indicates that absorption does not simply arise 



from randomly distributed uniform clouds. If absorp- 
tion arose in an ensemble of generic clouds, each with 
a fixed column density A^o, then the total column den- 
sity measured by intercepting p clouds would be pNo- 
For a random distribution of clouds the dispersion in 
the total column density should go as ±^/pNo, and the 
fractional error in A^(0 VI) would fall as 1 /^/p. Hence 
the dispersion in the observed column densities should 
decrease at larger distances, which is clearly not ob- 
served. 

To demonstrate this, we plot dashed lines in the 
bottom panel of Figure [15] to show how the error 
in A^(0 VI) would decrease with distance, assuming 
No = 5x 10'^ cm"-, which is typically the smallest col- 
umn density detected towards WDs. [Smaller values 
of A^o produce no dispersion in A^(0 VI) at any dis- 
tances, while larger values exceed the smallest column 
densities measured.] This suggests that the O VI ab- 
sorbing regions have deviations over a wide range of 
A'o (or scale sizes) compared to fluctuations caused by 
randomly situated clouds with a single value of N[). 

It follows that the largest of these clouds have a 
low space density, so that they do not appreciably in- 
crease the overall average density n^, but instead re- 
veal themselves by creating large relative deviations in 
A^(0 VI) for long sight lines. This is done in a man- 
ner that compensates for the expected 1 / ^/p decrease. 
Alternatively, the maintenance of a constant level of 
relative deviations in average space density at large 
distances could arise from a large-scale clustering of 
small clouds. 

Another way to measure the variation of dumpi- 
ness with distance is to repeat our fit of A^(0 VI) sin \b\ 
against |z| for sets of stars at different distances, and 
examine how the values of ucl change. We can, once 
again, divide our stars into two sets, those at dis- 
tances between 0.2-2 kpc and those beyond 2 kpc; 
we assume that there is no difference between sight 
lines in the northern and southern Galactic hemi- 
spheres, and that we can therefore use all our stel- 
lar sight lines. If we use column densities with a 
contribution from the LB removed, and all ROSAT 
classes, we find that at distances d > 2 kpc, (Tcl = 
0.28(+0.03,-0.02) dex, while for 0.2 < d < 2 kpc, 
<7ch = 0.29(+0.04,-0.03) dex. These results imply 
that (Tcl is indeed independent of distance. 

Could these similarities simply arise from the fact 
that the volumes sampled by the more distant stars 
have greater separations than those toward the nearby 
targets? That is, could large-scale regional differences 
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10 



16 



Fig. 15. — : Top: Plot of O VI column density 
against the effective distance d,, to a star, given by 
de = /![l-e(-''^'"l''l/'')] csc|Z7|, where h is the scale 
height, b is the Galactic latitude, and c/ sin j/jj = |z| is 
the height above (or below) the plane of the Milky 
Way. The values of h used are the same as those 



shown in Fig. 13 but since most of the stars are in 
the plane of the Galaxy, dg and d are virtually in- 
distinguishable, except for the halo stars. The size, 
shape and shading of the symbols used to plot the 
data are the same as those given in Fig. 13 The 



solid line shows the predicted value of the O VI col- 
umn density if A^(0 VI) = no d^, using the value of 
no shown in the bottom plot. (Note, this value of 
riQ is the value derived for = sight lines, which 
in this plot are shown by black points. Hence the 
fit should pass through the body of these data.) A 
difference of ±20 % in «o is shown by dotted lines. 
Bottom: The same plot, but with a contribution from 
the Local Bubble subtracted from every sight line 
[A?lb(0 Vl)= 1.11 X 10'^^ cm-2]. The two dashed lines 
demonstrate how the dispersion in n(0 VI) should 
change if A^(0 VI) simply increases by intercepting 
more clouds of a fi xed column density (in this case 
5 X 10'2 cm-2; see |l9i. 
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in the Galaxy, which are not evident in the short sight 
Unes, create variations that almost exactly compensate 
for the expected decrease in the dispersion in log« for 
a random distribution of uniform clouds? We can gain 
an insight on the possible importance of this effect by 
performing an analysis of variance (e.g. |Wolf|[T974] ). 
We divide all 56 cases for R = Q and d > 2 kpc into 



k = 7 Galactic longitude bin^ having indices /, with 
m = 8 samples (with indices j) per bin. We performed 
a Kruskal-Wallis //-test to see if the deviations of the 
results from one bin to the next are significant, and 
we found that they are to a 96% level of confidence 



'The bin boundaries were set to the following longitudes (in degrees): 
28, 73, 104, 210, 244, 292, and 314. The widths of the bins differ 
because of our requirement to have an equal number of samples in 
each bin. As before, we exclude sight lines to Carina and Vela. 
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Fig. 16. — : Average O VI volume density n [= A^(0 Vl)/c/] along sight lines as a function of Galactic longitude. Only 
stars more distant than 1 kpc are included. Sight lines toward the Vela SNR are indicated by star symbols. More distant 
stars are shown with smaller symbols, as indicated in the legend. Squares represent stars observed as part of our FUSE 
survey, while triangles represent stars observed with Copernicus (J78). Squares are colored black for R = sight lines 
and gray for 7J > sightlines. Upper limits to n are shown by open symbols with downward pointing arrows; again, 
symbols representing upper limits towards stars with R = Q have black outlines, while those representing upper limits 
towards stars with R > have gray outhnes. The extent of the Perseus Arm in the outer Galaxy is indicated (see 
Fig.|6}. 



(H = 13.3). Given that these deviations are probably 
real, we now must understand quantitatively what frac- 
tion of the overall sum of the squares of the deviations 
from the mean value of log n 



k m 



5o = 2^ 2_J}ogni,j - (log n) tot) 
,=1 j=\ 



(11) 



arises from deviations from one longitude group to the 



next, 

k 

5i=m^((log«),-(logn)tot)' (12) 

as opposed to the cumulative effect of deviations 
within each group, 

k m 

S2 = Y,Y.^\ogni,j-{\ogn)if (13) 
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where (log«)tot is the overall average logn and (logn), 
is the average logn within each group /. It can be 
shown that, in general, Sq = Si +82- For the samples 
defined above, we found the values. 



Sq = 5.22, Si = 1.30, and ^2 = 3.92 



(14) 



which shows us that only 25% of the variance of the 
entire sample arises from deviations from one longi- 
tude set to the next. We conclude that (Tcl is invari- 
ant with distance and that for the most part this ef- 
fect arises from the basic texture of the volumes that 
hold O VI absorbing gas, while a smaller contribution 
arises from large scale regional variations in the Galac- 
tic plane. 

Finally, we note that if the contribution from the LB 
is not removed, the values of ctcl change slightly, but 
there is still no difference between nearby stars and 
the stars further away: acL = 0.25(+0.03,-0.02) dex 
for > 2 kpc, and ctcl = 0.24(+0.03,-0.02) dex for 
0.2 < <i < 2 kpc. 

3.10. Variation of O VI volume density with Milky 
Way Spiral Arms 

On the premise that most of the O VI is gener- 
ated by blast waves from supernovae, we might ex- 
pect to find evidence in our data that the number den- 
sity of such events and the character of their devel- 
opment could be strongly influenced by whether their 
surrounding volumes are located within spiral arms 
or interarm regions. In order to explore this issue, 
we attempted to detect an arm vs. interarm contrast 
in the averages for the O VI volume density n. We 
defined spiral arms in terms of their electron densi- 
ties «e sensed by pulsar dispersion measures, which 
ultimately were assembled into the model created by 
|Cordes & Lazio| ( |2002| l shown in Figure [6] Our test 
was a simple one: we attempted to find a correlation 
between slight-line values of (n) and {n^), where the 
latter was expressed in terms of the model's disper- 
sion measure along the sight line divided by its length. 
We performed this comparison for all stars that had a 
ROSAT class R = and distance d > 0.2kpc, again ex- 
cluding Vela and Carina stars. The Pearson correlation 
coefficient between log(n) and log(ne) was 0.08 for 98 
sight fines, which means that no significant correlation 
could be found. (We recognize that errors in distance 
can induce a spurious correlation, but this is not an 
issue here since we are not claiming to see a correla- 
tion.) In short, we see no clear evidence that variations 



in sight-line values of (n) are strongly influenced by 
the presence or lack of spiral arms. 

3.11. Variation of O VI Velocity with Galactic 
Longitude 

Another question to consider in our quest to under- 
stand how O VI is distributed in the Galaxy is to ask 
whether the velocities of the O VI absorption follow 
that predicted from differential Galactic rotation. 

To construct the LSR velocities v that we would see 
along any particular line of sight, assuming the absorb- 
ing gas co-rotated with the Milky Way disk, we use the 
Galactic rotation curve obtained by Clemens ( 1985} , 
with R() = 8.5 kpc and Q = 220 km s~' . Since v is a 
function of the distance, we can calculate a minimum 
and maximum velocity where we might expect absorp- 



tion. In Figure 17 we plot bars that show the allowed 
LSR velocities as a function of Galactic longitude for 
each target. The variation of permitted velocities with 
distance is not linear — velocities often approach a 
constant value for a given distance, for example, and 
can sometimes backtrack at large distances to the same 
velocities predicted at smaller distances. To show this, 
each bar has a gray-scale in which darker gray col- 
ors mark velocities where the velocity of the gas is 
changing only very slowly with distance, i.e., darkness 
oc \/{dv/dr). [Note, each bar has grayscale coding 
scaled relative to the individual sight line, so a particu- 
lar gray value for one sight line does not represent the 
same \/{dv/dr) along another sight line.] Bars iden- 
tical to these are also shown in the bottom panels of 



Figure 24 although the velocities shown there are he- 



liocentric values. 

On top of each of the predicted-velocity bars, 
we plot the velocity and width of the O VI absorp- 
tion components measured towards a star either as a 
square, if the wavelength scale of the data was cor- 
rected using available STIS data, or a ti'iangl e, if the 
adopted +10 km s 
Appendix |A[ 



shift was applied (see 5 2.4 



and 



If no shift was thought to be appropri- 
ate, we plot a circle. We also plot the velocity of the 
H2 6-0 R(4) line as a cross. Regardless what kind 
of velocity correction was made, the location of the 
H2 velocity with respect to that of the O VI was al- 
ways correct, since the two were adjusted in the same 
manner. All O VI and H2 components are shown, and 
sight lines with no detectable O VI are still included 
(although obviously no O VI positions are marked). 
The figure shows that there is a considerable lack of 
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A +10 km s"^; 
• no shift 

+ H2 R(4) velocity 
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Fig. 17. — : In this figure, permitted velocities for gas rotating in the plane of the Galactic disk are shown as gray 
bars, sorted by Galactic longitude. The measured velocities of O VI absorption lines are overplotted with red triangles, 
squares, or circles (depending on how the offsets for the FUSE wavelength zero-point were calibrated — see legend 
box); R(4) H2 absorption lines are overplotted as blue crosses. (Blue crosses in circles indicate H2 measured in spectra 
which had no shifts applied.) The IDs assigned to the stars (see Table|2]i are given to the left of the Galactic longitude 
of each star, in square parentheses. To show the width of the O VI absorption, a red line extends a distance of ±bi/2 
from the centroid velocity of each absorption component (where bj is the Doppler parameter of component ;)■ 
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Fig. 17. — : continued. 



correlation between the velocities of O VI components 
and those predicted from differential Galactic rotation, 
even if one ignores regions where the gas is violently 
disturbed (i.e., Vela and Carina). The H2 lines often 
agree better with that predicted from rotation, but they 



too can be found outside the permitted range. At some 
longitudes (e.g. I ^ 343°), there appears to be some 
coherence where all the O VI is at the same velocity. 
In most cases, however, this is simply because the stars 
are all at nearly the same longitude. For the 147 indi- 
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vidual O VI components plotted in Figure 17 52 % 
have centroids outside the predicted range. 

Of course, the O VI absorption velocities only show 
the centroid of the absorption (values listed in Table |4] 
but converted to LSR velocities) and its width. As 



can be seen in Figure 24 there is often some overlap 
between velocities allowed by Galactic rotation, and 
some part of an O VI profile. Again, of the 147 in- 
dividual O VI components plotted in Figure 17 only 



16% have Doppler widths which do not extend into the 
region predicted for differential Galactic rotation. If 
the O VI line is comprised of many components, then 
some O VI absorbing regions will be correlated with 
the general kinematics of the ISM. However, Figure 17 
shows that the bulk of the O VI knows little about the 
rotation of the Galactic disk, or that a mechanism ex- 
ists to generate large peculiar motions in clouds that do 
follow Galactic rotation. 

3.12. Comparison of Low and High Ion Velocity 
Extremes 

Figure [Tt] shows that the velocities of the O VI cen- 
troids often fall outside the kinematic limits defined 
by differential Galactic rotation. Also, a direct com- 
parison of H2 and O VI velocities show that the two 
correspond to each poorly. This implies that the bulk 
of the O VI absorption does not mimic the behavior of 
the cool neutral interstellar medium of the Milky Way 
which makes up the Galactic disk. 

Measuring the centroids of the O VI absorption and 
comparing them with either the range of velocities pre- 
dicted from differential Galactic rotation or the veloc- 
ities of other interstellar features tells only part of the 
story, however. If one considers the outer edges of the 
O VI absorption line profiles and compares them to 



these velocities (see the bottom panels of Fig 24 1 the 
differences are even larger 

To see how the extremes of the O VI lines compare 
to the extremes of other lower ionization species, we 
have investigated the absorption line profiles of four 
other species: C III A977, O I A1039, Si III A1206 and 
C II A 1335 lines. All these transitions have large /- 
values and produce extremely strong absorption lines 
with sharp edges. These features are so strongly sat- 
urated that we can expect even relatively low column 
density gas that is well removed from the core veloc- 
ities to yield detectable absorption components. The 
first two of these four lines are recorded in our FUSE 
data, while the second two arise at longer wavelengths 



and were taken from the STIS data for the stars listed 
in Table [TOl Since STIS data were not available for all 
the stars in our FUSE sample, the number of Si III and 
C II lines studied is smaller than the number of C III 
and O I lines. 

The C III A977 line lies in one of the SiC channels, 
and for our analysis we usually used the SiC2A chan- 
nel. To ensure that the wavelength calibration of the 
SiC channel matched that of the LiF channel, we com- 
pared the H2 line nearest in velocity to the C III line 
[the 11-0 P(3) fine at 978.217 A] with the 6-0 P(3) 
fine next to the OVI absorption at 1031.192 A. Match- 
ing these fines in velocity provided the correct relative 
wavelength scales. No relative correction was made 
for the O I A 1039 since the line lies in the same chan- 
nel as the O VI line, and we assumed that the wave- 
length scale was sufficiently accurate at both wave- 
lengths. The absolute wavelength of the C III, O I, and 
O VI lines finally rested on how well the LiF channel 
zero-point was corrected usin g ava ilable STIS data or 
assuming a 10 km s"' shift (^4|. The Si III A1206 
and C II A 1335 lines in the STIS data were assumed to 
have negligible wavelength scale errors. 

To calculate the velocity extremes of all these 
species, vv, we first converted normalized absorption 
line profiles to AOD A^(v) column density profiles. For 

VI, we used the A^(v) profiles for which the H2 lines 
had been removed ( ^2.8.2[ ). We then fitted a high-order 
polynomial to the wings of each A^(v) profile and mea- 
sured the velocity extreme at a given N{v). The values 
of A^(v) used for the O VI, O I, C II, C III and Si III 
lines were: 2 x 10", 5 x 10'^ 1 x 10^^ 4 x 10", and 

1 X 10" cm"^ (km s~')"', respectively. (We comment 
on these values below.) In the few cases where indi- 
vidual absorption components could be seen resolved 
from the main body of the absorption, we used the 
outer edge of that component if it had A^(v) values that 
crossed the adopted A^(v) limits. 

The results are shown in Figure [TS] The error bars 
arise from re-measuring Vg for profiles normalized by 



upper and lower continua error envelopes (see §2.6 1. 



The figure shows that there is a good correlation be- 
tween the velocity extremes of all the species selected 
and those of the O VI absorption. In each panel, we 
show the Pearson correlation coefficients for the posi- 
tive and negative edges considered seperately. The fig- 
ure demonstrates that there is some form of coupling 
between the hot medium and the warm or cool phases. 
This correlation was seen earlier by |Cowie et al.| ( [T979] ) 
from the Copernicus O VI survey, and was thought to 
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Fig. 18.—: Velocity extremes for the profiles of the C II A1335, Si HI A1206, C IH A977 and O I A1039 lines plotted 
against the velocity extremes of the O VI A1032 lines. Circles represent negative edges, while squares represent 
positive edges of the absorption lines. For the positive and negative edges, the Pearson correlation coefficient r is 
shown at the top and bottom of each panel, respectively. The dotted line represents the line of equality for the velocities 
of any two ions. 



arise from conduction/evaporation interfaces between 
the two phases. The interface origin for some of the 
O VI absorption has been confirmed over short paths 
through the local ISM by [Savage & Lehner| ( ,2006; . 

This behavior is not unprecedented. Very strong 
lines from low-ionization species exhibit extraordinary 
dispersions of their velocity extremes, a phenomenon 
studied in detail by Cowie & York ( 1978a b ). One may 
imagine that very small fractions of the low-ionization 
material are stirred up by the same dynamical pro- 
cesses that produce the O VI absorbing material, such 
as the influences from supernova blast waves or mass- 



loss flows from stars, as we have discussed previously. 
In this case, we would expect little correlation between 
the velocity of gas in these structures and the bulk mo- 
tion of the cool ISM, as we discussed above in ^3.11| 



We note, however, that Figure 18 should be inter- 



preted with some degree of caution. The exact value 
of v'e depends on the value of A^(v) at which vv is evalu- 
ated, which in turn can be influenced by the resolution 
of the data. The Si III and C II lines in the STIS data 
have very sharp edges because of the high resolution 
at which they were measured. In these cases, the exact 
value of A^(v) adopted is less important, since the line 
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edges are well defined. For the O I and C III lines mea- 
sured at lower resolution in the FUSE data, the edges 
are smeared by the instrumental LSF. Inspection of the 
profiles used in this study show that these effects are 
Ukely to be small; the values of A^(v') we adopted to 
measure were set so as to ensure that the majority of 
Unes would have reliable measurements of Vg as close 
to Niy) = as possible, given the typical S/N of the 
FUSE and STIS data. 

3.13. Minimum O VI Line Widths 

The fraction of oxygen in the form of O VI 
peaks for coUisionally ionized gas at temperatures of 
logTmax = 5.45 if the gas is in collisional ionization 
equilibrium, at which point 20% of the oxygen is in 
the form of O"""^ ions. This temperature corresponds 
to a Doppler width of b = 0.0321\/r = 16.6 km s"'. 
What is the minimum reliable Doppler parameter in 
our sample? There are several sight lines for which 
we have measured b'^\6 km s"'. For these, how- 
ever, either the O VI line is weak, and the removal 
of strong HD 6-0 R(0) absorption has made the line 
profile uncertain (HD 1 10432), or the line is actually a 
separate component blended with other stronger com- 
ponents (HD 093205, HDE 303308). The sight lines 
with the narrowest, reliable values of b are towards 
HD 046202, with b = 16.4!4-^ km s"' and HD 093250 
with b = 16.7 ± 1 .6 km s"' . These are consistent with 
the smallest values of b found towards WDs (see, e.g. 
Fig. 10 of SL06). There are many more profiles with 
b > 20 km s"', a number which provides a reliable 
lower limit for b. Thus, our smallest value of b is 
consistent with detecting single thermally broadened, 
coUisionally ionized clouds. We cannot be sure how- 
ever, that our narrowest absorption lines actually arise 
from single clouds, since their column densities (0.65 
and 0.83 X 10^"* cm"^) are much greater than the weak 
lines seen towards local WDs and nearby Copernicus 
stars. These sight lines may intercept several clouds 
whose relative motions along the line of sight are at a 
minimum. Alternatively, if sight lines intercept clouds 
of very different sizes ( ^3.9| l the narrowest lines we 
detect may indeed represent absorption from single 
large clouds with moderate column densities. 

3.14. Distribution of O VI Column Density with 
Doppler Parameter, and with Distance 

A Gaussian profile well fits the distribution of 
Doppler parameters b found for the FUSE dataset. 



The profile is centered at (b) = 37 km s"' and has a 
standard deviation of 11 km s"' . However, the Doppler 
parameters are not entirely independent quantaties. A 
correlation between A^(0 VI) and b was first discussed 
by Heckman et al. (2002), who studied the trend using 



extragalactic sight lines and data from early investi- 
gations of O VI in our Galaxy. Although our survey 
does not include O VI absorption systems from clouds 
outside of the Local Group of galaxies, our increased 
sample size of Galactic systems permits a closer ex- 



amination of the correlation. Figure 19 shows that the 
correlation exists in the datasets adopted in this pa- 
per To understand better the origin of the trend, we 
have again indicated which sight lines have ROSAT 
categories R = ot R > 0, and which have distances 
d < 1 kpc or d > 1 kpc. The figure shows that there 
is no dependency on R; it also shows that most low- 
A^(0 VI), low-b points are towards nearby stars. 

We suggest that the N vs. b correlation has a sim- 
ple interpretation: we recall from our findings on the 



lack of change in ctcl with distance (§3.9i that stars 



that are further away intercept not only more O VI 
absorbing structures, but probably also larger, more 
sparsely distributed ones with higher values of O VI 
column densities. On theoretical grounds, it is plausi- 
ble to imagine that these larger clouds have large inter- 
nal motions because they have recently evolved from 
explosive events. The more numerous low-A'(0 VI) 
systems may represent pieces of much older regions 
that have been partitioned by incursions of the am- 
bient, cooler medium, but have not yet had a chance 
to cool appreciably below the temperatures represen- 
tative of O VI. These smaller regions may have had 
their internal velocities moderated because they have 
had a chance to couple dynamically to the surrounding 
quiescent gas phases. Alternatively, we can imagine 
that what we interpret to be large A^(0 VI) systems are 
actually clusters of small clouds, which have cloud-to- 
cloud velocity dispersions that are larger than that of 
the material inside individual clouds. 

4. SUMMARY 

1. We have used FUSE to observe early-type 
(mainly 02-B3) stars in the Galactic disk in or- 
der to characterize O VI A 1032 absorption lines 
arising from interstellar gas in the plane of the 
Milky Way. A total of 120 stars was observed 
as part of FUSE PI team programs PI 02 and 
P122, of which nine were eventually rejected 
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from our survey; the remaining 111 stars con- 
stitute a total investment of 754 ksec of FUSE 
time. We supplemented these sight lines with 
archival data available April 2003 or earlier; our 
final sample consists of 148 stars, most of which 
are at Galactic latitudes \b\ < 10° and distances 
of > 1 kpc. 

2. In order to derive accurate spectroscopic dis- 
tances to the stars, we have used several mod- 
ern databases to improve upon published dis- 
tances (Appendix [B). Hipparcos and Tycho cat- 
alogs have provided precise optical magnitudes, 
as well as identification of multiple systems and 
variable stars. 2MASS infrared magnitudes have 
enabled us to estimate the interstellar V-band 
extinction along a sight line, and the absolute 
magnitudes of stars — given a spectral type and 
luminosity class — have been constructed by 



collating available surveys from the literature. 
Despite these improvements, imprecise knowl- 
edge of stellar types and classes continues to 
have the greatest uncertainty in deriving a spec- 
troscopic distance, limiting the precision in a 
distance measurement to 10-30%. 



Where appropriate, we have added the O VI 
measurements from several other O VI stud- 
ies to our sample, in order to probe distances 
not covered by our FUSE stars. These include 
nearby (~ 20-100 pc) white dwarf stars (WDs) 
( Savage & Lehner|2006 1, intermediate distance 
(~ 100- 1000 pc) stars observed by the Coper- 
nicus satellite (Jenkins 1978a), distant halo stars 



( Zsargo et al.|2003), and extragalactic sight lines 



Savage et al.p003| [Wakker et al.|[2003] l. For 



many of the stars (but not including the WDs), 
we have rederived distances using the same set 
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of tools used to derive distances towards the 
FUSE sample of stars. 

ROSAT X-ray maps show that many stars lie to- 
wards regions of enhanced X-ray emission. This 
emission likely arises from hot bubbles blown 
out from the star or its stellar association. Other 
sight lines, however, are found at positions with 
very little X-ray emission, which we label as 
R = sight lines. In contrast, lines of sight in- 
tercepting X-ray emitting regions have designa- 



tions ofR > (5 3.4 and § 3.5 1. The average vol 



ume density towards stars in the Galactic disk, 
n = N(0 VYf/d, depends on whether we use 
R = Qor R > Q sight lines, and whether we sub- 
tract a contribution to the O VI column density 
from the Local Bubble (LB; ^TT) . When we 
subtract a value of A^lb(0 VI)= L 1 1 x 10'^ cm"^ 
arising from the LB, we find little difference in n 
between stars with distances 0.2 < d < 2.0 kpc 



and those with d > 2 kpc (§3.7.2 1; we also find 



that for stars d > 2 kpc, the R > sight lines 
have slightly higher densities than R = sight 
lines, even after excluding stars in the direction 
of the Vela SNR and Carina Nebula, which have 
the highest densities in our sample ( ^3.7.3| l. This 
is probably because the circumstellar environ- 
ment of the star contributes a small excess to the 
O VI column density. 

Our data well fit the relationship A^(0 VI)sin \b\ = 
no/i(l - exp"!'!'"'), which naturally arises from 
assuming that the density n of O VI absorb- 
ing gas falls with height |z| above the Galactic 
plane as n = noexp"!'''/''. Here, no is the mid- 
plane density and h is the characteristic scale 
height. We find no = 1.3 x 10"^ cm"^ [after 
subtracting a contribution to A^(0 VI) from the 
LB and using the R = sight lines; 5 3.8.3 1, but 



that h depends on whether we use extragalactic 
lines of sight in the northern (N) or southern (S) 



Galactic hemispheres (as first noted by [Savage 
let al.|2003| —see pXTj i. We derive values of 
h = [4.6,3.2] kpc when using the [N,S] sample. 
We also find it neccessary to include a dumpi- 
ness factor acL to represent the fact that the dis- 
persion in log[A^(0 VI)sin|fo|] around the theo- 
retical fit is much larger than can be accounted 
for by experimental errors alone ( ^3.8.2[ ). Our 
best estimate is that (Tcl — 0.26 ± 0.02 dex. 

6. A^(0 VI) and the distance to a star are correlated 



(i 3.9 1. This shows that the processes which give 
rise to O VI absorption are ubiquitous over the 
entire Galactic disk, and certainly do not arise 
primarily from circumstellar environments. The 
only inhomogeneity of the OV I volume density 
n can be found in a concentration of A^(0 VI) 
upper limits (10 out of 23 sight lines) towards 
d > 1 kpc FUSE sight lines in the region / = 
80° - 140°, the outer Perseus arm of the Galaxy. 
There is no evidence that the volume density of 
O VI is influenced by the presence of Galactic 
spiral arms along the stellar sight lines. 

7. Although O VI is distributed smoothly enough 
for N{0 VI) to correlate with distance, the rela- 
tive dispersion on either side of this correlation 
does not decrease with distance, as would be 
expected if more distant sight lines intercepted 
more clouds. Quantitively, we find no differ- 
ence in the dispersion (ctcl) between nearby 
(0.2-2 kpc) and distant (> 2 kpc) stars. This 
indicates that in addition to a random placement 
of small clouds, there may also be sparsely dis- 
tributed large clouds. Or, alternatively, some of 
the small clouds may simply be clustered over 
very large scales. 

8. The velocity edges of the O VI lines and those 
of the strong transitions of C III A977.02, 
O I A1039.23, Si III A1206 and C II A1335, 
are all correlated. This suggests that, at least to 
first order, the processes that move the hot gas to 
high velocity also affect the velocities of other 
species ( ^3.12[ ). 

9. The smallest reliable O VI Doppler parameter 
we find in our sample is b ~ 16-17 km s~', 
which is consistent with that expected for a ther- 
mally broadened cloud at a temperature close 
to the peak in the ionization fraction curves, 
logllnax — 5.45 (5 3.131. We also confirm ear- 



lier research which found that b and A^(0 VI) are 
correlated. We show, however, that the width of 
O VI lines cannot be accounted for by differen- 
tial Galactic rotation; the lines are far wider than 
expected if absorption came from a smoothly 
distributed ISM corotating with the disk of the 
Milky Way (§3.14i. The correlation suggests 



that the large clouds mentioned in Point|7]above 
have large internal motions, or, alternatively, 
that the small clouds within clusters are mov- 
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ing more rapidly than the gases within individ- 
ual clouds. 
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Table 2 

Sight Lines Observed with FUSE Sorted By Name, and Observation Log 



Star 


FUSE 


Unique 


I 


b 


Start 




Name 


Prog ID 


ID'' 


(deg) 


(deg) 


Date 


(ksec)*" 


BD+532820 


P12232 


34 


101.24 


-1.69 


2000 08 06 


11.3 


CPD-417712 


B02506 


134 


343.41 


1.17 


2001 08 09 


9.6 


CPD-592600 


P12214 


90 


287.60 


-0.74 


2000 03 21 


15.6 












2000 03 23 














2000 03 25 




CPD-592603 


P12215 


88 


287.59 


-0.69 


2000 03 21 


6.0 


HD000108 


Z90101 


43 


117.93 


1.25 


2002 08 31 


7.9 




D06401 








2003 08 05 




HD005005A 


P10201 


44 


123.13 


-6.24 


2000 08 14 


6.4 


HD012323 


P10202 


46 


132.91 


-5.87 


1999 11 25 


3.9 


HDO 13268 


PI 0203 


47 


133.96 


-4.99 


1999 11 24 


4.4 


HDO 13745 


P10204 


48 


134.58 


-4.96 


1999 11 24 


4.4 


HD014434 


P10205 


49 


135.08 


-3.82 


1999 1124 


4.4 


HD015137 


PI 0206 


51 


137.46 


-7.58 


1999 11 25 


2.3 


HD017520 


B09702 


50 


137.22 


0.88 


2001 09 11 


6.6 


HD030677 


PI 0208 


56 


190.18 


-22.22 


2001 01 01 


4.2 


HD034656 


P10113 


53 


170.04 


0.27 


2000 03 01 


4.2 


HD039680 


P10209 


57 


194.08 


-5.88 


2001 03 05 


3.9 


HD041161 


P10210 


52 


164.97 


12.89 


2003 09 25 


6.0 


HD042087 


P21603 


54 


187.75 


1.77 


2001 10 15 


2.9 


HD042088 


P10211 


55 


190.04 


0.48 


2000 1 1 05 


4.2 


HD042401 


P10212 


59 


197.64 


-3.33 


2000 1105 


1.4 


HD045314 


P10213 


58 


196.96 


1.52 


2000 03 16 


5.5 


HD046150 


P10214 


62 


206.31 


-2.07 


2001 03 05 


8.6 




CI 6802 








2004 02 25 




HD046202 


P11610 


63 


206.31 


-2.00 


2001 10 16 


4.9 


HD047088 


P15104 


61 


205.92 


-0.40 


2001 02 20 


4.6 


HD047360 


P10215 


64 


207.33 


-0.79 


2001 03 04 


5.0 


HD047417 


P10216 


60 


205.35 


0.35 


2000 03 15 


5.2 


HD052463 


P10218 


70 


238.90 


-10.52 


2000 01 24 


3.1 


HD058510 


P10219 


68 


235.52 


-2.47 


2001 04 06 


5.7 


HD060196 


P12209 


73 


242.95 


-4.65 


2000 04 16 


5.7 


HD060369 


P10502 


72 


242.68 


-4.30 


2000 04 17 


6.0 


HD061347 


PI 0220 


67 


230.60 


3.79 


2000 04 15 


6.8 


HD062866 


P12210 


69 


237.48 


1.80 


2000 04 02 


8.1 












2000 0411 




HD063005 


P10221 


71 


242.47 


-0.93 


2000 04 05 


5.3 


HD064568 


P12211 


74 


243.14 


0.71 


2000 04 03 


16.0 












2000 04 05 














2000 04 11 




HD065079 


P10222 


65 


217.87 


15.93 


2000 03 07 


4.4 


HD066788 


P10118 


75 


245.43 


2.05 


2000 04 06 


4.2 


HD069106 


P10223 


76 


254.52 


-1.33 


2002 12 31 


0.4 


HD074711 


P10225 


78 


265.74 


-2.61 


2000 01 26 


4.7 
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Table 2 — Continued 



Star 


FUSE 


Unique 


/ 


b 


Start 




Name 


Prog ID 


ID" 


(deg) 


(deg) 


Date 


(ksec)* 


HD074920 


P10226 


77 


265.29 


-1.95 


2000 01 26 


4.7 


HD075309 


PI 0227 


79 


265.86 


-1.90 


2000 01 26 


4.7 


HD088115 


P10123 


82 


285.32 


-5.53 


2000 04 04 


4.5 


HD089137 


P10228 


80 


279.70 


4.45 


2000 01 26 


4.6 


HD090087 


PI 0229 


81 


285.16 


-2.13 


2000 04 03 


3.1 


HD091597 


P10230 


84 


286.86 


-2.37 


2000 02 04 


14.3 












2000 02 06 














2000 02 08 




HD091651 


P10231 


83 


286.55 


-1.72 


2000 05 25 


18.1 












2000 05 27 




HD092554 


PI 0232 


89 


287.60 


-2.02 


2000 05 26 


14.7 


HD093129A 


PI 1702 


85 


287.41 


-0.57 


2000 01 27 


3.8 


HD093146 


P10233 


92 


287.67 


-1.05 


2000 02 01 


4.1 


HD093205 


P10236 


87 


287.57 


-0.71 


2000 02 01 


4.7 


HD093206 


PI 0234 


93 


287.67 


-0.94 


2000 02 01 


4.1 


HD093222 


PI 0237 


94 


287.74 


-1.02 


2000 02 03 


2.9 


HD093250 


P10238 


86 


287.51 


-0.54 


2000 02 04 


4.1 


HD093827 


P10239 


96 


288.56 


-1.54 


2000 02 02 


4.2 


HD093843 


P10240 


95 


288.24 


-0.90 


2000 02 04 


4.1 


HD094493 


P10241 


97 


289.02 


-1.18 


2000 03 26 


4.4 


HD096670 


P10242 


99 


290.20 


0.40 


2000 02 01 


4.3 


HD096715 


P10243 


100 


290.27 


0.33 


2000 02 01 


4.6 


HD096917 


P10244 


98 


289.29 


3.06 


2000 05 26 


8.0 


HD097913 


P12217 


101 


290.84 


1.41 


2000 05 26 


16.0 


HD099857 


PI 0245 


105 


294.78 


-4.94 


2000 02 05 


4.3 


HD099890 


P10246 


102 


291.75 


4.43 


2000 03 25 


4.6 


HD100199 


P12218 


104 


293.95 


-1.48 


2000 03 24 


4.6 


HD100213 


PI 0247 


109 


294.81 


-4.14 


2000 03 24 


4.6 


HD 100276 


PI 0248 


103 


293.31 


0.77 


2000 02 01 


4.1 


HD101131 


PI 0249 


106 


294.78 


-1.62 


2000 03 23 


4.3 


HD101190 


P10250 


107 


294.78 


-1.49 


2000 02 08 


3.4 


HD101205 


P10251 


110 


294.85 


-1.65 


2000 05 27 


0.1 


HD101298 


PI 0252 


111 


294.94 


-1.69 


2000 03 23 


4.6 


HD101413 


PI 0253 


112 


295.03 


-1.71 


2000 03 24 


4.2 


HD101436 


PI 0254 


113 


295.04 


-1.71 


2000 02 05 


4.2 


HD102552 


P10255 


114 


295.21 


1.35 


2000 03 19 


4.1 


HD 103779 


P10256 


115 


296.85 


-1.02 


2000 04 09 


47.5 




D07301 








2004 04 09 




HD 104705 


PI 0257 


116 


297.46 


-0.34 


2000 02 05 


4.5 


HD108639 


A12013 


117 


300.22 


1.95 


2000 03 20 


5.0 


HDl 10432 


P11614 


118 


301.96 


-0.20 


2000 04 04 


3.6 


HDl 12784 


Z90150 


119 


304.00 


2.26 


2002 04 21 


1.4 


HDl 14441 


P10258 


121 


305.80 


7.40 


2000 03 23 


4.4 


HDl 15071 


P10259 


120 


305.77 


0.15 


2000 05 14 


5.9 
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Table 2 — Continued 



Star 
N3111C 


FUSE 
Prog ID 


T Tniniip 
ID^ 


; 


1, 


Start 
Date 


T 


HD 115455 


A12007 


122 


306.06 


0.22 


2000 05 14 


5.7 


HD 116538 


PI 0260 


124 


308.23 


10.68 


2000 05 24 


5.9 


HDl 16781 


P10261 


123 


307.05 


-0.06 


2000 04 06 


3.0 


HDl 18571 


P12220 


126 


308.70 


1.35 


2000 03 25 


4.4 


HDl 18969 


P1222I 


125 


308.52 


-1.39 


2000 03 25 


11.2 


HDl 22879 


B07105 


127 


312.26 


1.79 


2002 03 03 


1.5 


HD124314 


PI 0262 


128 


312.67 


-0.42 


2000 03 22 


4.4 


HD124979 


PI 0263 


129 


316.41 


9.08 


2000 07 01 


5.1 


HD134411 


P10265 


131 


330.07 


15.53 


2000 07 10 


3.9 


HD148422 


P10150 


130 


329.92 


-5.60 


2000 04 08 


8.3 












2000 04 10 














2000 04 12 




HD151805 


PI 0266 


132 


343.20 


1.59 


2000 07 14 


2.9 


HD151932 


PI 1708 


133 


343.23 


1.43 


2000 07 12 


2.3 


HDl 52200 


B02501 


135 


343.42 


1.22 


2001 08 09 


5.4 


HD152218 


P10154 


140 


343.53 


1.28 


2000 09 08 


9.5 


HD152233 


PI 0267 


138 


343.48 


1.22 


2000 09 08 


4.1 


HD152245 


Z90160 


141 


344.45 


2.02 


2002 08 14 


1.6 


HDl 52248 


P10268 


136 


343.47 


1.18 


2000 07 12 


4.4 


HD152314 


PI 0269 


139 


343.52 


1.14 


2000 07 12 


4.0 


HD152623 


PI 0270 


142 


344.62 


1.61 


2000 07 11 


6.1 


HD152723 


P10271 


143 


344.81 


1.61 


2000 09 09 


5.0 


HD153426 


P10272 


145 


347.14 


2.38 


2000 03 31 


5.7 












2000 07 14 




HD156292 


P10274 


144 


345.35 


-3.08 


2000 04 03 


10.8 












2000 04 05 




HD157857 


P10275 


11 


12.97 


13.31 


2000 09 02 


4.0 


HDl 58661 


PI 2222 


7 


8.29 


9.05 


2000 09 08 


6.6 


HD161807 


PI 2223 


146 


351.78 


-5.85 


2000 08 17 


2.0 












2000 09 08 




HD163758 


P10159 


147 


355.36 


-6.10 


2000 08 16 


5.4 


HD163892 


P10276 


6 


7.15 


0.62 


2001 04 27 


12.0 












2001 04 29 




HDl 64906 


PI 0277 


3 


6.05 


-1.33 


2000 08 31 


5.3 


HD165052 


P10278 


4 


6.12 


-1.48 


2000 08 31 


4.9 


HDl 65955 


PI 0279 


148 


357.41 


-7.43 


2000 08 17 


4.1 


HD167287 


P12226 


8 


11.99 


-0.93 


2000 08 31 


5.8 


HDl 67402 


P10162 


1 


2.27 


-6.39 


2000 08 31 


3.9 


HDl 67659 


P10280 


9 


12.20 


-1.27 


2000 05 18 


5.8 


HDl 67771 


P10281 


10 


12.70 


-1.13 


2000 08 31 


3.8 


HD168080 


P12227 


12 


13.11 


-1.27 


2000 08 31 


5.7 


HD168941 


P10165 


2 


5.82 


-6.31 


2000 08 30 


8.1 












2000 09 01 




HD175754 


P10168 


13 


16.39 


-9.91 


2001 04 29 


1.9 
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Table 2 — Continued 



Star 


FUSE 


T Tninnf* 


/ 


I, 


Start 


T 

-'exp 


Name 


Prog ID 


ID^ 






Date 




HD178487 


P10172 


14 


25.79 


-8.56 


2000 08 27 


8.9 


HD185418 


P11623 


15 


53.60 


-2.17 


2000 08 10 


4.5 


HD 186994 


P21608 


27 


78.62 


10.06 


2001 07 02 


0.5 












2001 09 07 




HDl 87282 


P11712 


16 


55.62 


-3.79 


2000 08 1 1 


5.2 


HD 187459 


P10282 


19 


68.81 


3.85 


2000 08 10 


9.9 


HDl 90429 


P10284 


21 


72.59 


2.61 


2000 07 18 


5.4 


HD190918 


P10285 


22 


72.65 


2.06 


2000 07 18 


5.5 


HD191495 


P12229 


23 


72.74 


1.41 


2000 08 10 


12.2 


HD191765 


P11713 


24 


73.45 


1.55 


2000 07 19 


4.5 


HD191877 


P10287 


17 


61.57 


-6.45 


2000 06 05 


6.1 


HD192035 


P10286 


28 


83.33 


7.76 


2000 06 17 


17.6 












2000 06 19 














2000 06 21 




HD192639 


PI 1624 


25 


74.90 


1.48 


2000 06 12 


19.4 




C17101 








2002 09 04 




HD195965 


P10288 


30 


85.71 


5.00 


2000 06 20 


36.5 












2000 11 08 




HD199579 


P11625 


29 


85.70 


-0.30 


2000 07 19 


4.3 


HD201345 


P12230 


26 


78.44 


-9.54 


2000 06 13 


5.1 


HD202347 


P10289 


31 


88.22 


-2.08 


2000 06 20 


4.1 












2000 09 13 




HD207198 


P11628 


37 


103.14 


6.99 


2000 07 23 


10.2 


HD210809 


P12231 


32 


99.85 


-3.13 


2000 08 05 


12.1 












2000 08 07 




HD210839 


P11631 


38 


103.83 


2.61 


2000 07 22 


6.1 


HD212044 


P12234 


33 


100.64 


-4.35 


2000 07 23 


4.5 


HD216044 


P12238 


39 


105.94 


-3.64 


2000 07 23 


5.6 


HD218915 


P10188 


40 


108.06 


-6.89 


2000 07 23 


5.4 


HD224151 


P12241 


41 


115.44 


-4.64 


2000 08 1 1 


18.2 












2000 08 13 




HD224868 


P12202 


42 


116.87 


-1.44 


1999 11 28 


5.1 


HD326329 


B02505 


137 


343.47 


1.17 


2001 08 09 


5.5 


HDE225757 


P10177 


20 


69.64 


4.85 


2000 08 06 


12.1 












2000 08 09 




HDE232522 


P12201 


45 


130.70 


-6.71 


1999 11 28 


2.1 


HDE235783 


P12233 


35 


101.69 


-1.87 


2000 08 05 


11.3 


HDE235874 


P12237 


36 


101.97 


-5.93 


2000 07 23 


5.4 


HDE303308 


P12216 


91 


287.60 


-0.61 


2000 05 25 


13.0 












2000 05 27 




HDE308813 


P12219 


108 


294.80 


-1.61 


2000 03 23 


14.3 


HDE315021 


P12224 


5 


6.12 


-1.33 


2000 08 30 


5.4 


HDE332407 


P12228 


18 


64.28 


3.11 


2000 06 10 


8.5 












2000 06 12 
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Table 2 — Continued 



Star 


FUSE 


Unique 


I 


b 


Start 


Texp 


Name 


Prog ID 


ID" 


(deg) 


(deg) 


Date 


(ksec)'' 


LS277 


PI 2208 


66 


223.34 


1.57 


2000 03 31 


3.3 



"ID number corresponding to panel numbers shown in Fig. [24] as well as 
entries in Tables[T]|4]|5] and[8] 

''Total exposure time from coadded spectra actually used, summed over all 
dates and FUSE programs 

Note. — All observations were made with the LWRS aperture except for 
HD175754, HD04I161, and HD069106, which were observed using the 
MDRS aperture. 
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Table 3 

uv lines used in the analysis of fuse data 





Ao 




Line 


(A) 


/-value 


HD 7-0 R(0) 


1021.453 


0.024 


H2 7-0 R(6) 


1030.080 


not used 


H2 6-0 P(3) 


1031.192 


0.011 


ClI A1031.5 


1031.507 


0.151 


H2 8-0 R(8) 


1031.557 


0.019 


HD 6-0 R(0) 


1031.909 


0.023 


VIA1031 


1031.927 


0.133 


Hz 6-0 R(4) 


1032.354 


0.017 


VI A1037 


1037.617 


0.066 


H2 7-0 R(8) 


1042.745 


0.018 


H2 5-0 R(4) 


1044.542 


0.015 


HD 4-0 R(0) 


1054.286 


0.016 


H2 6-0 R(8) 


1054.520 


0.017 


H2 4-0 P(3) 


1056.472 


0.009 


H2 4-0 R(4) 


1057.376 


0.013 


H2 4-0 R(8) 


1079.932 


0.012 


CI I A1347'' 


1347.240 


not used 



''Measured in STIS data, used for cali- 
bration of zero point in FUSE data. 
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Table 5 

AOD Column Density Upper Limits Towards Stars With No Detected O VI 



ID 


Star 


Ni 


N 


M, 


a(N) 


a(N)T 


A'lirait'' 










(10" 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


15 


HD185418 


-3.11 


0.58 


3.95 


0.45 


3.56 


8.52 


29 


HD199579 


-0.17 


0.55 


1.37 


0.28 


0.82 


2.29 


31 


HD202347 


-0.94 


0.22 


1.29 


0.16 


1.13 


2.72 


37 


HD207198 


-2.01 


0.17 


2.28 


0.94 


2.34 


5.45 


38 


HD2 10839 


-1.50 


-0.67 


0.12 


0.34 


0.88 


1.59 


43 


HD000108 


-2.24 


-0.33 


1.43 


0.69 


1.96 


4.25 


44 


HD005005A 


-0.74 


0.23 


1.11 


0.44 


1.02 


2.49 


46 


HDO 12323 


-2.12 


-0.49 


1.10 


0.44 


1.67 


3.47 


47 


HDO 13268 


-1.18 


0.25 


1.62 


1.10 


1.78 


4.22 


48 


HDO 13745 


-4.00 


1.27 


5.94 


1.10 


5.09 


12.52 


49 


HDO 14434 


-9.44 


-2.17 


3.95 


1.72 


6.91 


14.30 


51 


HD015137 


-4.94 


-1.06 


2.26 


0.99 


3.74 


7.81 


52 


HD041161 


-1.78 


-0.52 


0.69 


0.27 


1.26 


2.54 


59 


HD042401 


-2.64 


0.32 


2.22 


0.70 


2.53 


5.98 


70 


HD052463 


-0.09 


1.06 


2.20 


0.31 


1.18 


3.53 


114 


HD102552 


-0.00 


2.85 


5.29 


1.12 


2.87 


8.81 


121 


HDl 14441 


-0.44 


1.38 


2.88 


0.41 


1.71 


4.96 



^2(7 confidence upper limit. See ^2.8.2| and Appendix|D]for details. 
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Table 6 

Rejected Stars from the FUSE Disk Survey 



Prog 


Star 


I 


b 




Start 




ID 


Name 


(deg) 


(deg) 


Sp/L" 


Date 


(ksec) 


P10207 


HD015642 


137.09 


-4.73 


09.5III:n 


2000 09 15 


5.1 


P10224 


HD074194 


264.04 


-1.95 


08.5Ibf 


1999 12 27 


3.4 


P10235 


HD093204 


287.57 


-0.71 


05Vf 


2000 02 01 


4.7 


PI 2205 


HD035215 


176.46 


-3.19 


BIV 


2000 02 29 


10.1 


P12212 


HD066695 


245.01 


2.21 


B0.5IV 


2000 01 24 


8.2 


P12225 


HD 166546 


10.36 


-0.92 


09.5II-III 


2000 09 03 


3.4 


P12235 


BD+523210 


102.44 


-3.37 


BIV 


2000 08 06 


15.7 












2000 08 09 




P12236 


BD+532885 


102.75 


-2.93 


B2III 


2000 08 02 


11.8 












2000 08 05 




P12239 


HD216438 


105.72 


-5.12 


Bin-in 


2000 07 22 


5.1 



^Spectral Type and Luminosity Class of star 
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Table 7 

Stars observed at more than one epoch 



Unique 










FUSE 








N(0 VI)" 


ID 


Star 


Date 




root 


Aperture 


Channel 


(ksec) 


(10'* cm-2) 


84 


HD091597 


2001-02- 


-04/06/08 


PI 023002/3/4 


LWRS* 


LiFlA 


14.3 


3.17±0.12 






2003- 


06 


-13 


D0640301 


LWRS 


LiFlA 


13.5 


3.14±0.13 


76 


HD069106 


2001- 


■04 


-04 


P1022301 


LWRS 


LiFlA 


0.06 


1.99 ±0.28 






2002- 


■12 


-31 


P1022302 


MDRS'' 


LiFlA 


0.4 


2.12±0.23" 


62 


HD046150 


2001- 


■03 


-05 


P1021401 


LWRS* 


LiFlA 


4.9 


0.96±0.16 






2004- 


■02 


-25 


CI 680201 


LWRS* 


LiFlA 


3.7 


0.80±0.13 






2004- 


■02 


-24 


PI 02 1402/3 


MDRS 


LiFlA 


16.3 


0.90 ±0.22' 


25 


HD 192639 


2000- 


06 


-12 


PI 162401 


LWRS* 


LiFlA 


4.8 


0.90±0.39 






2002- 


09 


-04 


C1710101 


LWRS* 


LiFlA 


14.6 


0.25 ±0.42 


43 


HD000108 


2002- 


08 


-31 


Z9010101 


LWRS* 


LiFlA 


4.8 


<0.63 






2003- 


08 


-05 


D0640101 


LWRS* 


LiFlA 


3.1 


<0.62 


115 


HD103779 


2000- 


■04 


-09 


P1025601 


LWRS* 


LiFlA 


4.3 


0.54 ±0.11 






2004- 


04 


-09 


D0730101/2 


LWRS* 


LiFlA 


43.2 


0.46 ±0.13 



"Column densities derived from profile fitting, with HD 6-0 R(0) line removed, as described in ^2.8. 1| 
''These datasets were eventually coadded for each star to produce a final spectrum. 
■^Assuming a MDRS instrumental resolution of 3 km s ' FWHM jFeldmanetal.|200l| . 
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Table 8 

FUSE Sight Lines Projected Through SNRs" 



Unique 


Target 


Star's distance 


SNR 


SNR's distance 


Ref.*" 


ID 


star 


(kpc) 


Went. 


(kpc) 




12 


HD168080 


2.7 


G13.3-1.3 


2.0-4.7 


1 


15 


HD185418 


0.8 


G53.6+5.7' 


6.7 


2 


29 


HD199579 


1.1 


G85.9 -0.6 


5. 


3 


60 


HD047417 


1.2 


G205.5 +0.5^ 


0.8 


4,5 


11 


HD074920^ 


1.5 


G266.2 -1.2 


0.3'' 


6 


79 


HD075309 


1.5 


G266.2-1.2 


0.3'' 


6 


128 


HD124314 


1.1 


G3 12.4 -0.4 


7 




131 


HD134411 


2.1 


G330.0 +15.08 


0.5: 


7 



^Excluding sight lines intercepting Vela 

''References for distance estimates: [1] Seward et al. 1995, [2] Rosado 1983, [3] 
Kothes et al. 2001, [4] Kirshner, Gull, & Parker 1978, [5] Lozinskaya 1981, [6] Red- 
man et al. 2002, [7] Milne 1979 

'^3C400.2 

""Monoceros Nebula (not to be confused with the much larger Monogem Ring) 

^The star is located on the edge of the SNR, so the sight-line interception, even if the 
SNR is in the foreground, is uncertain. 

^This SNR is presumed to be within the Vela SNR; see text for details. 

^The Lupus Loop 
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Table 9 

Line of sight volume densities of O VI in the Galactic Plane 







<log[N{OWT)/d]> loj 


I Y.N{0 viysc/ 




No. 


[log (cm-^)] 


[lo 


g (cm 3)] 


(1) 


(2) 


(3) 




(4) 


stellar distances > 2000 pc 


+LB, All R 


102 


-7.87 ±0.32 




-7.79 


+LB, R = Q 


56 


-7.92 ±0.31 




-7.79 


+LB, R>Q 


46 


-7.82 ±0.32 




-7.79 


-LB, All R 


96 


-7.92 ±0.36 




-7.81 


-LB, R = Q 


55 


-7.97 ±0.36 




-7.81 


-LB, R>Q 


41 


-7.87 ±0.36 




-7.62 


200 pc < stellar distances < 2000 pc 


+LB, All R 


88 


-7.75 ±0.31 




-7.69 


+LB, R = 


42 


-7.78 ±0.33 




-7.75 


+LB, R>0 


46 


-7.73 ±0.29 




-7.64 


-LB, All R 


81 


-7.91 ±0.45 




-7.72 


-LB, R = 


37 


-7.96 ±0.48 




-7.81 


-LB, R>0 


44 


-7.86 ±0.42 




-7.66 



Note. — Description of columns: (1) — "-LB" indicates that 
a contribution from the Local Bubble, [A^lb(0 VI)= 1.11 x 
10'^^ cm"^] has been removed from all column densities, while 
"+LB" indicates that no such correction has been made; R indi- 
cates the ROSAT class assigned in 5 3.4 and 5 3.5 for the -LB sam- 
ples, distances have also been reduced by the radius of the LB, 
which we take to be 100 pc; (2) — number of points used in sta- 
tistical test. (Note that the numbers in the H-LB sample for a given 
ROSAT class in each distance catagory are not quite the same as 
the number for the equivalent -LB class; this is because subtract- 
ing 100 pc for the -LB stars moves a few objects either from the far 
to near sample, or out of the near sample altogether); (3) — mean 
of logA^(0 y\)/d and the standard deviation of the distribution of 
column densities (see 5 3.7.1 1; (4) — distance weighted average of 
O VI column densities. Note: None of the statistics above include 
sightlines to Carina or Vela. 



A. THE ZERO-POINT OF THE FUSE WAVELENGTH SCALE 



As discussed in § 2.4 measuring precise velocities for O VI absorption lines potentially requires a correction for any 
errors which might exist in the FUSE wavelength-scale zero-point associated with CalFUSE v2.0.5. Fortunately, many 
of the stars in our sample have also been observed with the Space Telescope Imaging Spectrograph (STIS) aboard HST, 
which has a well determined wavelength scale. In particular, a substantial number of E140M and E140H echelle data 
have been recorded which cover the CI I A 1347 line. According to the STIS Handbook, the relative accuracy of the 
STIS echelle modes is 0.25-0.5 pixels across the full format of the MAMA detectors, and the absolute accuracy is 
0.5 - 1 .0 pixels ([Kim Quijan o|2004|l. Somew hat larger wavelength-scale errors have been found in some STIS echelle 
spectra (^ Jenkins & Tripp||2001 Tripp et aLl|2005 1 indicating relative wavelength scale accuracies « 1 pixel. To be 
conservative, we assume that the uncertainty in the STIS wavelength calibration is pixel. This corresponds to 1.5 
and 3.0 km s"' for the E140H and E140M STIS gratings, respectively. 

To derive corrections to the FUSE wavelength scale, we compared absorption lines of neutral chlorine with the 
molecular hydrogen lines flanking the O VI features, along individual stellar sight lines; H2 clouds are efficient pro- 
ducers of CI I (Jura 1974; Jenkins 1995), so a direct comparison between the H2 lines which flank the O VI A1032 
absorption in the FUSE data, and a strong (but unsaturated) CI I line whose velocity is known precisely, should provide 
a suitable anchor to the FUSE wavelength scale. While in principle we could have used the CI I A1031.5 A line (see 
Table|3]and Figure [TJ, it was often too weak, or absent altogether, for its wavelength to be measured. 

The STIS data are recorded at higher resolution than the FUSE spectra, so in order to compare the H2 and CI I lines 
directly, we smoothed the STIS data, scaled the CI I to match the strength of the H2 6 - R(4) and P(3) lines, and 
shifted the FUSE spectrum in wavelength space until the lines overlapped. The shape of the smoothed CI I line well 
matched the FUSE H2 line profiles in all cases, and we were able to record shifts to within ±0.2 pixels, or just over 
1 km s~' at the wavelength of the O VI A 1032 line. 

Apart from 46 sight lines in our Galactic disk sample which had also been observed with STIS, we were able to 
examine FUSE and STIS spectra of a few other non-survey stars. A total of 55 FUSE and STIS spectra enabled us 
to investigate whether the errors in the FUSE wavelength scale were dependent on the state of the spectrograph. We 
found a crude correlation between the more extreme shifts in the wavelength zero-point and the X-position of the Focal 
Plane Assembly, as recorded in the FPALXPOS keyword of the LiFl A headers. The data used in the study are listed in 
Table 10 while the measured shifts and their values of FPALXPOS are shown in Figure 20 We found no correlation 
between measured shifts and the Z-position of the FPA, as recorded in the data header keyword FPALZPOS. 

We note the following in Figure [20| The error bars in the ordinate direction are not statistical errors. As explained in 
|2.3| we used the sub-spectrum with the highest S/N ratio to act as a template against which the other sub-exposures 



were compared before sub-exposures were coadded. The bars in Figure 20 represent the minimum and maximum 
excursions in the values of the shifts measured for a set of sub-exposures. 

The majority of the FUSE spectra were taken with FPALXPOS in the range 117.0-175.0. For these data, a con- 
sistent shift of ^ +10 km s"' is required to correct the FUSE wavelength scale. (Again, this shift was not necessarily 
the shift of the first sub-exposure taken of a star, but was derived from the sub-exposure with the highest S/N.) The 
histogram on the right-hand side of Figure 20 shows the distribution of the shifts, along with a Gaussian fit to the 
distribution, which peaks at 9.2 km s"' and has a width of cr = 5.1 km s"'. This implies that a +10 km s"' shift to FUSE 
data which have no companion STIS data will provide an accurate re-calibration of the wavelength scale to within 
±5kms-'. 

With this information, we applied the following re-calibration to the survey stars. For the 46 stars also observed by 
STIS, the FUSE data were shifted by the numbers given in Table 10 Nine stars were observed either with the MDRS 
aperture, with a channel other than the LiFlA, or when the X-position of the FPA was > 200; no shifts were applied 
to these data. The remaining 93 spectra were shifted in wavelength by 0.0344 A, or 10 km s"' at 1032 A. 

Finally, we note that later versions of CalFUSE (v3.0) automatically applied a correction of 10 km s"^ to the zero- 
point of the wavelength calibration to correct for the offset described above. Versions of CalFUSE later than v3.2 used 
an improved scheme for the wavelength calibration of the raw data as a whole, reductions which are described in detail 
by |Dixon et aL] ( |2007l l. 
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No. 




Fig. 20. — : Plot of the shifts required to accurately calibrate the zero-point of the FUSE wavelength scale around 
the O VI A1032 line, as a function of X-position of the Focal Plane Assembly (as determined from the FPALXPOS 
keyword in the LiFlA headers). Values of FPALXPOS for our data fall close to one of four values: 117, 135, 175 
and 227. To better show the error bars in the Y direction (see text for a full explanation of the error bars), we have 
randomized the FPALXPOS values by small amounts to make the points slightly offset from each other. It should 
be understood that points within a particular gray box actually have FPALXPOS values almost identical to the value 
shown on the X-axis at the bottom of each gray box. The distribution of the shifts is shown as a histogram to the right. 
For spectra with FPALXPOS in the nominal range < 175, the shift is centered at 9.2 km s~', with an rms deviation of 
~ ±5.1kms-i. 
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Table 10 

Derived shiets eor FUSE data 



FUSE 




SllS 


SllS 


SllS 






Shirts 




Root 


Star 


dataset 


ID* 


grating 


FPALXPOS 


(pixs) 


(A) 


(km s ) 


B07101 


HD001383 


O5C07C010 


8241 


E140H 


117.0 


1.7 


0.035 


10.1 


B07104 


HD091983 


O5C08N010 


8241 


E140H 


117.0 


2.1 


0.043 


12.5 


P10165 


HD 168941 


O6LZ81010 


9434 


E140M 


117.0 


0.8 


0.016 


4.7 


P10179 


BD+354258 


O6LZ89010 


9434 


E140M 


117.0 


1.8 


0.037 


10.7 


P10188 


HD218915 


O57R05010 


7270 


E140H 


117.0 


0.5 


0.010 


3.0 


P10285 


HD190918 


O6359J010 


8662 


E140M 


117.0 


1.5 


0.031 


8.9 


P10286 


HD 192035 


O6359K010 


8662 


E140M 


117.0 


1.5 


0.031 


8.9 


P11631 


HD210839 


054304010 


8043 


E140H 


117.0 


2.7 


0.055 


16.0 


P12231 


HD210809 


O5C01V010 


8241 


E140H 


117.0 


2.3 


0.047 


13.6 


Z90101 


HD000108 


O5LH01010-40 


8484 


E140H 


117.0 


4.3 


0.088 


25.5 


B03001 


HD203374A 


O5LH08010-30 


8484 


E140H 


117.2 


1.3 


0.027 


7.7 


B03004 


HD208440 


O5C06M010 


8241 


E140H 


117.2 


1.3 


0.027 


7.7 


B03005 


HD209339 


O5LH0B020 


8484 


E140H 


117.2 


2.0 


0.041 


11.9 


B07105 


HD122879 


O5LH07020 


8484 


E140H 


117.2 


-0.5 


-0.010 


-3.0 


B07109 


HD206773 


O5C04T010 


8241 


E140H 


117.2 


1.5 


0.031 


8.9 


P10231 


HD091651 


O6LZ34010 


9434 


E140M 


117.2 


0.4 


0.008 


2.4 


PI 0232 


HD092554 


O6LZ36010 


9434 


E140M 


117.2 


1.8 


0.037 


10.7 


P10265 


HD134411 


O6LZ60010 


9434 


E140M 


117.2 


-0.2 


-0.004 


-1.2 


P10266 


HD151805 


O6LZ63010 


9434 


E140M 


117.2 


0.3 


0.006 


1.8 


PI 0279 


HD165955 


063599010 


8662 


E140M 


117.2 


-1.1 


-0.022 


-6.5 


P10289 


HD202347 


O5G301010 


8402 


E140H 


117.2 


0.9 


0.018 


5.4 


PI 1624 


HD192639 


O5C08T010 


8241 


E140H 


117.2 


0.9 


0.018 


5.3 


PI 1627 


HD206267 


O5LH09010-20 


8484 


E140H 


117.2 


1.9 


0.039 


11.3 


P12216 


HDE303308 


O4QX04020 


7301 


E140H 


117.2 


3.2 


0.065 


19.0 


PI 2230 


HD201345 


O5C050010 


8241 


E140H 


117.2 


2.6 


0.053 


15.4 


P12232 


BD+532820 


O6359Q010 


8662 


E140M 


117.2 


1.5 


0.031 


8.9 


P10275 


HD157857 


O5C04D010 


8241 


E140H 


117.5 


2.5 


0.051 


14.8 


P10123 


HD088115 


054305010-30 


8043 


E140H 


134.9 


2.4 


0.049 


14.3 


PI 0257 




O'S7R01fll0 


7270 


E140H 


134.9 


3.0 


0.061 


17.8 


PI 0262 


HD124314 


054307010 


8043 


E140H 


134.9 


1.9 


0.039 


11.3 


P10118 


HD066788 


O6LZ26010 


9434 


E140M 


135.1 


0.8 


0.016 


4.7 


P10221 


HD063005 


063531010 


8662 


E140M 


135.1 


0.6 


0.012 


3.6 


P10230 


HD091597 


O6LZ33010 


9434 


E140M 


135.1 


2.2 


0.045 


13.0 


P10236 


HD093205 


O4QX01020 


7301 


E140H 


135.1 


2.8 


0.057 


16.6 


PI 0237 


HD093222 


O4QX02020 


7301 


E140H 


135.1 


2.0 


0.041 


11.9 


P10241 


HD094493 


054306010 


8043 


E140H 


135.1 


1.6 


0.033 


9.5 


P10245 


HD099857 


054301010-30 


8043 


E140H 


135.1 


4.2 


0.086 


24.9 


P10246 


HD099890 


O6LZ45010 


9434 


E140M 


135.1 


0.7 


0.014 


4.2 


P10258 


HDl 14441 


O6LZ53010 


9434 


E140M 


135.1 


0.1 


0.002 


0.6 


P10261 


HD116781 


O5LH050 10-20 


8484 


E140H 


135.1 


0.9 


0.018 


5.3 


P12215 


CPD-592603 


O40P01D6Q 


7137 


E140H 


135.1 


1.3 


0.027 


7.7 


P12219 


HDE308813 


063559010 


8662 


E140M 


135.1 


0.9 


0.018 


5.3 


PI 0227 


HD075309 


O5C05B010 


8241 


E140H 


135.4 


1.6 


0.033 


9.5 


P10240 


HD093843 


O5LH04010-20 


8484 


E140H 


135.4 


1.5 


0.031 


8.9 



Table 10 — Continued 



FUSE 




STIS 


STIS 


STIS 






Shifts" 




Root 


Star 


dataset 


ID'' 


grating 


FPALXPOS 


(pixs) 


(A) 


(km s"') 


P10203 


HD013268 


063506010 


8662 


E140M 


175.1 


2.2 


0.045 


13.0 


P10205 


HD014434 


063508010 


8662 


E140M 


175.1 


1.3 


0.027 


7.7 


PI 0206 


HD015137 


O5LH020 10-40 


8484 


E140H 


175.1 


1.4 


0.029 


8.3 


P12201 


HDE232522 


O5C08J010 


8241 


E140H 


175.1 


2.3 


0.047 


13.7 


S3040201 


HD224151 


054308010 


8043 


E140H 


175.1 


1.6 


0.033 


9.5 


P10202 


HD012323 


063505010 


8662 


E140M 


175.4 


3.7 


0.075 


21.9 


P10204 


HD013745 


O6LZ05010 


9434 


E140M 


175.4 


1.9 


0.039 


11.3 


S3040202 


HD224151 


054308010 


8043 


E140H 


227.0 


-4.0 


-0.082 


-23.7 


PI 1623 


HD185418 


O5C01Q010 


8241 


E140H 


227.2 


-5.3 


-0.108 


-31.5 


P12241 


HD224151 


054308010 


8043 


E140H 


227.2 


-3.4 


-0.069 


-20.2 



""This is the shift which must be added to the FUSE data to aUgn the P(3) and R(4) 6-0 H2 lines with the CI I A1347 
line in the STIS data 

''Data from GO programs: 8241— "A SNAPSHOT survey of Interstellar Absorption Lines", Lauroesch PI; 8662, 
9434— "A SNAPSHOT survey of the Hot Interstellar Medium", Lauroesch PI; 8043, 8484 — "Physical Properties of 
H I and H II Regions", Jenkins PI. 
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B. STELLAR DISTANCES, THEIR ERRORS, AND CORRECTIONS TO AVAILABLE DATASETS 



A key set of parameters for our FUSE sample is the distance to each star, d, and its uncertainty. These quantities 
are important because one of our objectives was to measure as accurately as possible the mean and dispersion of the 
average density of O VI using all the available sight lines. To derive d in pc we need to calculate the spectroscopic 
distance (often referred to as the "spectroscopic parallax") given by 



5\ogd = imv-Av-Mv + 5) (Bl) 

where niy is the observed apparent magnitude in the V band. Ay is the interstellar extinction in the V-band, and My 
is the absolute V-band magnitude of the star of a particular spectral type and luminosity class (hereafter, 'Sp/L'). 
In the following subsections, we discuss how we derived each of these quantities. We show that it is possible to 
derive observed magnitudes and interstellar extinctions quite accurately using Hipparcos satellite data and infrared 
magnitudes obtained from the ground. We also show that the absolute magnitude of a star of known Sp/L is probably 
known to ^ 0.2-0.3 mags. What remains poorly known, however, is how well an Sp/L can be measured for a star, 
particularly for stars later in type than BO. 



B.l. Observed Magnitudes; Multiple Systems & Intrinsic Variability 



In order to use a set of magnitudes obtained in a uniform way and at a single epoch, we took stellar magnitudes 
from the Tycho Starmapper ( |ESA|[T997a| l catalogs. This instrument, part of the ESA Hipparcos mission, provides 
two-color photometry, Bj and Vj, in a magnitude system close to that of Johnson B and V bands (van Leeuwen et al. 



1997 1. The two magnitude bands are related via the transformation 



V = Vt-0.090(B-V)t 
(B-V) = 0.850(B-V)r 



(B2) 
(B3) 



provided -0.2 < (B-V)t < 1.8 (Eqn. L3.20 of 



ESA 



colors are given by equations L3.21 and L3.22 of |ESA ([ 



1997b 



The errors in the V magnitudes and in the B-V 



997b I and are reproduced here for convenience: 



cr(B-V) = G X cr(B-V)r 



= (1.09<+0.09<) 



1/2 



G X (al^+al^) 



1/2 



(B4) 
(B5) 



The factor G in the above equation provides the slope of the local (B - V) against (B - V)t relation and ranges from 
0.85 to 0.97. We calculate a precise value of G for each value of {B-V)t by extrapolating the values given in Table 
1.3.4 of ( ESA||T9"97b I. Using this scheme, the majority of the stars in our sample has B and V magnitude errors of 
< 0.04 mags. 

We note the following exceptions where we did not use these relationships: HD 326329 only has an entry for Vj 



(=8.90) in the Tycho catalog, with no Bj- For this star we used V an d B from Baume et al. ( 1999| ). HD 152314 is 
also not in the Tycho catalog, but does exist in the Tycho-2 catalog ( H0g et al.|2000 1, which lists both Bj, Vt and 
associated errors. Several stars have unusually large errors in the Tycho catalog, and their values are replaced by values 
from the Uterature. These stars are: CPD-417712 (Baume et al.|1999j ), and CPD-592600 & CPD-592603 ( |Massey &| 
|Johnson|1993] l. All the adopted magnitudes are listed in Table[T[ 

B.1.1. Photometric Binaries 

It is well known that stars are often members of multiple stellar systems. This presents the simple problem that 
a star's measured magnitude may not be a true measure of its distance if the measured flux really comes from two 



^"The "Guide to the Data" can be found at 
2'VizieR catalog 1/259 



http://www.rssd.esa.int/Hipparcos/CATALOGUE_VOLl/catalog_voll.html 
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(or more) unresolved stars. The Hipparcos data can again be used — at least for a subset of stars — to help identify 
multiple systems. 

Of the 148 stars in our sample, 106 are listed in the Hipparcos catalog and could therefore be analyzed in more 
detail than was possible from their Tycho data alone. Of these, 19 are listed in the Hipparcos Double and Multiples 
System Annex: Component solutions (DMSA/C) catalog (JESA 1997d I, which means that these system can be modeled 



as individual stars. (The Hipparcos catalogs labels these as "Component Solutions".) For a particular system, we adopt 
the magnitudes of the brightest component listed in the DMSA/C catalog; we thereby assume that the bulk of the flux 
which defines the Sp/L of the star comes from the brightest component. Six stars (HDs 099857, 101190, 101413, 
101436, 190918, 152248) are cataloged in the DMSA/C catalog, but the magnitude of only one component is given. 
We take this to mean that the secondary is too faint to be modeled; its presence is therefore not expected to alter the 
magnitude of the primary star. (We use the DMSA/C catalog magnitude nonetheless.) 

In the DMSA/C catalog, only Hipparcos magnitudes Hp are available, and not the Vt used to convert to Johnson V 
magnitudes for the other stars. A comparison between Hp and Johnson V for stars not labelled as multiple show that 
in all cases, differences in V (calculated from Vj) and Hp are small (< 0.05 mags). The transformations suggested by 
|Harmane^ ( | 1 998| l also show that for B-V = ±0.2 — the minimum and maximum values for our stars — the difference 
in Hp and V should only be =f0.06 mags. 

Two of the 106 stars (HD 185418 and HDE 232522) are listed in the Hipparcos catalog as being multiple due to 
their non-linear motion ("Acceleration Solutions" — similar to astrometric binaries) although they have no entry in 
the DMSA/C catalog. Both are actually spectroscopic binaries, and their V magnitudes are corrected for in same way 
as the rest of the spectroscopic binaries in our sample, discussed below. 

For the 42 stars not in the Hipparcos catalog, the Tycho catalog provides a flag to indicate which stars may be 
multiple system^ None of the 42 stars not in the Hipparcos catalog were flagged as multiples in the Tycho catalog. 

B.1.2. Spectroscopic Binaries 

An important class of multiple systems are the spectroscopic binaries (SBs). These stars are too close together to be 
resolved by Hipparcos (except for those which might also be eclipsing binaries — see below). Given the luminosity 
of the primary and secondary stars to be Li and L2, the correct distance to the system d' is 



d'=[\ + -^j d (B6) 

where d is the distance to the system estimated from the integrated (Li +L2) luminosity. In the worst case, L\ and L2 
are equal and d is incorrect by a factor ^/2. If the magnitudes of the primary and secondary, m\ and m2 respectively, are 
known, the ratio of the luminosities can be calculated in the usual way La/^i = lO*'^*'"'"'"^^ In many cases, however, 
the value of Am is unavailable. How then can we estimate what correction to make for these systems? 

For double-lined SBs, the ratio of the radial velocities for each component gives the mass ratio of the system 



M2/M1. The mass-luminosity relation for main sequence (MS) stars with mass > O.2M0 is given by Schmidt- Kaler 
([TM2}: 

log-^ = 3.81og-^-0.8. (B7) 
To convert between the bolometric luminosity L, and V-band luminosity Ly, we must apply a bolometric correction 



(BC) to equation B7 The BC is a function of luminosity, but can be derived for a given Sp/L; for MS stars, the BC 
can be shown to be -0.9 log(L/L0) + 1.2 ( |Schmidt-Kaler| 1 982 1 which gives 



-^Field T36 of the Tycho catalog, Source of photometric data; see ESA ■ 1997c 
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where the luminosities are those measured in the V-band. This can then be appHed to equation B6 to correct the 
distance to the system. 

To identify spectroscopic binaries in our sample of stars, we used The Ninth Catalogue of Spectroscopic Binary 
Orbits (SB9) by |Pourbaix et al.| ( |2004| , and the list of binaries in NGC 6231 by |Garcia & Mermilliodl ( |2001[ ), who note 



that HD 152314 (S161), HD 152200 (S266), HD 152233 (S306) and HD 152248 (S291) are probably spectroscopic 
binaries. We identify 21 of our stars to be SBs. For 11 of these stars, values of AL or AM are available, and their 
distances can be corrected directly from equation |B 6 1 or indirectly from equation |B 8 1 For the remaining ten stars, we 
note that iPinsonneault & Stanek (J2006) find that 50% of binaries have mass ratios > 0.87, and use this as a generic 



correction. This mass ratio corresponds to AL = 0.72, and hence the distance to the system is increased by a factor of 
1.3. A superscript next to the distance to a star in Table[T]indicates that the star was treated as a SB. 



As noted above, equation B7 is probably only applicable for MS stars, and evolved stars may well not obey the 
same M-L relationship. Unfortunately, only seven of our SBs are MS stars. It may be that equation |B7] does not 
represent a precise conversion between mass and luminosity, and, therefore, that the correction we make to the stellar 
distance is similarly imprecise. The distances to these systems should be treated with caution. 

B.1.3. Variable Stars 

Stars which are not multiple systems may still exhibit variability. These may be rotating, pulsating, eruptive or 
cataclysmic types of variability. B-type emission line stars may also show variations in luminosity. 14 stars in our 
sample which are not in SB9 or in the Hipparcos Double and Multiples catalog are listed in the Hipparcos Variability 
Annex ( ESA|[T997e l. For these stars, we take the maximum Hp magnitude (minimum luminosity) from the fitted 



Ught curve derived from the Hipparcos photometry. The assumption here is that the minimum luminosity of the 
star represents its 'true' luminosity for its given spectral type and luminosity class. This assumption is based on the 
idea that Be stars are brighter than B stars for a given spectral type ( |Briot et aL]|1997| l; obviously, for stars whose 
luminosities are not increased by mechanisms such as the emission from circumstellar disks which are thought to 
surround Be stars, this assumption may be invalid. However, of the 14 stars corrected, all but one have a difference 
between the faintest magnitude and the mean magnitude of less than 0.06 mags. (HD 042401 shows a difference of 
0.17 mags). Although these corrections are therefore largely irrelevant, they have been applied to our data. Distances 
adjusted this way are noted with a superscript in Table [T] 

Variable stars may, of course, be (partially) eclipsing binaries, more akin to the SBs discussed above. However, 
most of the stars in the Variability Annex which are labelled as likely (or already known) eclipsing binaries (field 
'P5' in the Annex) are already identified in SB9, and their distances are corrected in the way described above. (Only 
HD 042401 is identified as a /? Lyrae type of eclipsing binary and is not in SB9.) Hence the luminosities for these stars 
are not corrected in the same way that the spectroscopic binaries were corrected above. 

Unfortunately, for the 42 stars of our sample not in the Hipparcos catalog, the Tycho analysis provides no clear 
information on which stars are intrinsically variable. Although a variability flag does exist, the flag appears unreliable 
(failing to note variability when detected by Hipparcos, and flagging variability when Hipparcos detects none) and we 
have not included the information in Table [T] Instead, we searched for our stars in the Combined Ge neral Catalog ue 
of Variable Stars (GCVS) and the GCVS Suspected Variable stars and Supplement catalogs (Samus~et al.||2004p ^ 



This catalog contains a very inhomogeneous set of data, but can indicate that a star was recorded to have a fainter 
magnitude than was seen during the ~ 3 year Hipparcos mission. 

For stars which have no entries in the Hipparcos catalogs, are not SBs, or which have distances that are not taken 
from the literature, we find that: a) only one has a minimum-maximum variation listed in the GCVS catalogs greater 
than 0.06 mags | |Lipunova & Putilina| (|1985) noted that HDE 332407 fell from V = 8.51 (close to the Hipparcos value 
of y = 5.82) by ~ 0.1 mags. Hipparcos found no variation, but we have made this correction in Table [T]; and b) six 
stars have variations less than 0.06 mags[^, which we consider small enough to ignore. 



-^vizieR catalog 11/250 

^■^Using a value of 0.06 mags to separate 'weak' from 'strong' variability is purely arbitrary, but the value is used in the Hipparcos catalog to provide 
a course variability flag (Field H6 of the main catalog), which is why we use that same value here. Stars in our sample which show a variability 
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Given the short Ufetime of the Hipparcos mission, it is also interesting to consider whether the maximum magni- 
tudes of the variable stars that are indeed in the Variability Annex are significantly different from the maxima given 
by the GCVS catalogs. In fact, for the six stars for which there is a difference > 0.06 mags (excluding SBs), the 
distances we use are taken from the literature, measured by Hipparcos parallax, or the GCVS magnitude is not from 
V-band measurements. There is therefore no evidence that our distances are affected by not considering a long enough 
timeline for the variable starf^ 



B.2. Interstellar Extinction 



The relationship between interstellar extinction Ay and the color excess of stars E{B — V) is usually taken to be 
Ay = Rv E(B - V), where Ry is the ratio of total to selective visual extinction. In the Milky Way, a typical value of 
Ry through the diffuse ISM is ~ 3.1. In a compilation of Ry values for 417 03- to B5-type stars, Valencic et al. 
( |2004 ) found a Gaussian distribution of Ry centered at 3.0 (largely as expected), but with a la dispersion of 0.35. 
Assuming that a large fraction of this dispersion arises from real line-of-sight differences in Ry, and not just from 
measurement errors, then using a value of 3.1 for Ry to calculate Ay might be lead to errors of at least 0.35£'(B-y) 
mags in calculating the distance to a star in equation |B 1 [ even if the errors in E{B-V) were zero. For values of Ry 
very different from ~ 3 the errors are even larger. 

We would therefore like to have a way to measure Ay along a sight line directly. Whittet ( 1992 1 describes a method 
which uses infrared (IR) magnitudes: Ry ~ I.IE(V-K)/E{B-V). Since Ay = /?v'£(fi-V'), we can directly measure 
the extinction from 



Ay = l.l E{V-K) 



(B9) 



P6 



The 



To obtain A'-band magnitudes, we used the Two Micron All Sky Survey (2MASS) ( Skrutskie et al.|[2006f 
2MASS magnitude was converted to Johnson K using the relatiorj^/T = /r2MASS + 0. 027 + 0. 00787(72mass-^2mass) 
unless the 72mass was flagged as unreliable, in which case the color term was neglected. Values of ^r2MASS were only 
used if an 'A' photometric quality flag (Qflg) was given (see the 2MASS catalog for a full explanation of this flag). 



Intrinsic IR colors for early type stars were taken from Wegner ( 1994) . 

We calculated the error in the IR color excess as aiEy-^)^ = aiV)^ + a(K)^ + 0.03^, where the last term is adopted 



as a typical error in the intrinsic V-K colors given by Wegner ( 1994 1; (V-K\) appears to vary slowly with luminosity 
class, but strongly with stellar type. With this error, we were able to reject poorly determined IR color excesses 



when E(V -K) < 2a{Ey-x), and did not use equation B9 E{V -K) was also not used if more than one stellar type 
or luminosity class was given for a particular star, or if the star was in the Hipparcos Variability Annex or DMSA/C 
catalogs, since we have no information on how a star might vary in the IR from 2MASS. We allowed E{V-K) to be 
used for the SBs; here we assume that AL in the optical is similar to AL in the IR. 

If we were unable to use E(V -K) to estimate Ay, we fell back to the more usual practice of calculating Ay = 
3.1 E(B-V). E(B-V) is given for most of our stars in TablefT] except when distances were taken from the literature. 



The error in the Johnson B-V color is given by equation |B5| above; we again used Wegner ( 1994 1 for the intrinsic 
(B-V) colors of the stars, and assumed a value of ~ 0.01 mags as a likely error, since the intrinsic color varies very 
slowly with spectral type and luminosity class for early type stars. 

B.3. Absolute Magnitudes of O-type and Early B-type Stars 

Beyond the problems in measuring the observed magnitudes of, and the extinction towards, a particular star, the 
largest uncertainty in measuring its distance comes from a poor knowledge of its absolute magnitude My in equa- 
tion |BT] There are two problems: knowing the absolute magnitude of a star given its MK-system spectral type and 



< 0.06 mags usually do not show up in the Variability Annex. Hence we cannot coirect the magnitudes of these stars anyway. 
^^Such variations may exist of course, but they have not been recorded in the literature. 

^''vizieR catalog 11/246 

27 



From the Explanatory Supplement to the 2MASS All Sky Data Release at http://www.ipac.caltech.edu/ 2mass/releases/allsky/doc/sec6_4b.html by 
R. M. Cutri et al.. 
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03 07 09 BO Bl B3 

Spectral type for luminosity class V 

Fig. 21. — : Plot of absolute V-band magnitude My against spectral type (03-B3) from the authors listed at the top 
left of the figure. Our adopted relationship between My and spectral type is shown as a solid black line; the gray region 
either side marks a difference of ±0.25 mags which we take as the error in knowing My for any given stellar type. 



luminosity class [at least, the extension developed by Walborn ( 1971[ l, since the original MK classification by Mor- 



gan et al. ( 1943| l did not cover stars earlier than 09]; and, more importantly, knowing precisely the spectral type and 



luminosity class of any given star 

The first of these problems is largely tractable, thanks to surveys which use stars with well-defined Sp/Ls, carefully 
measured observed magnitudes, and distances measured using methods other than the spectroscopic parallax. Al- 
though the measurements of all these quantities are fraught with their own uncertainties, there appears to be consensus 
on how My varies with Sp/L for early type stars. To produce our own My-Sp/L scale, we collated data from Walborn 



1972) 1, IWalboriil ( [T973| ), |Leshl ( [T979] l, [Schmidt-Kaler| ( [T9821 ), [Humphreys & McEiroy] ^19S4\ , and |Howarth & Prinja 



1989|l. We also added, when an Sp/L was available, magnitudes from models by Vacca et al. ( |1996| l and Martins et al. 



( 2005| l for comparison. Not every available published scale was included: for example, absolute magnitudes for OB 
stars have been given by Wegner ( 2000[ l based on Hipparcos parallaxes, but these values differ so significantly from 
the other authors just cited that we did not include them. The final My magnitudes for a given Sp/L are tabulated in 
Table [n] 

An example of one My-S^fL scale is shown in Figure[2T| which plots My against spectral type for luminosity Class 
V stars. (These, along with Class III stars, are the most well studied of the luminosity classes.) In general, authors 
tend to agree on the relation between My and spectral type to within several tenths of a magnitude for this luminosity 
class. In Figure 21 we show our adopted scale as a straight line, along with an error of ±0.25 mags (shown by gray 
boundaries). Theoretical values are also in good agreement with the data: the models by |Martins et al.|(2005| l appear 
to be too faint by only w 0.25 mags, while those from Vacca et al. ( 1996 1 agree well with observations for the early 
stellar types until BO. 
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Figure 21 also shows the values of My derived by Cramer ( 1997 1 for B-type stars whose parallaxes were measured 
by Hipparcos. We did not use these data in setting ourMy-Sp/L scale as no O-type stars were included, and a complete 
analysis remains unpublished. The variation of My with Sp is not entirely smooth, presumably due to the small number 
of B-type stars close enough to the Sun for an accurate parallax to be determined. Nevertheless, although the stars as 
a group seem to be slightly fainter than our adopted scale, the values of My are largely consistent with our adopted 
scale, to within an error of « 0.25 mags. (We return to the question of how well our spectroscopic distances compare 
to Hipparcos parallaxes in Appendix |C.l. 3 1 ) 

Plots of My and spectral type for other luminosity classes show similar results, although the Class I stars have no 
'break' around BO but have roughly constant absolute magnitudes for all spectral types, as expected. 

One other problem which remains is that of the Be stars, those which presently, or in the past, have shown emission 
lines. There are nine Be stars in our sample; the adopted Sp/L for these should probably be used with caution, since 
the classification of Be stars is known to be problematic, due to emission lines filling-in the diagnostic absorption 
lines, and due to the large rotation speeds which are usually associated with these systems. Unpublished studies of 
the difference in My between B and Be stars using Hipparcos parallax distances by [Briot et al.| ( |1997| and Briot & 



Robichon ( 2000| suggests that although Be stars of a given Sp/L are indeed brighter than B stars of the same Sp/L, the 



differences for the stars of interest in this paper, B0-B3, are negligible. Hence we make no additional corrections to 
My for Be stars. 

Given the relative small dispersions in My and the general agreement with theoretical values (when available), we 
assume that My can be known for a star of a particular Sp/L to within ±0.25 mags. 

The second problem in determining My for a given star — that of knowing its Sp/L precisely — is not so straight- 
forward. Stellar classification depends not only on the quality, resolution, and wavelength extent of the spectra used 
to make the classification, but on the experience of the observer making the classification. (Many stars were assigned 
their Sp/Ls several decades before the development of automatic classification of digital spectra, for example.) More- 
over, although a particular set of stars can be used to define a class and type (e.g.l Garrison|1994| ), Sp/L categories are 
discrete, and any given star may in reality fall between the designated category. For example, the absolute magnitudes 
in Figure 21 decreases rapidly for stars between types Bl and B2, though only very small changes are seen in line 
ratios used to define the Sp/Ls. Hence a star may have a magnitude a few tenths different from the My assigned from 
the adopted Sp/L. More serious, however, is the fact that the use of a poor quality spectrum of a star can lead to a 
complete mis-classification of a star, an occurrence which appears more frequently for the B-type stars. 

In most cases, we have no access to the original data used to determine the Sp/L of the stars in our survey. We have 
therefore searched through the available literature to find suitable Sp/Ls for our stars. The references for the Sp/Ls we 
adopt are given in the final column (16) of Table [T] 

Given the difficulties which exist in assigning an Sp/L to a star, we have developed a simple procedure for measuring 
the likely range of plausible distances to an individual star. Our primary goal is to quantify a generic error in My for 
each star. Assigning only a single number for this error is inappropriate, since uncertainties in My are larger for, e.g. 
B-type stars than for O-type stars of the same luminosity class. Instead, we calculate an error in three parts: we first 
take the difference in magnitude between the adopted value and the magnitude of the next spectral type in the series; 
we then calculate the difference in magnitude between the adopted spectral type and that of the same spectral type in 
the next luminosity class. Finally, we add these two values in quadrature with the error we adopt for assigning My to 
a given Sp/L, which we take to be a(My)c = 0.25 mags, as discussed above. In notation form, we calculate an error 
in My for possible positive deviations from the adopted My thus: 



a{My,{Splh)if = AMy(Spi^^-Spf + AMy(Li^i-Lif + a(My)l (BIO) 

So, for example, the positive error in a B3 111 star is formed from taking Mv(B4 III) - My{B3 III) in Table 11 
Mv(B3 IV) -Mv(B3 III), and adding these in quadrature, along with 0.25 mags for (j{Mv)c- The negative deviation is 
calculated in the same way, only taking differences between Sp, and Sp,-i, and L, and L,_i. Using these values of My 
leads to an upper and lower bounds for the distance, du and di, respectively. 

Even this method is imprecise however, since its is easier to classify, e.g. an 04 star than a Bl star. In this sense. 
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the ranges in My determined for O type stars may be somewhat over-estimated, while for B type stars the errors may 
be too small. However, the magnitude of the errors are likely to be roughly correct, since in the scheme outlined in 
equation [BlOj O stars have smaller errors than B stars. 

From estimating distance 'errors' in this way, we conclude that the errors we derive in My dominate all other errors 
discussed in this Appendix. For example, for the Class V dwarfs and Class IV subgiants, errors range from ~ 0.3-1.0 
mags for 03-B3 type stars, or a 14-46% error in the distanc^^ The largest errors in this scheme occur for the B3 
Class II giants, since their magnitudes are furthest from both the la supergiants as well as the Class III giants; errors 
range from 0.4- 1 .5 mags for 03-B3 type stars, 20-70 % in distance. Clearly, the errors associated with the extinction 
Ay (— 0.1-0.2 mags) and those with the Tycho and Hipparcos catalog magnitudes (<< 0.1 mags), are insignificant 
compared to the errors in My- 

Further, in some cases, the stars in our sample have no published luminosity class, or only a range of luminosity 
classes are given. In these cases we calculate the distances for the brightest and faintest of the specified classes (or 'V 
and Ta' if no Class is given), and use the average value for the distance. The minimum and maximum error are simply 
the minimum and maximum values calculated for each class in the way described above. The value of E{B - V) in 
Table[T]is also an average value (although the average was not used in the calculation of the distance.) 

B.4. Distances from the Literature 



There are some exceptions to this scheme for calculating the distance to our survey stars. Our sample contains four 
Wolf-Rayet stars; distances and values QiE{B-V) for three of these, HD 151932 (WR78), HD 191765 (WR 134) and 



HD 190918 (WR 133) are taken from the catalog of [van derHucht ( 2001 1, while data for HD 187282 (WR 128) come 
from |Conti & Vacca| ( |T990| . HD 93129A is an 02 If star ( jWalbom et al.|2002b| in Trumpler 14, one of the ionizing 
clusters in Carina, and since 02 stars have no defined values of My, we adopt the distance given by Tapia et al. (2003 1. 
The Hipparcos DMSA/C catalog lists four components to the multiple system of which HD 005005A is a member, 
but their membership of the IC 1590 cluster seems unequivocal — unlike the Sp/L of HD 005005A — and an accurate 
distance is given by Guetter & Turner| ( |1997) , which we use. LS 277 has no published Sp/L, and its distance is taken 
from Reed ( 1993 1. Finally, HD 1 10432 is close enough that the parallax distance from Hipparcos (300 ±50 pc) should 



be more accurate than its distance determined from its Sp/L (see Appendix C. 1 



Note that in Table [T] use of a distance from the literature is indicated with a superscript next to the distance. If no 
errors for a published distance are available, we adopt a value of ±20 % of the distance. 



'u{d)ld = OA(,u(Mv) 
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Table 1 1 

Adopted My mags for given Spectral Types and Classes 



Luminosity Class 



Type 


V 


IV 


III 


II 


lb 


lab 


la 


03 


-5.80 


-6.00 


-6.30 


-5.95 


-6.30 


-6.65 


-6.90 


04 


-5.65 


-5.88 


-6.20 


-5.94 


-6.29 


-6.64 


-6.91 


05 


-5.50 


-5.76 


-6.10 


-5.93 


-6.28 


-6.63 


-6.92 


05.5 


-5.35 


-5.64 


-6.00 


-5.92 


-6.27 


-6.62 


-6.93 


06 


-5.20 


-5.52 


-5.90 


-5.91 


-6.26 


-6.61 


-6.94 


06.5 


-5.05 


-5.40 


-5.80 


-5.90 


-6.25 


-6.60 


-6.95 


07 


-4.90 


-5.28 


-5.70 


-5.89 


-6.24 


-6.59 


-6.96 


07.5 


-4.75 


-5.16 


-5.60 


-5.88 


-6.23 


-6.58 


-6.97 


08 


-4.60 


-5.04 


-5.50 


-5.87 


-6.22 


-6.57 


-6.98 


08.5 


-4.45 


-4.92 


-5.40 


-5.86 


-6.21 


-6.56 


-6.99 


09 


-4.30 


-4.80 


-5.30 


-5.85 


-6.20 


-6.55 


-7.00 


09.5 


-4.15 


-4.68 


-5.20 


-5.84 


-6.19 


-6.54 


-7.01 


09.7 


-4.00 


-4.56 


-5.10 


-5.83 


-6.18 


-6.53 


-7.02 


BO 


-3.85 


-4.44 


-5.00 


-5.82 


-6.17 


-6.52 


-7.03 


B0.2 


-3.70 


-4.32 


-4.90 


-5.70 


-6.10 


-6.51 


-7.04 


B0.5 


-3.55 


-4.20 


-4.80 


-5.50 


-6.05 


-6.50 


-7.05 


B0.7 


-3.42 


-4.08 


-4.60 


-5.30 


-6.00 


-6.49 


-7.06 


BI 


-3.04 


-3.96 


-4.25 


-5.10 


-5.95 


-6.48 


-7.07 


BI.5 


-2.66 


-3.50 


-3.90 


-4.90 


-5.90 


-6.47 


-7.08 


B2 


-2.28 


-3.10 


-3.55 


-4.70 


-5.85 


-6.46 


-7.09 


B2.5 


-1.90 


-2.70 


-3.20 


-4.50 


-5.80 


-6.45 


-7.10 


B3 


-1.52 


-2.30 


-2.85 


-4.30 


-5.75 


-6.44 


-7.11 


Refs:" 


1-7 


1,3-5 


1-7 


1,2,4,5 


1,2,4,5 


1,2,4,5 


1,2,4,5,7 



"Collation of data from: (1)— }Walbom|{l972| ; (2)-{Walborn ( 1973|; (3)— |Lesh| 
1979)); (4)— |Schmidt-Kaler| ( [T982); (5)— [Humphreys & McElroy] ( [T9g4| l; (6)— 
IHowarth & Prinja| ( |1989t ; (7 HVacca et al.H1996> 
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C. REVISITING EARLIER SURVEYS 



C.l. Corrections Made to Copernicus Data 

C.1.1. Notes on Stellar Parameters 



In §3.1.2 we introduced the stars studied by J78. Since they are both bright and comparatively nearby, there is a 
considerable amount of data available in the literature for these stars. In deriving distances to the Copernicus stars, 
we have tried to treat them in exactly the same way as our FUSE sample. All the Copernicus stars are listed in the 
Hipparcos main catalog, which means that they can be treated even more uniformly than the FUSE sample of stars. 



The stellar data for the Copernicus stars are given in Table 12 below, we outline some of the changes we have made 
to the stellar data since the stars were studied by J78. 

Sp/Ls: these are taken directly from Table 1 of J78. Many of these stars have been studied since the advent of stellar 
spectroscopy, and their Sp/Ls are well understood (and often actually define a particular Sp/L). Small changes have 
been made to the Sp/Ls of 7 Cas, e Per, C, Ori, C, Pup, (3 Cen, (3 Lup, 5 Sco, and <j Sgr since improved classifications 
are available in |Garrison & Beattie| ( |1996| . 

Magnitudes: V and B magnitudes are converted from Bj and Vj given in the Tycho catalog in the same way as 
described above. We have again tried to use E(y -K) to calculate Ay, but we find that many of the stars are saturated 
in the 2MASS data. 

Emission line stars: more information is now available for whether emission lines have been detected in the stellar 
spectra of the Copernicus stars, and we have appended 



suffixes to the Sp/Ls in Table 12 for the following stars: 
a Eri ( [Hanuschik et al.|1996[ ), 77 Ta u ([Tycner et al.|2005| l, a Cam ( [Morel et al.|2004| l, A Eri ( Hanuschik et al..l996), 
25_Ori ( [Banerjee et al.|2000^ e Ori ([Morel et al |2004pID 051283 (|Irvine|1990| ), HP 052918 (|Cote & van Ke'rkwijk 



Kerkwijk|1993l ), 60 Cyg (Koubsky et al. 200 0), and 8 Lac ( Chauville et al.|2001| l. 



19931 ), ^ Sco ( [Galazutdin ov & Krelowski''2006\ A Pav (■Mennickent"1991j), 2 Vul ( |Barker|1983| l, k Aql ( [Cote & van 



Spectroscopic binaries: we have again used the SB9 to flag which stars are SBs, and have adopted values of AL or 



AM when available (Appendix B.1.2 1. For HD 057060 (29 CMa), AL and its Sp/L are taken directly from Bagnuolo 



et al.[ ( [1994| l. HD 037303 is cited by MorreU & Levato[ ( [T99T) as being an SB in the Orion OBI association, though 
the star is not listed in SB9. We assume this designation is correct. For HD 047839 (15 Mon) Hipparcos gives the 
magnitude of the primary component {V = Hp = 4.61), and we take Am = 1 mag between primary and secondary from 



Gies et al. ( 1993 1. HD 158926 (A Sco) is composed of two early-type B stars and a low-mass pre-main-sequence star. 



according to Uytterhoeven et al. ( 2004[ l, but Hipparcos finds no multiplicity in the DMSA/C catalog, and the system 
is not listed in SB9. We use these authors' estimate for the mass ratio of the two stars, AM ~ 0.84. 

Variables: in general, stars flagged by the GCVS as variable are also found in the Hipparcos catalog. There is 
Uttle published evidence that any of the Copernicus stars have been substantially fainter over the last century than the 
faintest magnitudes found by Hipparcos. We know of two clear exceptions: the long period V-band light curve for 
HD 005394 (7 Cas) has a minimum luminosity of V = 2.75 ( [Doazan et al.|198 3D, while HD 200120 is known to show 
a minimum of at least V = 5.05 ( Harmanec et al.|2002 1. We adopt these magnitudes as the 'true' values for these two 
stars. 

Multiple systems: HD 36486 {5 Ori) is a triple system. Hipparcos resolves the eclipsing binary, but the primary 
is an SB. Harvin et al. ( 2002] l show that the flux from the SB is dominated by the 09.511 primary, and that AL is 
small, only 0.09. We use the magnitude of component Aa, and correct that using AL = 0.09. For HD 57061 (r CMa), 
Hipparcos finds two components; using the magnitude of the primary, and assuming this is composed of two stars, 
we obtain a distance of 1224 pc. However, the system appears to be quite complicated [see the discussions by van] 
[Leeuwen & van Genderen ( [1997[ ) and Stickland et al. ( 1998| l]. Although we could use the distance to NGC 2362, [e.g. 
1480 pc by Moitinho et al.[ ( 2001| l], [van Leeuwen & van Genderen[ ( [ 1997 1) note that HD 57061 may no longer be near 
the center of the cluster. Given the lack of information, we simply retain the derived distance of 1224 pc. 
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C.1.2. Errors on Copernicus O VI Column Densities 



The Copernicus stars tend to probe distances in the Galactic disk in the range between the WDs and the stars studied 
in the present survey. As such, their inclusion in an analysis of how, e.g., A^(0 VI) varies with distance, is valuable. 
Unfortunately, no errors on the column density a(N) were given by either J78 or | Jenkins & Meloy ( 1974it. Yet a{N) 
for these stars is an important quantity, since we need to know how much statistical weight to give A^(0 VI) when 
deriving quantities such as the average volume density in the Milky Way disk (as discussed in ^3.8[ ). 

We can estimate these errors by using the coadded scan data of the stars from the Copernicus MAST archive. These 
data include a 2a noise error array, which can be used to measure a{N). 

Deriving A^(0 VI) from Copernicus MAST archive data is complicated; estimating the continuum of the star around 
the O VI line is often difficult, and the precise subtraction of the HD absorption (when other HD lines are available to 
produce a satisfactory model of the absorption) is challenging. We believe we are unlikely to improve upon the actual 
values of A^(0 VI) given by J78, and so concentrate on how best to obtain a(N). The Copernicus data are of variable 
quality, so a fixed percentage (for example) of the measured A^(0 VI) would be an inappropriate measure of the error 
in the column density. Since we only require approximate estimates of (t{N), our approach is to produce a very crude 
re-measurement of A^(0 VI) and extract a value of a(N) associated with that column density. To first order, this value 
of a(N) should be close to the real value. The values are likely to be less precise than those derived for the FUSE 
sample, but they can be used to correctly weight the Copernicus values when data from all samples are combined. 

We fitted the available Copernicus coadded scans made around the O VI A 1032 line in the same way as described 
for the FUSE sample in 5 2.6 and measured A^(0 VI) with the AOD method (i 2.8.2 1. The velocity range used by J78 to 
measure A^(0 VI) was variable, anywhere between [-50,+50] km s"' to twice that range. For simplicity, we measured 
a(N) between -70 and +70 km s"' . We calculated a final value of (7{N) by adding in quadrature the errors from the 
continuum uncertainties and the noise errors given by the coadded scan error arrays. 

Despite the simplicity of this approach, our measurements of A^(0 VI) and those published by J78 are in reasonable 
agreement. More importantly, J78 plotted profiles of the optical depth at each pixel, r,, along with 2a error envelopes, 
2ct(t,), for about half of the Copernicus sample. This means we can directly compare the values of ct(t,) measured 
from the Archive data with those shown in Figure 2 of J78 for many of the stars. In general, the agreement was very 
good. Only three stars showed apparent problems: HD 064740, HD 106490 and HD 121263 had unusually large 
values of cr(T,) compared to the envelopes shown in J78, all because of inexplicably large noise errors cr, in the Archive 
data error arrays. We estimated that the errors were too large by approximately a factor of 4-5, and subsequently 
reduced all values of cr, by a factor of 5 before calculating cr(T,). 

Compared with the dataset as a whole, (t(t), for another three stars, HD 52918, HD 157246 and HD 158926, also 
seemed unreasonably large given the apparently good S/N of the data. These stars were originally discussed by Jenkms] 



& Meloy ( 1974 1 but neither they or J78 showed plots of r(v) against which we could compare the new values. (The 



original analysis of these stars is no longer available.) We therefore again reduced the values of cr, by a factor of five, 
to give errors in cr(T,) apparently commensurate with the S/N of the data. 

Values of a(N) are given in Table 12 [along with the original values of A^(0 VI) given by J78]. Again, we emphasize 
that these values of a{N) are only approximations and may only be useful for the statistical weights they provide for 
the Copernicus column densities. 



C.1.3. Comparison Between Sp/L and Hipparcos Parallax Distances for Copernicus Stars 

While we derived the distances to the Copernicus stars in exactly the same way as we did for FUSE stars, it is 
also true that reliable Hipparcos parallax distances are also available for many of the Copernicus stars. The parallax 
distance d^^ can be derived from the parallax tt given in the Hipparcos catalog from the relation d.^ = IOOO/tt pc, 
where tt is the parallax measured in mas, which has an error of cr(7r). 

Apart from providing a direct measure of the distance to a star, these data also enable us to compare distances with 
those derived from given Sp/Ls, and examine how good the My and color calibrations discussed above actually are. In 
Figure|22]we plot the Hipparcos distance against the distance derived from spectral typing. In the figure, we distinguish 
between those stars with ir/ain) > 5 (black circles) and the rest of the sample (gray circles). Within 100-200 pc. 
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Fig. 22. — : Comparison between parallax distances estimated from Hipparcos data and those derived from adopted 
spectral type and luminosity classes, for stars observed with Copernicus. One star in the sample (HD 149881) has a 
negative parallax, and is not plotted. 



the correlation between the distances derived in the two different ways is quite good, and the spectroscopic distances 
are generally equal to the Hipparcos distances to within at least twice our estimates of the error range determined 
by equation BIO This substantiates our use of the errors from equation BIO as a useful 'lu' estimate. Beyond 
200-300 pc, parallaxes become too small to be measured accurately by Hipparcos. We note that the Be stars do not 
have errors in spectroscopic distances much different from B stars. 

Given that the Hipparcos distances for the nearest stars have much smaller errors than the spectroscopic distances, 
we adopt the Hipparcos distance for any star with 7r/cr(7r) > 5. (We note our adoption of the Hipparcos parallax 



distance with a superscript 'c' against the value of d in column 11 of Table 12 ) For these stars, the errors in their 
distances should be < 5 % ( Brown et al.|1997 1. In keeping with the range in errors adopted for the FUSE targets, we 
take the range in possible distances, di and du, to be 1000/[7r + cr(7r)] and 1000/[7r-(T(7r)], respectively. 



C.2. Notes on Corrections Made to Data from Z03 



All the stars observed by Zsargo e t al.| ( |2003 1) can be found in the Hipparcos catalog, so we have again used the 
cataloged values for the stellar magnitudes. In many cases, however, Sp/Ls are either not known, are only poorly 
determined, or are irrelevant — as in the case of the small number of post-asymptotic giant branch (PAGB) stars. For 



these stars, we have taken distances directly from the literature. Details are given in Table 13 
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D. Upper Limits 

Astronomical observations sometimes exhibit measurement outcomes that are not significantly above the noise. 
Ideally, in such circumstances one should report the formal measurement and the associated error, as we have done 
in columns 3-7 of Table [5] This procedure conveys the most information. Nevertheless, it is a common practice to 
define a minimum measurement level above which one can claim a "detection" (with an associated uncertainty), and 
this comfort-level threshold is usually expressed as some factor times the standard deviation of the noise a. On the 
occasion that a measurement is found to fall below this threshold, it is customary to state that only an upper limit at 
the Na level of confidence for the true signal is warranted. In defining this upper limit, many investigators arbitrarily 
discard the actual measurement outcome and claim that the upper limit is simply equal to Na. When this is done, there 
is a complete loss of potentially useful information about whether the measurement came out at just barely below the 
threshold, or alternatively, below zero by some one or two standard deviations of the noise (which can happen with 
reasonable frequency for null signal levels in the absence of any noise). The former represents a weaker case for an 
upper limit than the latter, since we know that it is much more difficult for a noise deviation to drive any moderately 
positive signal to a negative outcome than to cancel it just enough to go somewhat below the defined threshold level. 
This difference is not reflected in the claimed degree of confidence for the upper limit. Moreover, there is an untidy 
discontinuity in the upper limits inferred from the "upper limit only" cases for outcomes just below the threshold, 
compared with "detections plus Na" for outcomes just above the threshold. For these reasons, we adopted a more 
refined procedure. 

In §2.8.2 we made use of a scheme described by |Mar shall ( 1992 1 for creating upper limits at a given level of confi- 



dence a (we chose to use the "2cr" level of significance, making a = 0.97725) for marginal or negative measurement 
outcomes, under the condition that we have a priori information that the quantity being measured must be positive. 
Here, we derive the formula for evaluating these upper limits and show a plot that allows one to find quickly what the 
upper limits would be at a few levels of significance that differ from the "2(t" value. 



Our derivation differs somewhat from that given by Marshall ( 1992 1, and makes use of an extension of Bayes' 



Theorem for determining the probability of condition A, given conditions B and C, 

Prob(A|B,C) = Prob(A|B)Prob(C|A,B)/Prob(C|B) . (Dl) 

If we designate the true value of the quantity as f , our stated upper limit as m, and have a measurement outcome datum 
d with an associated noise amplitude n, we can make the assignments for the conditions A = t>u, B = t>Q and 



C = d,n. Restating Eq. Dl in terms of our explicit conditions yields 



Prob(f > u\t > 0,d,n) = Proh(t > u\t > 0)Prob(c/,n|f > u,t > 0)/Pvoh(d,n\t > 0) , (D2) 

We set the first term on the right-hand side of this equation equal to 1 (this clearly must be true if m < given that 
t >Q; for M > there is no constraint on the probability). The second term represents the small probability {I -a') that 
we could have obtained an outcome d with noise n when the true value f violates our declared upper limit u. This term 
is given by the area under the one-sided tail of a Gaussian distribution with unit variance and zero mean integrated 
from (u-d)/n to infinity, i.e.. 



u — d 



Pmhid,n\t > M, f > 0) = ierfc ( ] (D3) 



where erfc is the complementary error function, defined as 



Similarly, 



erfc(x)= l-erf(x)= ^ / e\p(-t2)dt (D4) 



Prob(J, n |f > 0) = ierfc ( — £- ) = 1 - ^erfc ( ) (D5) 
\V2nJ \v2n. 



so that our final result for the probability that there will be a violation of the upper limit is given by 



1 - a = erfc 



/ u-d^ 



2 - erfc 



/ d 



\V2n 



(D6) 
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Fig. 23. — : The family of curves for different levels of confidence a, plotted in terms of the declared upper limits di- 
vided by the measurement uncertainties (ordinates), u/n, vs. the measurement outcomes divided by their uncertainties 
d/n (abscissae). 



for which a represents our level of confidence that m is a valid upper limit for t. 

Figure 23 shows the solutions to equation D6 for various confidence levels a. This figure should help one to convert 
our raw measurements of barely or undetected O VI into upper limits of one's choice. The values that we had listed in 
Table |5]corresponded to a = 0.97725, i.e., a "2cr" upper hmit. 



E. Final Spectra 



The spectra obtained from our FUSE survey are shown in Figure 24 These include all the data we deemed usable 
from the PI Team programs P102 and P122 as well as sight lines taken from the FUSE Archive prior to April 2003 
(see 52.1 1. Several stars from programs P102 and P122 were not used, and are not shown here (see 52.10 and Fig. [5] 



instead). For each spectrum in Figure|24j there are three panels which depict the data analyzed in the ways discussed 
throughout ^ Below, we describe the information presented in each panel. 

Top panel: This panel shows the coadded spectra for each sight line. All data come from the LiFl A FUSE channel, 
unless indicated at the bottom right of the panel, where the channel used is indicated. Spectra are plotted in specific 
flux units, along with the continuum used to normalize the data (solid line). The ±lcr error 'envelopes' to the fitted 
continuum are shown as dotted lines about the adopted best-fit continuum. The star's name. Galactic co-ordinates [Ub], 
spectral type, and ROSAT designation R (see ^3.4|and|3.5|l are shown at the top left of the panel. The distance to the 



star d, along with the difference between the distance and the upper and lower bounds to the distance (Appendix B.3 1, 
are also given at the top left of the panel. The FUSE program identification number is given at the top right. A unique 
identification number is shown above the top left corner of the top panel. This number corresponds to the ID numbers 
given in Tables [l] |2] |4) and|5] 

If the wavelength scale has been corrected using data from HST, then an "S" is marked beneath the program number 
Similarly, an "X" indicates that no shift has been made (for the reasons given in Appendix[A|, and the velocities of the 
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absorption should be considered uncertain. Otherwise, a +10 km s"' shift (+0.0344 A) has been applied to the data, 
which is indicated by an "F". (This is the same nomenclature used in column 14 of Table|4]) Stars that lie against the 
X-ray emission of the inner regions of the Carina Nebula or the Vela SNR are noted as such with a label at the top of 
the panel. Note that these are included only to indicate the region of interest, and do not imply that the star is at the 
same distance as Carina or Vela. Sight lines which likely pass through the SNRs listed in Table |8]have a label 'SNR' 
(except for Vela). The flux is shown with the HD 6-0 R(0) line already removed; the region shaded in gray near the 
bottom of the O VI profiles represents the difference in the spectrum before and after the removal of the HD line. 

Middle panel: The data normalized by the best-fit continuum are plotted in wavelength (bottom x-axis) and ob- 
served heliocentric velocity for the O VI A1032 transition (top x-axis). The wavelengths of five lines are marked as 
follows. "1": the predicted position of the H2 7-0 R(6) line at 1030.08 A, based on the measured velocity of the H2 
6-0 R(4) line; "2": the measured position of the H2 6-0 P(3) line; "3": the predicted wavelength of the CI 1 A1031 
line, based on the velocity of the H2 6-0 R(4) line. "4": the measured position of the HD 6-0 R(0) line which is 
removed from the data as discussed in § |2.7[ "5": the measured position of the H2 6-0 R(4) line. 

Line profiles corresponding to the best-fit solutions discussed in ^2.8.1 over-plot the data (histogram) with curved 
lines. If more than one component is fitted to the O VI line, individual profiles are shown with dotted lines. The 
best-fit values of b, v and for the O VI feature are listed at the bottom left of the panel, in units of km s"', km s"', 
and cm"^, respectively. The total column density is given at the bottom right of the panel, along with two 'errors': the 
first two values in parentheses represent the range in A^(0 VI) from the upper and lower continuum fits, N„-N and 
N-Ni, where is the O VI column density measured from the best-fit continuum; the second value in parentheses is 
the error from Poisson statistics alone. The central velocities/wavelengths of O VI absorption components are marked 
with filled triangles. The b and values which are used to model the H2 lines shown, are fixed by fits to other H2 lines 
in different regions of the FUSE data. Two components are sometimes fitted to the H2 lines in these other regions; 
these are indicated with two tick marks for an individual H2 line. 

Bottom panel: The variation of A^(v) as a function of heliocentric velocity for the O VI A 1032 line derived from 
the AOD method is shown in the bottom panel. All A^(v) profiles are plotted at the same scale except for HDs 167402 
(star #1), 168941 (#2), 178487 (#14), and HDE 225757 (#20); flie O VI absorption towards fliese stars is so strong that 
the y-axis is compressed to accommodate the larger profile. 

To enable calculation of an AOD column density even when the O VI line is contaminated by one or both of the H2 
lines, the P(3) and R(4) lines have been removed from the profile (see 5 2.8.2 for more details). The dotted lines show 
the positions and strengths of the removed H2 lines, scaled as if they had the same /-value as O VI A 1032. They are 
included to illustrate how the O VI AOD profile was affected by the H2 lines before they were removed. 

A^(v) profiles derived using the continuum-fitting error envelopes (yielding values of Nu and Ni) are shown above 
and below the best-fit A^(v) (which forms the middle, shaded histogram). The la Poisson noise errors in A^(v) are 
shown as normal vertical error bars for each pixel, although they are only plotted for the middle (adopted) A^(v) profile. 

The total column density measured over all velocities using the AOD method is quoted at the top right of the bottom 
panel. The velocity range used was -120 km s"' < V0 < +120 km s~', although this was modified for profiles that 
clearly extended beyond this interval. The velocity range used can be inferred from the total range covered by gray 
pixels; these are the A^(v) values that contribute to the calculation of the total column density. As with the total column 
densities that we listed for the profile fits (middle panel), the differences in the column density from the upper and the 
best continuum fit, and the best and the lower continuum fit, [A^„ -A^, N-Ni], are shown in the first set of parentheses. 
The value in the second set of parentheses gives the total error from Poisson noise. 

Below the A^(v) profile, a gray bar indicates the range of velocities expected from differential Galactic rotation 
for absorbing gas along the line of sight. (These are the same bars that are shown in Fig. [T7]) Darker colors mark 
velocities where the velocity of the gas is changing only very slowly with distance. This is discussed more fully in 
§3.11 A small '0' above the bar marks the value of given in Table [T] to indicate which end of the bar corresponds 
to distances ofd = Q. As in the middle panel, small triangles along the baseline show the central velocities of the O VI 
components derived from profile fitting (for comparison with the AOD profile), along with the velocities of the H2 6-0 
R(4) components, which are shown as crosses. 
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Fig. 24. — : Spectra of Milky Way disk stars observed with FUSE around the wavelength region of O VI A1032. See 
text in Appendix E for full details. 



The remaining 144 panels showing all the spectra in this 
survey can be found in the full-resolution preprint, 

available at: 

http://www. astro.princeton.edu/^dvb/o6 
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